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I. Cosmic-ray evidence shows that the soft component, 
which persists even under great thicknesses of matter, 
consists, primarily at least, of electrons and y-rays. The 
observations of the magnitude and variation with depth 
of the soft component and showers, and of the size and 
material dependence of bursts, lead to the following rough 
phenomenological description: In addition to ordinary 
ionization losses, the mesotrons have an appreciable chance 
of transferring a considerable fraction of their energy to 
the soft component. For transfers above 2 X 10" v this prob- 
ability is roughly independent of mesotron energy. Under 
10” v the probability of large transfers is much greater, 
roughly 20 times as great as at higher energies. II. The 
production of secondaries of energy <10" v can be ac- 
counted for by a familiar process: the Lorentz contracted 
Coulomb field of the mesotron, in sweeping over an atom, 
ejects a high energy electron. This mechanism, however, is 
inadequate to explain the very large energy transfers, 
>10" v, involved in bursts. The theory which has been 
developed to describe the mesotron of spin one associates 
with transitions in which the direction of the mesotron spin 
changes an intrinsic dipole moment; this dipole field favors 


high energy transfers, and, roughly, can account for the 
observed bursts. However, the dipole field also gives 
probabilities of high energy bremsstrahlung, > 10" v, much 
too large to be compatible with the observations. III. The 
cross sections for the production of electron secondaries 
and bremsstrahlung by the intrinsic mesotron dipole field 
have been estimated under the assumption that the coup- 
ling of mesotron and electromagnetic field is small; when 
this is not so the formula cannot be right. Examination of 
the approximation made indicates that the estimate of the 
cross section for electron secondaries should be valid up to 
mesotron energies of ~10" v, and is thus applicable to 
the bursts. The bremsstrahlung formula, on the other hand, 
fails at 10” v; it thus cannot be used for bursts, while for 
energies 10°< E<10" it leads to no contradiction with the 
experimental evidence. The problem of extending the 
formulae above these critical energies probably goes beyond 
the framework of present theory. The evidence indicates 
that it is the largeness of the coupling, and not the oc- 
currence of lengths smaller than the critical h/uwce~2X10-" 
cm, that limits the applicability of the quantum me- 
chanics. 


I. INTRODUCTION: OBSERVATIONS 


HE persistence of the soft component of 

cosmic radiation under great thicknesses of 
matter shows that this is secondary radiation 
produced by the mesotrons of the penetrating 
component. At sea level the cascade radiation 
which comprises the soft component still has 
three sources of roughly equal importance :! 


1 We are grateful to Professor E. Amaldi for informing 
us of unpublished experimental results. For a general dis- 
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degraded primary cascade radiation of high 
initial energy; secondaries produced by the 
mesotrons in matter; and (less certainly) second- 
aries from the disintegration of mesotrons in the 
atmosphere. But after a few meters water 
equivalent below sea level only the second group 
of secondaries is of importance. 

From this point on, the fraction of soft com- 


cussion, see Euler and Heisenberg, Ergeb. d. exakt. 
Naturwiss. 17, 1 (1938). 
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ponent increases with depth.” Of the order of 10 
percent of the ionizing particles are soft at 10 m 
under sea level ; the proportion is about doubled 
at 30 m. The number of showers and of bursts 
also tails off less rapidly with filtration than the 
number of mesotrons. At 75 m below sea level 
the relative number of showers is about four 
times that at 10 m. This soft radiation has the 
typical behavior of cascade radiation of electrons 
and y-rays, and there is no evidence that heavy 
particles or low energy mesotrons play an im- 
portant part in it. 

These results may be simply interpreted: 
filtration increases the mean mesotron energy, 
and for a power law energy distribution this 
increase is linear with depth. Thus the increase in 
the soft component with filtration means that 
the mean energy transfer from mesotrons to soft 
component must increase with mesotron energy. 

Some information on the energy transfer cross 
section in the region of very high energies is 
given by what is known of the variation of 
burst frequency with burst size and with ma- 
terial. A most striking result here* is that the 
distribution of bursts in size, above roughly 


10! vy, follows the same law as the energy 


distribution of the penetrating mesotrons: the 
number of bursts greater than S in size is given 


Nss~S~’, y~1.8, (1) 


by whereas the number of mesotrons of energy 
greater than E is given by the law, obtained from 
the absorption curve on the assumption that most 
of the absorption comes from ionization losses, 


Nys~E-1, 7~1.9. (2) 


This ‘parallelism of the burst frequency and 
absorption curves follows immediately if for high 
energies the probability of a given fractional 
transfer of energy from a mesotron to soft 
radiation is independent of mesotron energy. 
Such a law of energy transfer does, in fact, give a 
mean energy loss just proportional to the meso- 
tron energy. 

Additional evidence on the form of the cross 
section is provided by observations on the rela- 


2 W. M. Nielsen and K. Z. Morgan, Phys. Rev. 54, 245 
1939)° P. Auger and T. Grivet, Rev. Mod. Phys. 11, 232 

2 A. Sittkus, Zeits. f. Physik 112, 626 (1939); M. Schein 
and P. S. Gill, Phys. Rev. 55, 1111 (1939). 
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tive burst frequencies in different materials. As 
is discussed in more detail elsewhere,‘ a cross 
section depending only on fractional energy loss 
gives, for the frequency of bursts in substances 
of different atomic number, a Z dependence 


(3) 


where g(Z) gives the atomic number dependence 
per atom of the large energy transfers, and J(Z) 
is the critical energy of the cascade theory, which 
increases with decreasing Z. It seems probable 
that, by taking into account the effect of scat- 
tering on the low energy cascade radiation, this 
result, if g(Z)=Z, can be reconciled with the 
observed approximate equality of burst fre- 
quencies in different materials. 

These considerations suggest that, in addition 
to the ordinary ionization losses, mesotrons can 
transfer to the soft component an appreciable 
fraction of their total energy, and that the 
probability of this is roughly independent of the 
mesotron energy. It is, however, not possible to 
fix the constant determining this probability to 
give agreement both with the increase in soft 
radiation on filtration and with the absolute 
value for the probability that a high energy 


mesotron make a burst. To see this, let us write 


for the cross section per atom for a mesotron of 
energy E and mass yu to transfer an energy fE to 
the soft component 


do=aox(f)df, (4) 


Then the probability that in a centimeter of 
material a mesotron of energy greater than E 
makes a burst of energy greater than E is just 


t=ooN f K(f (5) 


where N is the number of atoms per cc of the 
material. According to Schein and Gill® the 
probability that a mesotron will make a burst of 
energy greater than 2X10" v is about &=2X10-° 
per cm. 

An independent check on this figure we may 
obtain by the following argument : the maximum 
of the transition curve for bursts in Pb, at ~5 
cm, is about twice the value at greater thick- 


‘ de Oppenheimer, Rev. Mod. Phys. 11, 264 (1939). 
5 Schein and Gill, reference 3. 
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nesses. Thus the bursts produced by cascade 
radiation give a contribution equal to that from 
the mesotrons: decay electrons will not give an 
appreciable contribution. The total number of 
cascade electrons is about 1/20 the number of 
mesotrons. Of these the theory of showers shows 
that (4X107/e)'* have an energy greater than e. 
On the other hand, if > 10°, roughly (2 10°/e)!-* 
mesotrons have an energy greater than e. The 
bursts they make have a “range” of about 5 cm 
in Pb. Thus we must have, for the probability & 
that a mesotron of energy >« make a burst in a 
cm of Pb, 


(1/20)(4107/6)"?, 
t= (1/100)(1/50)"®~0.5 10-5. 


The value gives 


Now this value is not great enough to account for 
the increase in relative soft radiation with 
filtration. To see this we have to compare ‘the 
total energy lost to the soft component per cm by 
(4) with the mean ionization energy loss of the 
secondaries, which gives then directly the number 
of soft particles per mesotron. This ratio is 


1 
R=mB fdf/4y?c? In (1(Z)/ZRh). (7) 


Here m is the electron mass and E the mean 


mesotron energy. This gives, for the increase in 
ratio of soft to hard component on filtration by 
30 m water equivalent, if we use (6) to estimate 
x,° only 0.5 percent, whereas the experimental 
value is more nearly 10 percent. 

Since the mesotron energies and energy trans- 
fers involved in the increase of the soft compo- 
nent are primarily less than 10"° v, and the bursts 
involve energies greater than 2X10" v, we are 
led to the following rough phenomenological 
description: In addition to ordinary ionization 
losses, the mesotrons have an appreciable chance 
of transferring a considerable fraction of their 


We have supposed for a distri- 


bution «(f) which greatly favored the emission of slow 
electrons this would no longer be right. For example 
«(f)~1/f? would increase the fraction of soft component 
of energy >Emin by a factor In (B/Emin). 


energy to the soft component. For transfers 
above 2X10" v this probability is roughly 
independent of mesotron energy.and is given by 
(4) and (6). Under 10" v the probability of large 
transfers is much greater, and in this range it is 
roughly 20 times as great as at higher energies. 


II. ELECTRON SECONDARIES AND 
BREMSSTRAHLUNG 


The elementary processes by which mesotrons 
transfer energy to the soft component may con- 
veniently be classified as nuclear or electromag- 
netic, according to whether the coupling of 
mesotrons with the heavy particles or with the 
electromagnetic field plays the primary part in 
them. Because of the large nuclear coupling and 
high mesotron mass it has usually been assumed 
that the nuclear processes would alone be im- 
portant for high energy transfers. In fact, calcu- 
lations of nuclear effects, such as mesotron 
absorption and scattering by nuclei, calculations 
based on the perturbation theoretic treatment of 


-the coupling as small, lead to cross sections so 


large, for high energies, that they completely 
contradict the high penetrating power of the 
mesotrons. It has been emphasized especially by 
Heisenberg’ that the prediction of these large 
cross sections rests on the essentially erroneous 
treatment of the interactions as small; despite 
several attempts no reliable estimate of them has 
been given, and this problem probably goes 
beyond the framework of present theory. Nor is 
it sure in what way such impacts will give 
appreciable amounts of soft radiation: the emis- 
sion of y-rays from nuclei excited by such impacts 
hardly seems a very plausible mechanism; the 
radiative capture of a mesotron by a nucleus, 
inverse to the mesotron photo-effect, seems more 
probable. 

Under these circumstances we have thought it 
profitable to re-examine the electromagnetic 
effects a little more closely. For the problem of 
energy transfers to the soft component two types 
of collision are of first importance : elastic impacts 
with free electrons, and bremsstrahlung: to this 
latter corresponds, for the inverse problem of the 
creation of mesotrons by the soft component, pair 
production. 


7W. Heisenberg, Zeits. f. Physik 113, 61 (1939). 
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If mesotrons satisfied the Dirac equation, cross 
sections for these two processes could be obtained 
from the familiar calculations for electrons. For 
the cross section per atom that a mesotron of 
energy £ gives a fraction f of its energy to an 
electron, 

af 


= 200(u2c?/mE)df/f?. (8) 


Similarly, for the cross section per atom that a 
y-ray of energy fE be emitted 


4aZoo 2E(1—f) 4 (1-f) 


a=e?/he. 


do= 


do= 


In this result atomic screening has been neg- 
lected, but the Coulomb field has been smoothed 
out for r<Zih/uc; and the Born approximation 
has, of course, been used. 

The cross section (8) gives a logarithmic 
increase of the soft radiation with filtration; the 
fraction of the electrons of energy greater than 
Emin at depth ¢ is 


2me? /I(Z) 

In —— 

In 10(€/10°)? 


36—41n Z 


where ¢ is the ionization energy loss of the 
mesotron in penetrating to the depth ¢. This 
gives, at sea level, R=13 percent, at 30 m below 
sea level, with Z=12, R=25 percent, at 80 m 
R=31 percent, in rough agreement with the 
observed values, although perhaps somewhat low 
at 80 m. From (9) we obtain R~3 X 10-5Z(e/10%), 
a much smaller contribution than (10). 

For showers of energies X10° v <E<10" v, 
(8) gives, instead of (1), the law Nys~S. 
Both the size dependence and absolute number 
of showers predicted by (8) agree well with the 
experimental results.* 

However, (8) has the wrong energy dependence, 
(9) the wrong Z dependence, to account for the 
observations on bursts. But it is interesting to 
note that both (8) and (9) give for bursts 


*See A. C. B. Lovell, Proc. Roy. Soc. 172, 583 (1939). 


for Emin=10’v, (10) 


78 J. R. OPPENHEIMER, H. SNYDER AND R. SERBER 


>2X10! v, a probability in rough agreement 
with observation; thus (9) gives 


xk~3aZ~2 for Pb, 


whereas the experimental value is x~}. Also (8) 
gives, for all Z and bursts >2 X10", a burst 
frequency corresponding to x~}. 

There is, however, valid ground for doubting 
the applicability to the burst processes of (8) and 
(9). For Yukawa’s theory requires an integral 
mesotron spin, and the spin dependence of 
nuclear forces shows. that this spin must be one. 
The theory which has been developed to describe 
these particles? makes the mesotron electro- 
magnetic current density depend on derivatives 
of the mesotron field, and associates, with transi- 
tions in which the direction of the mesotron spin 
changes, current distributions more singular than 
those of a Dirac electron. The Fourier com- 
ponents of current corresponding to a momentum 
hk > uc behave like k!; for the Dirac electron, and 
for mesotron transitions not changing the spin, 
they behave like k-}. For a classical point charge 
they behave like k°, and for a point dipole like k; 
thus one may roughly ascribe these singular 
currents to an intrinsic mesotron dipole moment. 
The Fourier components of current corresponding 
to this dipole moment are larger by a factor k/yc 
than those corresponding to the electric charge. 

These singular currents of course interact very 
strongly with high frequency radiation fields. 
They radically alter the high energy cross 
sections, giving much larger values than (8) and 
(9). At the same time they introduce couplings 
so large that the question of the perturbation 
theoretic estimate of the cross section requires 
re-examination. 

This behavior of the cross sections is well 
illustrated by the work of Laporte® on the 
Coulomb scattering of mesotrons. Treating this 
problem by the Born approximation, one finds 
for scattering without spin change the Rutherford 
result 


do =4Zoo(uct/E)? csct }0d(cos 0), E>yc?. (11) 


8 Yukawa, Sakata, and Taketani, Proc. Phys.-Math. 
Soc. Japan 20, 319 (1938); Yukawa, Sakata, Kobayasi, 
and aketani, Proc. Phys. -Math. Soc. Japan 20, 720 
(1938); N. Kemmer, Proc. Roy. Soc. Al66, 127 (1938); 
Frohlich, Heitler, and Kemmer, Proc. Roy. Soc. A166, 
154 (1938) ; H. J. Bhabha, Proc. Roy. Soc. A166, 501 (1938). 
QO. Laporte, Phys. Rev. 54, 905 (1938). 
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For scattering with spin change 


da = $6d(cos @). (12) 


For elastic impacts with electrons one finds,'° 
again using the Born approximation, an analo- 
gous result: for no spin change we again get (8) ; 
for spin change 


This result is indeed of the form (4), and thus can 
account roughly for the size and material- 
dependence of large bursts. It gives 


«=0.22. 


(13) 


This is of the order of magnitude observed, and 
together with (8) can probably account for the 
greater part of the high energy bursts. 

For the bremsstrahlung, the high probability of 
large angle scattering given by (12), and the 
large coupling with high frequency radiation 
fields, give a cross section increasing rapidly 
with energy. For an unscreened Coulomb field 
the cross section takes the asymptotic form 
o~aZoo(E/yuc*)*; most of these are close impacts 
and correspond to large angles of scattering and 
radiation. It is therefore essential to take into 
account the modification of the Coulomb field 
within the nuclear radius ~Z'h/yc. If we do this 
by eliminating all Fourier components of the 
Coulomb field with k >uc/hZ', and thus ignoring 
nuclear collisions, we get for high E, still using 
the Born approximation, a cross section of the 
form 


(14) 


It is clear that neither the E nor Z dependence 
of this cross section is in agreement with the 
observations on bursts. For E>2 v, (14) 
cannot be right. However, in the range 10° v 
<E<10" v, (14) is in no contradiction to the 
experimental evidence: its contribution to the 
soft component would be only R~1 percent. 

We see, so that (13) can account roughly for 
soft radiation with E>2 v; that (14) can 
surely not be right for E>2X10" v. On what 


10 This result has been independently derived by Corben 


and Massey, Proc. Phil. Camb. Soc. in press. We are 
indebted to Dr. Corben for telling us of his results. 


grounds can we regard (14) as invalid; and do 
these apply, and if so in what range of energies, 
to (13)? 


III. 


Two distinct but related criteria" are involved 
in the question of the validity of both our results 
(13) and (14) for secondaries and bremsstrahlung: 
Formulae for mesotron charge and current 
density have been used for arbitrary small wave- 
lengths; one may question, as has often been 
done, whether they are right for distances 
<h/wc. On the other hand, we have throughout 
treated the coupling of mesotron and electro- 
magnetic field as small, have used, that is, a 
Born approximation: under what conditions is 
this treatment justified? The formulae (13) and 
(14) will surely have no validity when it is not. 

The former condition is not in Lorentz- 
invariant form, nor can it be reformulated to 
apply only to the rest system of the mesotron, 
since in the impact the mesotron momenta 
change. A natural extension of the criterion is 
that only for those impacts in which momentum 
transfer + and energy transfer « satisfy 


r—eé/ep pe (15) 


can charge and current expressions of mesotron 
theory be used. This condition is thus clearly 
equivalent to the rejection of relativistic mesotron 
theory. It is incisive; there is no evidence that it 
is right. 

To apply this condition we note that in the 
rest system of the center of mass of electron and 
mesotron the energy of the electron does not 
change, and the transverse and longitudinal 
momentum changes are of the same order of 
magnitude. Since the transverse momentum 
transfer is invariant, and is of the order (2mE)!, 
we get from (15) 


E> pc?/2m~10" v. 


It is, as we have seen, only above this energy that 
the terms (13) become of importance compared 
to (8). 

In the same way we may apply (15) to the 
bremsstrahlung. If the Coulomb field were not 


1! W. Heisenberg, Zeits. f. Physik 110, 251 (1938). 
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cut off at the nuclear radius, the scattering and 
radiation would be at large angles in the nuclear 
rest system, and (15) would give 


pct. 


For the screened field the scattering and radiation 
is roughly isotropic in a coordinate system S 
where the mesotron momentum is (2uE)!, so 
that wE>y’c?; in either case the bremsstrahlung 
formulae would have no range of validity at all, 
since they are for the case E>yc?, and (15) 
denies that mesotron theory is applicable to such 
a problem. 

For nuclear problems, where the relative 
magnitude of the coupling energy compared to 
the kinetic energy of the mesotrons is measured 
by gE/uc*, with g a constant of order unity, the 
smallness of the interaction energy essentially 
requires (15). For electromagnetic effects the 
relative magnitude of the coupling is a'E/yc’, 
a=e?/hc; and the coupling may remain small 
even when (15) is violated. 

It is clear from (13) and (14) that for high 
enough energies the coupling terms are not small. 
If we use for each problem the coordinate 
systems in which scattering is roughly isotropic, 
in which therefore for secondaries and brems- 
strahlung, respectively, the mesotron momentum 
P is (2mE)' and (2uE)!, the validity of the Born 
approximation at least requires that the cross 

. section be small compared to the summed areas 
of the partial waves involved. For the second- 
aries, with nearly spherical scattering, this means 
that o<(h/P)*, and thus that 


E<pe?/mZc? ~4X10"/Z v. (16) 


For the bremsstrahlung, the half-width, 50, of the 
angular distribution of the scattered mesotron, 
for fixed direction of emission of the y-ray, is of 
the order (60)?~ ; thus par- 
tial waves up to L?~Z!E/yc? are involved; we 
must have o<(AL/P)? and 


v. (17) 


To see whether these conditions, which are 
certainly mecessary for the smallness of the 
coupling, are also sufficient, let us look more 
closely at the coupling energy density and the 
kinetic energy density. This is simplest in the 
coordinate systems S, where one may use wave 


packets of approximate dimensions //P through- 
out. The kinetic energy density is then of the 
order (P/h)*Pc whereas the coupling energy, for 
those transitions involving a change of mesotron 
spin, is ~(P/h)*a!(P/uc)Pc, and the ratio is of 
the order ° 


oP /uc. (18) 
Then we must have 
E<wc?/ma~2X10" v for the validity of (13), 


and 
E<yuc?/a~10" v for the validity of (14). 


These conditions are more incisive than (16,17), 
because the coupling of the mesotron with its 
radiation field becomes large before the coupling 
with the weaker fields of electron and nucleus do. 
Thus in the problem of the bremsstrahlung the 
maximum value of the electric field of the 
nucleus (in coordinate system S) acting on the 
mesotron is /h?, 
whereas the zero point electric field fluctuations 
of wave-length #/P, which induce mesotron 
radiation of this wave-length, are of the order 
(hc)! P?/h?=(he)pE/h?, and thus (E/yucaZ!)! 
larger. 

The derivation of the results (13) and (14) can 


certainly then not be justified for energies . 


greater than 10° v for the bremsstrahlung, 
2X10” v for the secondaries. For to extend these 
results to higher energies would involve not 
treating the coupling of mesotrons and radiation 
as small. It is well known that this program 
leads to divergent results, and that only the 
roughest correspondence-theoretic analogies offer 
any guide to their interpretation. It should be 
emphasized that it would be of no use to treat 
more strictly, say in the problem of brems- 
strahlung, the motion of the mesotron in the 
Coulomb field ; it is not primarily the magnitude 
of this coupling which makes the trouble. 

This discussion therefore suggests that the 
bremsstrahlung formula (14) may well be right up 
to 10'° v; at much higher energies the observa- 
tional material on burst frequency shows that it 
must be wrong. On the other hand (13), for high 
energy secondaries, can be justified up to 2K 10" 
v, and, as we have seen, can explain at least a 
good part of the bursts observed. 
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It may be mentioned that for the inverse 
problem of mesotron production, pair production 
in a “cut-off’’ Coulomb field again gives a cross 
section of the form ¢=aZ!ooE/yuc?, and again a 
limit ~10' v for its validity. At this highest 
energy the cross section per atom of oxygen is 
10-28 cm?, which is about 300 times too small to 
account for the observed high intensity of 


secondary mesotrons of this energy. A cross 
section 300 times larger would, as a matter of 
fact, be inconsistent with the high penetrability 
of mesotrons. The large observed discrepancy 
between mesotron production and absorption 
probabilities has been repeatedly emphasized.” 


12 See for example, L. W. Nordheim, Phys. Rev. 56, 502 
(1939). 
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The Proton-Deuteron Transformation As a Source of Energy in Dense Stars* 


W. A. Wi_pHack 
George Washington University, Washington, D. C. 
(Received November 4, 1938) 


The rates of energy evolution due to the transformation to helium, starting with the reaction 
H+H=D-+et, in hydrogen at densities of 10* to 10® g/cm*, were calculated on the basis of 
complete degeneracy, and the assumption of a crystal-like spacing of the protons. The results 
indicate that any considerable amount of hydrogen in white dwarf stars would lead to much 
higher luminosities than those observed. Thus the low effective molecular weight (1.5) as 
calculated for some of these stars from the accepted white dwarf model, cannot be due to a 
high content of hydrogen. It might be explained as due to very large content (~100 percent) 
of the helium isotope He® but it is very difficult to see how such large amounts of this isotope 
could be present in these stars. It appears that the paradox can be removed only by revision 
of the observational data concerning the white dwarf radii. 


1. INTRODUCTION 


HE importance of the reaction H+H =D-+et 
as a major source of stellar energy has been 
established by Bethe and Critchfield,' at least for 
stars lighter than the sun. In view of prevailing 
theories of stellar composition, which indicate a 
considerable proportion of hydrogen in most 
stars, it is desirable to investigate the rate of 
energy production due to this H—H reaction in 
dense stars. 

Since the zero-point energy of matter at high 
stellar densities (p= 10° or 107 g/cm*) approaches 
that due to stellar temperatures (10°10? °C), it 
is permissible, for a conservative approximation, 
to neglect temperatures in computing the rate of 
combination of protons. The results may then be 
interpreted as indicating a minimum value for 
the rate of energy production for assumed com- 


* The preliminary results of this work were referred to 
by G. Gamow, Phys. Rev. 55, 723 (1939). 
1H. Bethe and C. Critchfield, Phys. Rev. 54, 248 (1938). 


positions, or as indicating a maximum hydrogen 
content for a given rate of energy production. 

We will first consider an electron-proton gas, 
corresponding to a pure hydrogen composition 
and consider later the modifications necessary to 
take account of heavier particles. We assume 
closest cubic packing of the protons. Then, if 27 
is the distance between nearest protons, the 
volume per proton is 


4(v2)r> and 2r=(v2M/pu)', (1) 


where M is the mass of the proton (1.66 X 10-* g) 
and px is the density of protons in g/cm’. 


2. THE POTENTIAL FUNCTION 


In deducing the form of the potential function 
governing the relative motion of two protons, the 
influence of neighboring protons will be con- 
sidered first. The electrical potential at the 
equilibrium points may be taken as zero. Each 
proton, in its equilibrium position, is more or 
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TABLE I. Numerical results for different densities. 


eH 104 105 106 107 108 
rX10% 3.1 1.4 0.66 0.31 0.14 
1 2.1 4.7 10 21 

0.8 0.75 0.7 0.64 0.57 

E,fev) 22 70 230 750 2500 

dz ‘f(x)dx 1.29 1.23 1.14 1.02 0.81 
6X10-* 2X10-® 1.9X10-* 3.3104 

§ 4.4 3.7 3.1 2.7 2.2 

Y 0.013 0.018 0.026 0.034 0.052 

6X10-* 2.1X10- 5.8X10? 1.1108 4.7x108 


less effectively surrounded by a shell of other 
protons. According to theory, the potential would 
be the same at all points within a uniform shell. 
Actually, the 12 neighboring protons do not form 
a complete shell, but to a first approximation, 
there will be no force acting to resist an initial 
displacement from an equilibrium position. 
Imagine two protons to approach each other, 
one at a time, each moving a distance d from its 
equilibrium position. In accordance with the 
foregoing, there will be no change of potential 
due to the motion of the first. The motion of 
the second involves a change in potential since 
its shell is distorted because of the initial dis- 
placement of the first proton. The work done in 
moving the second proton in the field of the first 
would be e?/(27r—2d) —e?/(2r—d) in the absence 
of the shells. To allow for the effect of the shells, 
we subtract from the above expression the work 
which would have been done in moving the 
second proton in the field of the first proton, 
had the first not been displaced. This gives the 
potential due to the protons, as a function of d: 


1 1 1 
= 
—d 2r—d 2r 


V,(d) =e? 2 
2r 


Substituting (r —d)/r= x simplifies the expression 
t 
3) 


The kinetic energies of the electrons are very 
much greater than their interaction with the 
protons. The electron distribution will therefore 
not be influenced greatly by the displacement of 


2G. Gamow and E. Teller, Phys. Rev. 49, 895 (1936). 
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3G. Gamow, Atomic Nuclei and Nuclear Transformations (Oxford, 1937), p. 94. 


the protons and the effect of the electrons on the 
individual protons will be nearly the same as that 
of a uniformly ‘‘smeared out” charge distribu- 
tion, of density e/4v2r*. For a displacement of a 
proton by a distance d, there will be a force of 


4v2r* 


dynes acting to restore it to its equilibrium 
position. 

The potential due to electrons, for a displace- 
ment of two protons by a distance d, is, in- 


tegrating the force expression : 
V.(d) = (v2/3) (e*d?/2r'). (4) 


Since (d/r)?=(1—<x)*, this expression may be 
combined with (2) to give the total potential 


(5) 


x(i+x) 3 


The approximation involved in assuming that 
the proton shells and the electrons are not 
influenced by the motion of an individual proton 
is on the conservative side, that is, the potential 
(5) is slightly too large. This will give values of 
penetration and energy evolution which are 
slightly low, which is consistent with considering 
them as lower limits. 


e? 
Viz =—(1—x) 
(z) 


3. THE PENETRATION FACTOR 


To obtain the probability of two protons 
coalescing, we may determine their zero-point 
energy as an oscillator and use this value of 
energy in combination with the potential func- 
tion (5) to calculate the penetration factor. Then, 
from the probability of the two protons being 
together, the probability of deuteron formation 
may be found from Fermi’s theory of 8-emission, 
as modified by Gamow and Teller,? and Bethe 
and Critchfield| 

The quantum-mechanical formula for the 
probability of penetration of a potential barrier,’ 
when adapted to the notation of Eq. (5), is: 


vir 


(6) 
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where E’, the total energy, is in units of e?/2r. In 
this expression, (E—V)p is the kinetic energy 
outside the barrier ; (E— V); is the kinetic energy 
inside the barrier; the lower limit of integration 
x, is the ratio of the width of the potential well 
to the equilibrium separation ; the upper limit x2 
is that value of x for which E= V,.), correspond- 
ing to the classical limit of approach. 

The proper frequency of vibration in a sym- 
metrical potential of the form (5) is not directly 
calculable. However, the form of the potential 
function (5) is sufficiently well approximated, 
from x=0.6 to x=1 by the parabola 


V =2.3e(1—x)?/2r, (7) 
for which the proper frequency is 
vo = (8) 


The zero-point energies Ey)=}hv for the values 
of r corresponding to the different densities, are 
given on Table 1, as are the values of x2, where 
V=Eo. 

The factor [(E— V)o/(E—V);]* in (6) may be 
evaluated from the value of Ep just obtained. 
The term (E— V), is approximately the depth of 
the potential hole, 10 Mev, or 16X10-° erg, 
while (E—V)p is given by the value of Ep in 
Table I. The value of the integral was deter- 
mined by mechanical integration of the curve 
plotted from computations. Its values, as well as 
the values of G finally obtained for various 
densities are listed in Table I. Each of the 
variable terms in (6) except the integral have 
been determined as function of 7, and therefore, 
of px (by Eq. (1)). To obtain an approximate 
analytical formula for G, the integral may be 


replaced by an expression of the type a(px)?, 


with a and 6 determined from the values in 
Table I for pp=10* and 10%. After making the 
indicated computations, the penetration proba- 
bility in terms of px is obtained: 


G=6.2X10-'pxt exp (9) 


4. ENERGY EVOLUTION 


The probability of reaction of the two protons 


(the probability of 8-emission per second) is 


given by Bethe and Critchfield as 
B=1.2X10-°M? sec.—, (10) 


SOURCE OF ENERGY IN DENSE STARS 


where 
M= | Ppadr| (11) 


with y, and ya, the normalized wave functions 
of protons and deuteron, respectively, and dr the 
space element. Their detailed calculations of y, 
and ya are not directly applicable here. As a 
conservative approximation, we assume the sim- 
plified radial wave functions: 


=([r(1—x2) (12) 


inside the deuteron potential hole, with ya=0 
outside, and integrate over the radius of the 
potential hole 79, which is taken as 3X10-" cm. 
This gives: 


va=ro%, 


B=1.2X10-'Gro/r(1 —x2) (13) 


for protons in an S state, which is the state 
implied by the form of Eqs. (12). 

The maximum angular momenta computed 
from the zero-point energies and the proton 
separations, which are given in units of # as J in 
Table I, indicate that not all pairs of neighboring 
protons can be considered as in S states. Since 
the probability of penetration is much less for 
other states, we neglect it entirely. We take the 
probability of a given pair of protons having 
relative angular momentum less than 4/2 as the 
probability of their being in an S state. This is 
approximately 


Y=(1/2J)?=h?/16r7MEo, (14) 


the values of which are listed in Table I. 

Since there are twelve nearest neighbors for 
each particle in the close packed cubic lattice, the 
probability of reaction will be increased, but only 
by a factor of six, since in summing, each pair 
of neighboring protons must be counted only 
once. Because of the influence of the spin on 
penetration probabilities, only those protons 


TABLE II. Minimum energy evolution for different densities 
and hydrogen concentration. 


Density (G/cm®) 108 5 X<10¢ 107 5 X10" 
eforCu=1, K 14X10-* 0.26 82 1.3X10* 
(erg/gsec.) W 6X10? 1x10® 3x10’ 


efor Co =0.3, K 1.2X10—* 1.3X10-7 4.9X10~* 4.310? 
(erg/gsec.) W8 2x10" 6x10 
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with opposite spins may be expected to react, 
and this reduces the probability by a factor of }. 
In view of the approximations above, these 
opposing factors might be neglected, but we 
retain them for completeness. 

Each proton-proton combination initiates a 
chain of reactions leading probably to the 
formation of a sHe* nucleus with the evolution 
of 4.3 10-5 erg.! However, Alvarez and Cornog* 
have found that :He® is a stable isotope of 
helium. If :He* is the end product, the energy 
evolution will be somewhat smaller. From the 
value given with Eq. (13), the energy production 
per gram (6X10* protons for a composition of 


100 percent hydrogen) is 
= 4.6 X 10"G Yro/r(1 (15) 


There seems to be no reason to expect that the 
reaction mechanism will be much changed by the 
inclusion of other nuclei besides protons. The 
probability of a given proton reacting might be 
slightly reduced by heavier nuclei among its 
nearest neighbors, but this effect would reduce 
the energy evolution only in proportion to the 
number of protons displaced, say by 34 for 
‘reduction of the hydrogen content from 100 
percent to 40 percent. On the other hand, 
because of the stronger repulsion by the heavier 
nuclei, zero-point energies of the protons would 
increase. From this point of view, (15) may be 
expected to give a minimum value for eq inde- 
pendent of total density as long as the hydrogen 
content is at all appreciable, and is to be in- 
terpreted as giving the rate of energy evolution 
per gram of hydrogen. The energy production 
rate per gram of total stellar material would be 


€=4.6X10"G YroCu/r(1 —x2), (16) 


where Cy is the percentage, by weight, of hydro- 
gen, r is given by Eq. (1) with pp=pCu and G 
and x2 correspond to this value of r. The rates 


‘of energy production computed on this formula 


are given in Table I for different values of px. 
Substituting the expressions (1) (9) and (14) for 
the various factors, in terms of pp=pCu, we 
obtain the approximate analytical formula: 


€=4.3X 105 p!Cy5 exp (17) 


* Alvarez and Cornog. Phys. Rev. 56, 379 (1939). 
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5. KoTHARI’s FORMULA 


While this investigation was in progress, a 
note appeared by Kothari,’ on the lower limit of 
energy evolution in stellar matter. He applied the 
Bethe-Critchfield formula to degenerate protons, 
indicating that proper changes had been made to 
take account of the new velocity-distribution 
law. His final formula was 


€9 = 2.4X 10° (18) 
where Cy is the hydrogen concentration, and 


The rates of energy production calculated by 
the two methods are at extreme variance. 
Kothari’s results are given in Table II, with our 
(approximate) corresponding values. 

It may be instructive to consider briefly in a 
semi-quantitative fashion the modification of the 
Bethe-Critchfield treatment necessary to take 
account of the Fermi-Dirac distribution and to 
compare the resulting energy formula with (17) 
and (18). Bethe and Critchfield' developed their 
formula from the general equation for the number 
of combinations per cm* per sec., p’, (their 


Eq. (30)): 
p’ (19) 


where N is the number of protons per cm*, (V) 
the distribution function for the relative velocity 
V, v the scalar velocity, and o the cross section 
for combination of two protons of relative 
velocity V. From this point on, their develop- 
ment included the evaluation of the integral for 
a Boltzmann distribution function. Presumably 
Kothari diverges at this point to obtain an 
evaluation of the integral for the relative veloci- 
ties in a Fermi-Dirac gas, using the same ex- 
pression for ¢ obtained by Bethe and Critchfield 
on the basis of Coulomb potential barriers. 

The distribution function in (19) may be re- 
placed by a delta-function, which is unity for the 
average effective velocity and zero for all other 
values of V. Without deriving the value of the 
average effective velocity here, it will obviously 
be less than the maximum relative velocity, 
which is, from the Fermi-Dirac formula: 


2Vo=2(h/m)(3N/8r)}, (20) 


5D. S. Kothari, Nature 142, 916 (1938). 
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TABLE III. Observed and com rates of energy production 
in white dwarf stars. 


CENTRAL 


STAR DENsITY G/cm® oss. (CALC.; =0.3) 
Sirius B 10° 10 
O, Eridani 5x 105 3x107? 107! 
Van Maanen’s 10’ 1x10-4 104 
A C 70° 8247 5107 105 
so that 
(21) 


where o is the cross section for the relative 
velocity 2Vo. The energy evolution per gram per 
second is given by 


«=4.3X10-5p'/p. (22) 
From Bethe and Critchfield, for V=2V> 
o=(3X10-/4V;2) exp [—1.45X10°/2Vo]. (23) 


Since the velocity appears in the exponent and 
the faster protons are the more numerous in the 
Fermi-Dirac distribution, the average effective 
velocity will not be greatly different from 2V> 
and the right side of (21) may be expected to 
give p’ to the correct order of magnitude. 

Substituting (20) and (23) in (22) with 
N=pCuX6X 10", we have 


€=3.4X 10° exp [ — 4.26 (24) 


as an expression setting an upper limit for the 
energy evolution on the assumption of Fermi- 
Dirac energy values and cross sections based on 
a Coulomb potential. It may be noted that the 
exponent in (24) is precisely the same as in 
Kothari’s equation, (18). The coefficient is 
greater than Kothari’s by a factor of 1400/(pCu)! 
which for pCy=10' is only about 7. 

Comparing (18) and (24) with (17) it is seen 
that the coefficient in (17) differs from Kothari’s 
by a factor of 178 (pCu)-', or 1.78 at pCy=10°. 
In the exponent, Cy enters to the 4 power on the 
basis of the Fermi-energy and the Coulomb 
potential, but only to the 0.22 power in our 
development. Not only does the exponent vary 
more slowly with pCy in (17) than in (18) or 
(24), but it is only about one-third as large 
at pCy=10°. 

The Bethe-Critchfield formula was derived on 
the basis of a free path long in comparison with 
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the effective radius of the protons—i.e., by using 
a Coulomb potential. But (5) shows that the 
potential differs radically from the Coulomb 
form when the nearness of other protons at high 
densities is taken into account. If we equate the 
Coulomb potential energy e?/r = 3.1X10-"(pCy)! 
to the maximum value of the zero-point energy 
for the Fermi-Dirac proton gas, 


Eo=2X10-(pCu)!, 
we obtain 
pCu>3X10° (25) 


as the condition for proton energies to be great 
enough to allow two protons to approach each 
other within their actual separation distances. 
Thus at pCy=10°, for instance, if we assume a 
Coulomb potential between two protons and 
neglect the effect of electrons, we have the 
paradoxical situation that the Fermi energy will 
permit an approach:to only 5X10-'® cm, while 
the average separation distance is 1.2 cm! 
This paradox disappears at high temperatures 
and low densities; the Bethe-Critchfield assump- 
tions are quite valid for the densities and 
temperatures which they considered (p~10*; 
T~10"). However, it indicates strikingly the 
burden placed on the Fermi-Dirac energies by 
using the unmodified potential at p= 10° to 10’. 

It seems also that our procedure of quantizing 
the motion of the protons in the modified poten- 
tial (5) is preferable to the Fermi-Dirac formula 
for point particles which takes no account of 
the effective finite size of the particles due to the 
thickness of the electrostatic potential barrier. 
The effect of the barrier thickness is to decrease 
the effective separation, or free path, and there- 
fore also the wave-length, thus increasing the 
energy and greatly increasing the penetration. 

However, the Fermi energy increases with 
(pCu)!, while the zero-point energy of an oscil- 
lator between potential barriers of the form (5) 
increases roughly as (pCy)** over the range of 
densities considered. Thus it appears that at 
higher densities (pCy~ 10°), the greater increase 
of the Fermi energy will cause the values of « 
computed by (18) to equal those computed by 
(17). As indicated by the inequality (25), it is 
just in this range that the zero-point energies 
become so great that the form of the lower part 
of the potential barrier is of no great consequence. 
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At even greater densities, it is probable that 
Kothari’s formula is to be preferred, but these 
are beyond the range which is of special physical 
interest. 


6. ASTROPHYSICAL CONSEQUENCES 


The significance of the results is apparent 
in the determination of the composition of white 
dwarf stars. Relevant data on a few of these 
stars are given in Table III, based on values 
given by Chandrasekhar.* The rates of energy 
production given in the table were computed on 
the basis of the central densities. It is seen that 
the computed values are much greater than the 
values obtained from observational data. The 
conclusion is that only a very inconsiderable part 
of most of these stars can be composed of 
hydrogen. 

This conclusion lends support to the theories of 
stellar evolution which place the white dwarfs, at 
the end of the evolutionary scale. It also necessi- 
tates a critical review of either the observational 
data on these stars or their theoretical interpre- 
tation. Without considering in detail the equa- 
tions of structure of degenerate stars, it will 
suffice to recall that certain relations (depending 
somewhat on the model chosen) subsist among 
the variables: mass, luminosity, radius, and 
effective molecular weight. For degenerate stellar 
material, the effective molecular weight, u, is the 
mass, in atomic units, per electron. For hydrogen 
Me is unity; for zHe* and heavier light elements, 
2.0; for heavy elements, 2.0 to about 2.5. The 
heavier elements presumably occur only in small 
quantities, so that it seems to have been assumed 
that there are essentially only two different 
values for u., namely, 1.0 for hydrogen, 2.0 for all 
other elements. Thus, when the other variables 
are determined from observational data, the 
value of yu, obtained from the equations of struc- 
ture is used to determine Cy, the relation between 


ue and Cy being approximately 
Cu=2.0—we. (26) 
The values of Cy so obtained for Sirius B, (the 


white dwarf for which the observational data are 
the most reliable) are rather large ; Chandrasekhar 


*S. Chandrasekhar, An Introduction to the Study of 


Stellar Structure (University of Chicago Press), Chap. 11. 
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(reference 6, p. 433) gives Cy=0.52, and 
Eddington’ gives Cy =0.66. 

It seems clear that these values are untenable 
in the light of the preceding calculations; the 
energy production for such concentrations would 
be vastly larger than those corresponding to the 


observational data. On the other hand, it does 


not seem that the observational data on this star 
are likely to be greatly enough in error to over- 
come this difficulty. 

One possible escape from the paradox is sug- 
gested by the work of Alvarez and Cornog® who 
separated a beam of :He’ from ordinary helium, 
using the 60-inch cyclotron. While the abundance 
ratio of this isotope is apparently very small 
compared to He‘, the fact that it is stable under 
terrestrial conditions suggests the possibility that 
it may be stable under stellar conditions, also. 


The reaction 


H+H=D-+et (27) 
is quickly followed by the reaction 
D+H=:He'+y. (28) 


However, the further reaction assumed by Bethe 
and Critchfield : 


eHe*+ sHe*+ (29) 


might proceed slowly enough for fairly large 
amounts of He’ to be present. 

The value of u, for sHe’ is 1.5. Thus even if it is 
assumed that no hydrogen remains, values of y, 
less than 2.0 (down to 1.5) may still be explained 
as due to appropriate amounts of light helium. 
There are serious objections to this possibility 
for it is difficult to see how the star could arrive 
at the white dwarf stage without passing through 
an earlier stage of higher temperatures which 
would have furthered the reaction (29). It ap- 
pears, therefore, that the difficulties can be 
overcome only by the revision of the values of the 
white dwarf radii, as obtained from the obser- 
vational data. 

In conclusion, I am happy to express my great 
indebtedness to Dr. G. Gamow, for suggesting 
the problem, to Dr. E. Teller for valuable sug- 
gestions and careful criticism of the manuscript, 
and to both for stimulating discussions. 


7A. S. Eddington, Monthly Notice 99, 595-606 (1939). 
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Relation Between Actual and Virtual Ionospheric Height 


H. G. Booker AnD S. L. SEATON 
Department of Terrestrial Magnetism, Carnegie Institution of Washington, Washington, D. C. 
(Received August 24, 1939) 


Translation of virtual height of the ionosphere into 
actual height is achieved by fitting a parabolic maximum 
of electron density to observations of variation of virtual 
height with wave frequency. Methods for measuring the 
thickness of a region and the height of maximum electron 
density are described and applied. The simplest measure 
of the actual height of maximum electron density is the 
virtual height at five-sixths of the penetration frequency. 


This measure is found to be remarkably reliable for F 
region at night. During the daytime it is often necessary 
to correct F, observations for presence of E region, and F, 
observations for presence of F; region. A simple technique 
of correction is described and applied. The disadvantages 
of minimum virtual height for quantitative study of di- 
urnal variation of ionospheric height are emphasized. 


I. INTRODUCTION 


HE object of this paper is to describe a 
method for translating virtual height of the 
ionosphere into actual height. The method is 
quite simple, and yet sufficiently accurate for 
many practical’ purposes. It applies primarily to 
the F region, but difficulty arises when this 
region is separated into F; and F» regions. It is 
assumed that the observational material consists 
of variation with frequency of virtual height of 
the ionosphere for radio waves reflected vertically. 
Virtual height is proportional to time delay. 
Observations of this type over practically the 
entire range of wave frequency for which echoes 
are normally received from the F region are made 
automatically every 15 minutes at the ob- 
servatories maintained at Huancayo, Peru, and 
Watheroo, Western Australia, by the Depart- 
ment of Terrestrial Magnetism of the Carnegie 
Institution of Washington. The problem of trans- 
lating virtual height into actual height is one to 
which thought has already been given.'-* One 
method is to fit to the observations a parabolic 
maximum of electron density (see Fig. 1). It is 
this method which will be considered here. 
In Section II we treat an idealized problem. It 
is supposed that: (i) The ionosphere consists of a 
parabolic maximum of electron density ; (ii) prop- 
agation of radio waves is unaffected by presence 
of the earth’s magnetic field; (iii) there is no 


1E. V. Appleton, Proc. Phys. Soc. 41, 43-59 (1928). 
as 37) V. Appleton, Proc. Roy. Soc. A162, 451-479 


3 FH. Murray and J. B. Hoag, Phys. Rev. 51, 333-341 
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Lorentz polarization correction ;':* (iv) obser- 
vations have been made of variation of virtual 
height with wave frequency from zero frequency 
up to the penetration frequency. The problem is 
to find convenient methods for deducing from 
the observations the constants which define the 
parabolic maximum of electron density. The 
extent to which the results of Section II are 
applicable in practice is discussed in Section III, 
and the actual distribution of electron density 
with height is deduced for two sample ionospheric 
records. In Section IV the results of Sections II 
and III are considered in connection with the 


ELECTRON DENSITY,N 


Fic. 1. Parabolic maximum of electron density. 


*C. G. Darwin, Proc. Roy. Soc. Al46, 17-46 (1934). 
5H. G. Booker and L. V. Berkner, Terr. Mag. 43, 427- 
450 (1938). 
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PENETRATION FREQUENCY, fo 


WAVE FREQUENCY, 
0.4% 0.6% 0.8% 


Fic. 2. Relation between virtual height and wave fre- 
quency for parabolic maximum of electron-density. 


problem of routine scaling of a mass of ionospheric 
records. 


II. PARABOLIC MAXIMUM OF ELECTRON DENSITY 


An ionosphere consisting of a parabolic maxi- 
mum of electron density is depicted in Fig. 1. 
The distribution of electron-density N with 
height h is defined by three quantities: (i) The 
maximum electron density, No; (ii) the height ha 


TABLE I. Tabulation of ¢(x) == tog, 
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at which the electron density is maximum; 


(iii) the height A, at which the electron density 
begins to increase from zero. 
The relation between N and h is 


No—N hu—h\? 
hu—h» 


The quantity 


t=hy—hy (2) 


will be called the thickness of the parabolic 
maximum of electron density. 

An ionosphere of the type shown in Fig. 1 has 
a penetration-frequency fo. Only radio waves of 
frequency less than fo are reflected at vertical 


incidence. It is well known [compare reference 2} 


that, on the assumptions mentioned in Section I, 
the variation of virtual height h’ with wave 
frequency f for the parabolic maximum (1) of 
electron density is given by 


Au—hn 2fo 


This variation is depicted in Fig. 2. 

The problem is: Given the variation of h’ 
with f shown in Fig. 2, to determine the corre- 
sponding variation of N with A shown in Fig. 1. 
This is equivalent to determining the three 
quantities No, hy, and h,, or alternatively the 


(3) 


its 


R 


| 
| 88 6. 
q 
| 
q 
i | me 
| | 
| 
| 
qi 
| 
a 
in 
| 
id . 0.00 0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09 
1,000 -—1.000 -—1.000 -—0999 -0.998 -—0.997 —0.996 -0.995 -—0.994 —0.992 
—0.990 -—0.988 —0.986 —0.983 -0.980 -—0.977 -—0.974 -—0971 -0.967 —0.963 
—0.959 -—0.955 -—0.951 —0.946 -—0.941 -0.936 -—0.931 -—0.925 -0.919 —0.913 
q —0.907 —0.901 —0.894 -—0.887 -—0.880 -—0.872 -—Q864 -—0.856 -—0.848 —0.839 
—0.830 —0.821 —0.812 -—0.802 -—0.792 -—0.782 —0.771 —0.760 —0.749 —0.737 z 
| —0.725 —0.713 —0.700 -—0.687 -—0.674 -—0.660 -—0.646 —0.631 -—0.616 —0.600 b 
—0.584 —0.568 -—0.551 —0533 —0.515 -—0496 -—0477 -—0457 -0436 —0.415 
—0.393 0.370 —0.347 -—0322 -—0.297 -—0.270 -—0.243 -—0.214 -0.185 —0.154 
—0.121 0.087 -—0.051 -—0.014 +0.026 +0.068 +0.112 +0.160 +0.211 +0.265 
+0.325 +0390 +0462 +0542 +0634 +0.740 +0.868 +1.029 +41.252 +1.620 E 
+ +1.678 +1.354 +1.170 +1.045 +0.950 +0874 +0812 +0.759 +0.714 
4 +0.674 +0.639 +0.608 +0580 +0.554 +0531 +0510 +0490 +0471 +0.455 th 
+0.439 +0424 +0410 +0.397 +0385 +0.373 +0.362 +0.352 +0.342 +0.333 
+0.324 +0316 +0.307 +0.300 +0.292 +0.285 +0.279 +0.272 +0.266 +0.260 
J +0.254 +0.249 +0.243 40.238 +0.233 +0229 +0224 40.220 +0.215 +0.211 
+0.207 +0.203 +0.199 +0.496 +0.192 +0189 +0186 +0.182 +0.179 +0.176 
+0.173 +0.167 +0165 +0162 +0.159 +0.157 +0155 +0.152 +0.150 
40.148 40.145 +0143 +0141 40.139 +0137 +0135 +0133 40.131 +0.129 
4 40.127 +0.126 +0.124 +0122 +0.120 +0.119 +0118 +0.116 +0114 +0.113 
+0.112 +0.110 +0.109 +0107 +0.106 +0105 +0104 +0.102 +0.101 +0.100 
a 
q 
7 


* IONOSPHERIC HEIGHT 


three quantities No, hy, and +r. N is derived 
directly from the observed penetration frequency 
fo. If No is measured per cc and fo in mc/sec., then 


N=kfe, (4) 


(5) 


where 
k=1.24X10'. 


The problem is therefore reduced to the following : 


Given the variation of h’ with f shown in Fig. 2, 
to deduce hy and h», or hy and r. 

Define a function 
1+x 


—-—1. 


—x 


¢(x) “ log. (6) 


g(x) is plotted in Fig. 3 and tabulated in Table I. 
Using (2) and (6), one may write (3) 


h'=hu+r- o(f/fo). (7) 


Given the variation of h’ with f, we may deduce 
hy and 7 as follows: (i) Measure the values of h’ 
corresponding to any two values of f/f) lying 
between 0 and 1; (ii) from Table I find the 
values of ¢ corresponding to these two values of 
f/fo; (iii) substitute the two pairs of corresponding 
values of h’ and ¢ into (7), and solve for hy and r. 

In practice it is convenient to choose the two 
values of f/fo, not at random, but so that the 
corresponding values of g are numerically simple. 
Values of f/fo corresponding to a series of useful 
values of ¢ are tabulated in Table II and 
indicated in Fig. 2. The values in Table II are 
deduced from Table I. 

Consider the. wave frequency 


f=0.834fo. (8) 


From Table II the corresponding value of ¢ is 
zero. Hence, at the wave frequency (8), (7) 
becomes 

h'=hy. (9) 


Eqs. (8) and (9) state that, at the wave frequency 
equal to 0.834 times the penetration frequency, 
the virtual height is the actual height of maxi- 
mum electron density. Thus, if h’s34 is the virtual 


10 “4s 


Fic. 3. Graph of ¢(x) =1/2(og. -1). 


1—x 


height at the wave frequency equal to 0.834 
times the penetration frequency, then 


(10) 


hu =h's34 


(see Fig. 2). This remarkable result is valid no 
matter what may be the values of the thickness 
7 and the penetration frequency fo. The peculiar 
property of the wave frequency (8) is demon- 
strated in Fig. 4, which shows the effect upon the 
(h’,f) curve of varying r keeping hy constant. We 
may note that 0.834 (or more accurately 0.8336) 
is practically five-sixths. 

Measurement of h’ at any two of-the wave 
frequencies mentioned in Table II gives im- 
mediately the values of both ky and +r. For 
example, let h’e4s and h’s25 be the virtual heights 
at the wave frequencies equal to 0.648 and 0.925 
times the penetration frequency. From Table II 
the values of ¢ corresponding to these wave 
frequencies are —} and +43, respectively. Hence 
(7) gives 


(11) 
(12) 


/2)r, 
h' 925=hu +(1/2)r. 


TABLE II. Useful wave frequencies. 
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Solving (11) and (12) for hy and 7 we obtain 
= (1/2) (h’s25+h’ as), (13) 
=h'525—h' gas. (14) 


The validity of (13) and (14) is also obvious 
from Fig. 2. 

Formulas for hy and 7 involving values of h’ 
at more than two wave frequencies are easily 
deduced. For example, it is obvious from Fig. 2 
that 


hu =(1/5) { (h’s25-+h’ cas) 
+h’ ssa}, (15) 


t= (4/5) (h’s25—h’eas) 
(16) 


III. DistRIBUTION OF ELECTRON DENSITY 
WITH HEIGHT 


Figure 5 shows a typical example of an 
ionospheric record obtained at night by the 
automatic multifrequency technique.’ Between 
sunset and sunrise the ionosphere consists nor- 
mally of a single main ionized region—the F 
region. We shall show that the results of Section. 
II may be applied directly to records such as 
this to quite a high degree of accuracy. To do 
this we must examine the validity of the as- 
sumptions upon which Section II is based. 


Me 
ACTUAL HEIGHT OF 
MAXIMUM ELECTRON DENSITY 

© 

2 
N 
wave FREQUENCY, 
0.2% 04% 06% 08% 


Fic. 4. Effect upon (h’,f) curve of varying 7 keeping 
hu constant. 


*L. V. Berkner, H. W. Wells and S. L. Seaton, Trans. 
Edinburgh Meeting 1936; Internat. Union Geod. Geophys., 
Ass. Terr. Mag. Electr., Bull. No. 10, pp. 340-357 (1937). 


SCALE OF VIRTUAL HEIGHT IN KM 


Fic. 5. Observed variation of virtual height with wave 
frequency, Watheroo magnetic observatory, March 14, 
1939, 20* 15” to 20* 30” local time (120° East Meridian). 


It is assumed that the ionosphere consists of a 
parabolic maximum of electron density. Ac- 
cording to Section II hy and r may be deduced 
from the observed values of h’ at any two wave- 
frequencies. Table III shows the results obtained 
for the record reproduced in Fig. 5 with pairs of 
wave frequencies mentioned in Table II. Heights 
were scaled by use of the second multiple echo 
(2F). The consistency of results obtained with 
the ordinary wave is quite remarkable. As might 
be expected, results were not quite so consistent 
for the extraordinary wave. There is little doubt 
that, for the record reproduced in Fig. 5, the 
distribution of electron density with height is 
substantially parabolic as shown in Fig. 1. This 
distribution is maintained from the maximum 
down to an electron density, say, 20 percent of 
the maximum. However, as might be expected, 
it is not maintained down to zero electron 
density. This is demonstrated by the fact that 
the minimum virtual height of the F region for 
the record reproduced in Fig. 5 is only 225 km, 
whereas h,, is 247 km. 

It is assumed in Section II that propagation of 
radio waves in the ionosphere is unaffected by the 


TABLE III. Consistency of deduced heights. 


RATIO OF FREQUENCIES} ORDINARY WAVE |EXTRAORDINARY WAVE 
Usep TO PENETRATION 
FREQUENCY hy hy hm 


0.648, 0.925 359 247 112) 350 240 110 
0.725, 0.901 358 247 111 | 350 230 120 
0.757, 0.887 358 248 110 355 235 120 
0.834, 0.969 358 246 112/355 220 135 
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Fic. 6. Observed variation of virtual height with wave 
frequency, Watheroo magnetic observatory, December 18, 
1939, 16* 15" to 16* 30” local time (120° East Meridian). 


earth’s magnetic field. At Huancayo, close to 
the geomagnetic equator, this assumption may be 
regarded as fully satisfied for the ordinary wave. 
But at Watheroo, where the magnetic dip is 
about —64°, the earth’s magnetic field affects 
propagation of both magneto-ionic components. 
However, the work of Goubau’ shows that the 
correction required on this account is at least 
twice as great for the extraordinary wave as for 
the ordinary. It follows that, first, it is better to 
apply the methods of Section II to the ordinary 
wave rather than the extraordinary wave and, 
second, the error in so doing should be less than 
the difference in results obtained from the two 
magneto-ionic components. The difference be- 
tween the virtual heights for the two components 
at the wave frequencies equal to 0.834 times 
their respective penetration frequencies is usually 
between five and ten km when the ionosphere is 
substantially a parabolic maximum of electron 
density. We conclude that effect of the earth’s 
magnetic field upon results obtained from the 
ordinary wave is small. Nevertheless it is 
systematic. 

Similar remarks apply to the assumption in 
Section II that there is no Lorentz polarization 
correction. The only important effect of including 
this correction would be to replace (5) by 


k=1.86X10*. (17) 


7G, Goubau, Hoch: tech. u. Elek: akus. 44, 17-23, 138- 
139 (1934). 


We may sum up by saying that, for records of 
the type reproduced in Fig. 5, the technique of 
Section II applied to the ordinary wave gives 
results sufficiently accurate for most practical 
purposes. 

Between sunrise and sunset, records are not 
usually as simple as that reproduced in Fig. 5. 
Fig. 6 shows a record (the trace at the top of this 
record shows variation of virtual height with 
time for a fixed wave frequency of 4.8 mc/sec.) 
more typical of conditions existing during day- 
light. Three main regions are apparent—E, F,, 
and F:. In these circumstances it is clearly 
impossible to represent the ionosphere by a 
single parabolic maximum of electron density. 
It is possible, however, to represent it by a 
triplet of parabolas as shown in Fig. 7. Fig. 7 was 
deduced from the trace corresponding to the 
ordinary wave on the record reproduced in Fig. 6 
in the following way: 

(i) The technique of Section II was first 
applied to E region. 

(ii) The E region was assumed to consist of a 
complete parabolic maximum of electron density 
and observations of virtual height for F; region 
were corrected for presence of E region. This was 
done by subtracting from the virtual height of F; 
region at wave frequency f an amount 2r- o( f/f»), 
where + and fo are, respectively, the thickness and 
penetration frequency deduced for E region, and 
the function ¢ is given by Table I. 

(iii) The technique of Section II was applied to 
F, region corrected as above. 

(iv) The F,; region was assumed to consist of 
half of a parabolic maximum of electron density 


Fic. 7. Distribution of electron density with height 
corresponding to record reproduced in Fig. 6. 
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Fic. 8. Diurnal variation of h’ and h” for F region, Huancayo magnetic observatory, May 31, 1939. 


Fic. 9. Diurnal variation of square of penetration frequency of F region for ordinary wave, Huancayo 
magnetic observatory, May 31, 1939. 


and observations of virtual height for F: region 
were corrected for presence of F; region. This was 
done by subtracting from the virtual height of F»2 
region at wave frequency f an amount r- ¢(f/fo), 
where 7 and fy are, respectively, the thickness and 
penetration frequency deduced for F; region. 
Correction of F, observations for presence of E 
region was found to be negligible. 

(v) The technique of Section II was applied to 
F, region corrected as above. 

It is believed that Fig. 7 is quite a good 
representation of the distribution of electron 
density with height corresponding to the record 
reproduced in Fig. 6. Note the low actual 
heights associated with F; region. 


IV. RouTINE SCALING OF IONOSPHERIC RECORDS 


In mass analysis of ionospheric records it has 
been the practice to make scalings for Z, F;, and 
F, (or F) regions as follows. Records made at 


hourly intervals are scaled for penetration fre- 
quency fo of the ordinary wave and for minimum 
virtual height. It would now be desirable to 
replace scaling of minimum virtual height by 
scaling of the height 4y of maximum electron 
density and the thickness r by the technique of 
Sections II and III. This would provide practi- 
cally a complete record of the distribution of 
electron density with height in the ionosphere up 
to the level of maximum electron density. 
Moreover, it would then be possible to study 
variation of the quantity (2/3)rkf.? which is 
representative of the total number of electrons in 
a region below the level of maximum electron 
density. Such a program may not be practical 
because of the labor involved. In the meantime 
it is desirable that a less ambitious program be 
drawn up. 

Consider the following program. Scale for the 
F, and Fy; (or F) regions: (i) The penetration 
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Fic. 10. Diurnal variation of h’ and h” for F; and F; regions, Watheroo magnetic observatory, December 18, 1938. 
Fic. 11. Diurnal variation of squares of penetration frequencies of F, and F; regions for ordinary wave, 
Watheroo magnetic observatory, December 18, 1938. 


frequency fo for the ordinary wave; (ii) the 
virtual height h’’ for the ordinary wave at the 
wave frequency 0.834f); (iii) the minimum 
virtual height h’. The result of scaling in this 
way the records from Huancayo for a single day 
is shown in Figs. 8 and 9. Separation of F region 
- during daylight into F; and F: regions was not 
sufficiently marked to take account of. At night 
h’ does not differ greatly from h,, (Fig. 1), and 
gives a useful indication of what may be called 
the bottom of F region. During the day, however, 
h’ is appreciably greater than h,, owing to pres- 
ence of E region. h’’ is quite a satisfactory meas- 
ure of the height 4, of maximum electron density 
in F region throughout this particular day. A 
check on this was obtained by comparing two 
values for the intensity of the earth’s magnetic 
field. One is the intensity deduced for height h”’ 
from the intensity at ground level on the as- 
sumption that the strength of the earth’s mag- 
netic field varies inversely as the cube of distance 
from the center of the earth. The other is the 


intensity deduced for the actual level of maxi- 
mum electron density from the penetration 
frequencies of the region for the ordinary and 
extraordinary waves.’ The average of the 96. 
ratios of the second intensity to the first obtained 
for the day represented in Figs. 8 and 9 is 
0.99+0.01. This is to be regarded as a satis- 
factory check. It may be mentioned that it is 
difficult to obtain a similar check at Watheroo 
for the following reason. Because of the incli- 
nation of the earth’s magnetic field at Watheroo, — 
the ordinary wave is deviated in the magnetic 
meridian towards the south pole and the ex- 
traordinary wave towards the equator. Conse- 
quently, the two magneto-ionic components do 
not penetrate the ionosphere at exactly the same 
points. There is a small systematic difference 
between the maximum electron densities at the 
two points of penetration, and this is sufficient 
to ruin the check unless a reliable correction can 
be applied. 


*E. V. Appleton, Nature 133, 793 (1934). 
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Figures 10 and 11 show a more complicated 
state of affairs than that represented in Figs. 8 
and 9. Figs. 10 and 11 refer to a magnetically 
disturbed day at Watheroo. During the night 
hr’ gives a useful indication of the bottom of F 
region, and hp” is quite a satisfactory measure of 
the height of maximum electron density. During 
the day, however, the situation is very different. 
At sunrise F region separates into F,; and F2 
regions. fre’ increases rapidly to enormous 
virtual heights. Between 06’ 15" and 06* 30”, 
hr’ becomes equal to hr2’’. This is because the 
wave frequency at which virtual height of F: 
region is minimum has increased to 0.834 times 
the penetration frequency of F: region. In these 
circumstances hr,” clearly exceeds the actual 
height of maximum electron density in F: region 
due to presence of F; region. Between 06* 30” and 
11* 00" the wave frequency at which virtual 
height of F; region is minimum actually exceeds 
0.834 times the penetration frequency of F: 
region. Consequently even the minimum virtual 
height of F, region exceeds the height of maxi- 
mum electron density. During such an interval it 
is useless to scale hr2’’ unless it is corrected for 
presence of F; region. Between 11* 00" and 
16* 15" the wave frequency at which virtual 
height of F: region is minimum is roughly equal 
to 0.834 times the penetration frequency of F2 
region. Consequently, hr2’ is approximately equal 
to hp2’’, and only one curve is shown. In this 


-interval the common value of Are’ and hr” is 


necessarily greater than the height of maximum 
electron density, but not to the enormous extent 
prevailing between 06* 30" and 11* 00". The 
broken curve shows the actual height of maxi- 
mum electron density in F: region between 
05* 45" and 17* 45" deduced by correcting for 
presence of F; region by the method outlined at 
the end of Section III. 

The height hp,” in Fig. 10, is appreciably 
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greater than hp,’ throughout most of the daytime 
and gives a useful measure of the height of 
maximum electron density in F; region so far as 
this region has individuality (see Fig. 5). Cor- 
rection of hp,’ for presence of E region would 
decrease values by only about ten km. hip’, on 
the other hand, does not give a useful indication 
of the bottom of F; region due to presence of E 
region ; see Fig. 5, which refers to epoch 16* 15” 
in Figs. 10 and 11. 

It is worth emphasizing that minimum virtual 
height is illadapted for quantitative study of 
solar diurnal variation of the height of F region. 
At night minimum virtual height is a measure of 
the height of the bottom of the region. But 
during the day minimum virtual height usually 
measures the height of a level appreciably above 
the bottom of the region, and not uncommonly 
above even the level of maximum electron 
density. This arises partly from variable re- 
tardation in a lower region, but mainly from 
variation in the ratio of the wave frequency at 
which minimum virtual height is measured to the 
penetration frequency of the region under 
examination. 

Such applications as have so far been made of 
the technique described in this paper indicate 
that it should prove a powerful method of 
approaching a number of ionospheric problems. 
Berkner, Wells, and Seaton have already used 
the technique in connection with an intense 
magnetic storm.? In studying regular diurnal 
variation of the ionosphere it is desirable to use 
observations averaged for a number of days. 
Monthly averages will shortly be available. 

Throughout the conduct of this work we have 
had the benefit of helpful discussion with L. V. 
Berkner, to whom we desire to express our 


thanks. 


*L. V. Berkner, H. W. Wells and S. L. Seaton, Terr, 
Mag. 44, 283-311 (1939). 
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For many years it has been generally known that strong, 
sharply defined radio echoes occasionally return to the 
sending station, or to points within the normal skip-zone, 
after having traversed a path which greatly exceeds the 
round trip distance to the F layer of the ionosphere. Previ- 
ous investigators tentatively ascribed such effects to reflec- 
tion from distant mountains or from concentrations of ions 
in the polar regions. Studies of numerous reflection patterns 
strongly indicate that the delayed echoes are returned from 
regions where there is marked curvature of the F layer. 


A region of this sort normally occurs at the edge of the 
sunlit zone and can turn back a ray which may have 
traveled many thousands of kilometers around the dark 
side of the earth. Small nighttime variations in the 
curvature of the F layer are of very common occurrence 
and are believed to explain several phenomena, including 
much long period, long distance fading. These variations 
are frequently of cyclic character and may tentatively be 
attributed to wave motion in the upper atmosphere. 


INTRODUCTION 


ONTINUOUS records of ionospheric layer 

heights frequently exhibit a phenomenon 
which may be described as a sudden variation in 
the number of multiple reflections returned from 
the F layer. A record may be obtained showing, 
say, three multiples which have intensities de- 
creasing with increasing path length. After a 
time the second and third multiples will disappear, 
leaving a record typical of rather ‘“‘weak”’ re- 
flections. This condition may persist for fifteen or 
twenty minutes and be followed by a period of 
extremely “‘strong”’ reflections, when as many as 
eight or ten multiples are observed. Since the 
coefficient of reflection at the F layer, for the 
condition of total internal reflection, is always 
very nearly unity, it is necessary to explain this 
behavior in terms of large variation in the 
absorption of energy in the lower layers or to 
assume that the reflecting surface is not (ap- 
proximately) concentric with the earth. 

We find many records which exhibit little 
variation in the strength of the echoes over 
periods of hours before or after the occurrence of 
this phenomenon. Examination of such records 
leaves an impression that the bursts of multiples 
indicate an intensity far greater than normal, 
while the energy returned in the intervals be- 
tween bursts is below that which we would 
expect. This fact can be explained on the 
assumption that a series of waves appears at 


times on the surface of the F layer. At a moment 
when the surface over the transmitting (and 
receiving) location is concave downwards, the 
transmitted energy is focused upon the receiver ; 
and when the layer, at the point of reflection, is 
concave upwards the received energy is sub- 
normal. It should be noted that, if the radius 
of curvature is considerably greater than the 
equivalent height of reflection, the focusing effect 
will become stronger as the total path length 
increases, thus partially offsetting the inverse 
distance effect. With a very moderate assumption 
of curvature, in fact, we may expect the sixth or 
eighth multiple reflection to be returned as 
strongly as, say, the third. 

It will be noticed that the illustration of the 
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. 1. “Staggered’’ multiple reflections from the F layer, 
recorded at 3.5 megacycles, 
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Fic. 2. Another example of ‘‘staggered”’ multiple reflections. 


phenomenon shown in Fig. 1 exhibits a remark- 
able modification of the effects to be expected 
under the simple theory outlined above. It is 
characteristic, at least in the early evening, that 
the onset of one of these bursts of reflection is 
indicated by the appearance, in addition to the 
first two or three multiples, of two or three 
multiples of high order, such as the eighth and 
ninth. With the passage of a few minutes the 
highest multiples will disappear while those next 
below appear in sequence. Five or ten minutes 
after the onset of the high order reflections, we 
find that the entire sequence of multiple reflec- 
tions has been run through in reverse order, 
leaving only the lowest (continuous) orders. An 
interesting modification of this pattern is shown 
in Fig. 2, where each multiple remains on the 
record until all have appeared, and shortly 
thereafter the abnormal multiples fade out 
simultaneously. 

We can conceive of no disposition of low layer 
absorption which could cause the attenuation of 
certain multiple reflections while permitting 
those of higher order to return to the region near 
the transmitter. We must, therefore, see if our 
hypothesis of varying curvature of the surface of 
the equivalent reflector can explain the phe- 
nomenon. Unfortunately our knowledge of the 
geometry of such second-order curvature of the 
layer is indirect and inaccurate. We hope, how- 
ever, to show that curvature of a type which may 
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reasonably be expected to exist can explain 
records of the sort we have obtained, and to 
suggest possible effects of such curvatures upon 
communication circuits. 


A QUALITATIVE EXPLANATION OF 
THE PHENOMENA 


Figure 3 is a cross section of an equivalent 
reflector upon which we have indicated a number 
of ray paths. This diagram represents a condition 
we believe to exist and is shown in exaggerated 
form in order that we may trace a number of ray 
paths with a minimum of confusion. We have 
constructed this equivalent reflector by super- 
imposing one wave-length of a cosine wave upon 
a plane, the cosine curve having a wave-length 
of 500 km and an amplitude of 5 km. Since this 
surface is large and at a small angle to the 
horizontal, we can always find a ray from the 
transmitter position O (which we have chosen to 
locate underneath a node of the curve) which is 
perpendicular to the reflector. This perpendicular, 
OA, represents the path of a ray which will be 
reflected directly to a receiver which shares the 
transmitter location. We can trace a very large 
number of rays which leave O at different angles 
to the vertical. These rays will be reflected back 
and forth between the earth and the layer and, in 
general, will not be returned to O. A few rays, 
however, will sooner or later strike a reflecting 
surface normally and will then be sent back 
along their outward courses and returned to O 
after having traversed the region between the 
earth and the layer a number of times. Consider, 
for instance, a ray which departs from O along 
the path OC, which happens in this case to be 
nearly vertical. This ray is reflected in sequence 
at Ci, C2, C3, and so on until it meets the earth 
normally at Cs. It is then returned exactly along 
its former path, reaching O where it can be 
recorded (if it be not too weak) as the eighth 
multiple reflection. 

Figure 3 has been constructed so that we may 
use its nomenclature to discuss this problem 
qualitatively. Consider the rays which first 
reflect from the layer to the left of A, the point 
marked by the directly reflected ray. Some ray 
such as OB, will be reflected vertically downward 
to Bz and then returned to O as the second 
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multiple. Another ray, not shown, will, at its 
third reflection, meet the layer normally and 
become the third multiple. This will continue 
until some ray, such as OM,, is lost because it has 
not been turned back towards O before escaping 
from under the curved region. We find, thus, 
that, if we examine in sequence all rays to the 
left of OA, among a multitude of rays which are 
lost to us are a few which are returned in order as 
the first, second, third multiples, and so on. This 
is exactly what we would find in the case of a 
plane reflector, except that in our case the 
intensity of the multiple reflections is decreased 
because of their repeated reflection from a convex 
surface. 

For rays reflecting originally to the right of A 
the situation is changed in several respects. We 
have found above that the ray OC, returns 
finally as the eighth multiple. Similar tracing of 
the rays OD;, OE,, and OF, shows that they 
return, respectively, as the seventh, eighth, and 
sixth multiples. Other rays not shown are, of 
course, also returned. There is a limiting angle in 
this case, as in the other. Rays which leave O 
at too great an angle to the vertical will escape at 
the right-hand side of the diagram before being 
turned back. Such a ray is shown as ON,. In 
general, the farther to the right the original ray is 
reflected, within the limit, the sooner it is likely 
to return because it cannot turn back until it has 
reached a negative slope to the right of O. At the 
same time the probability of its being returned to 
O is decreasing as its original angle to the vertical 
increases. Two factors operate to increase the 
intensities of multiple reflections returned to O. 
In general, these rays are being reflected from a 
concave surface and are consequently focused 
back, to some extent, upon O. The other reason is 
that a number of rays may be returned to O after 
having made the same number of traverses of the 
region between the earth and the layer; a 
condition which has no parallel in the case of a 
plane reflector. For example, the ray paths 
labeled C and F, in Fig. 3, each contribute to the 
8th-order multiple reflection. 

We may summarize this case by stating that 
the normal set of multiple reflections which we 
would expect from a plane reflector are now 
represented by the reflections from the left-hand 
part of the diagram, these being weaker than in 


the plane case, especially in the higher multiples. 
In addition to this set of reflections we have an 
entirely new family reflected from the right-hand 
side of Fig. 3. In this family the higher multiples 
are more probable and there are cases of several 
rays combining to form a single, intense reflection 
of a high order, as well as a general increase in 
intensity due to a focusing effect. There is a limit 
set by the geometry of the layer which operates 
in such a way that no multiples of less than a 
certain order can be present in this family. If the 
reflections from the layer shown in Fig. 3 were to 
be recorded with an instrument having a finite 
threshold sensitivity we might well expect to see 
a record showing the presence of the first, sec- 
ond, third, sixth, seventh, and eighth multiples 
only. 

It is obvious that the sixth multiple (made up 
of OF,, FiFs2, etc.) in this case would have a 
total path length much greater than six times 
OA, while in the actual case path lengths for the 
higher orders do not differ by more than one or 
two percent from multiples of the layer height 
indicated by the first reflection. This means that 
the layer we have sketched is indeed an exagger- 
ated example of any similar formation which 
may appear in the F region. The same conclusion 
follows from the fact that if the layer height 
varied by ten kilometers from point to point and 
minute to minute direct evidence of the height 
variations would be obvious. We must, therefore, 
restrict ourselves to postulation of the existence 
of ripples in the F layer of not more than one or 
two kilometers amplitude and consider that the 
phenomena of Fig. 3 are compressed into a 
lateral space of not more than fifty or a hundred, 
rather than five hundred, kilometers. Even with 
this restriction, however, it seems reasonable to 


Fic. 3. A cross section of an assumed equivalent reflector 
showing a number of possible ray paths. 
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assume that this simple theory of a “dimple” in 
the layer can explain not only the observed 
variations in the number of multiple reflections, 
but also the fact that, from time to time, certain 
intermediate multiples are missing. 

If we consider that a series of waves of this 


Fic. 4. A straight line  “< e to the reflector of 
ig. 3. 


form roll across the surface of the F layer it is 
easy to see certain variations which are to be 
expected in the reflection pattern. When the 
transmitter is under the center of the part of the 
curve which is convex downwards, all the multiple 
reflections will be present in their normal order 
but with intensities which decrease more rapidly 
than the increase in the squares of their path 
lengths. When the layer is symmetrically concave 
towards the transmitter, all the multiples will be 
present but with increased intensities. At inter- 
mediate positions certain multiples may be 
returned with great intensities and other multi- 
ples of lower order may be forbidden or present 


only very weakly. 


REFLECTIONS From A LAYER OF SYMMETRICAL 
TRIANGULAR FORM 


- This problem does not lend itself to mathe- 
matical treatment, and the complexities of a 
graphical solution for all possible layer curvatures 
and transmitter positions are insuperable. We 
can, however, reduce our dimple in the layer to 
the simple triangular form shown in Fig. 4 and 
draw some useful conclusions with the aid of the 
nomenclature we have indicated. We have given 
a region within which the layer has a slope of 
tan @ and, to the right, a similar zone where the 
slope is tan(—a). From the transmitter at O, on 
the left-hand side of the diagram, a ray departs 
at an angle @ to the vertical. So long as this ray is 
being reflected between the part of the layer 
having a positive slope and the earth, its angle of 
incidence upon a reflector increases by a at each 
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reflection. Similarly, when the ray is in the right- 
hand part of the diagram, the angle of incidence 
decreases by a at each reflection. Since the angle 
of incidence upon either the layer or the earth 
must be zero in order to have the ray return 
along its outward path it follows that @ must 
equal ma before a ray can be of interest to us, n 
being a positive or negative integer, or zero. This 
greatly limits the number of rays we must 
consider for any given extent and slope of the 
layer and any transmitter position, but it is still 
uncertain whether each of these possible rays will 
return to O, or, if so, after how many reflections. 

If a is so small that we may take tan na=na, 
it is easy to express the distance OS: as 


se= 2h(0+a), (1) 


where h is the average height above ground of the 
points of reflection on the layer. This may be 
extended to: 


=n(6/a-+n/2)ha, (2) 


where m is the number of times the ray has 
traversed the space between the earth and the 
layer, and @ and a are to be taken as positive for 
the departing ray shown in Fig. 4. This equation 
breaks down at the first reflection from a slope 
other than a, but with its aid a tabulation may 
be made showing the horizontal distance trav- 
ersed by each of a family of rays leaving O at 
angles to the vertical of 0, +a, +2a, etc. The 
distances from one point of reflection at the 
earth to the next of course decreases, after the 
rays pass under the region BA’, in the same way 
that they had increased before. 

It is easy, after these preparations, to trace 
each possible ray for a given set of values of h, a, 
the distances AB and BA’, and for each of a 
number of transmitter positions. We find that, as 
O moves from A toward B, the first returning 
rays (as we examine those with @ positive in 
sequence) appear as high multiples. As O ap- 
proaches B’, the energy returning to O is con- 
tained in lower and lower multiples until, when 
O is at B’, all multiples are present as though 
reflection were taking place from a plane surface. 
As O moves from B’ toward the right, the 
procedure repeats itself in the reverse direction. 

We have traced the rays corresponding to the 
case of a dimple 2 km deep and 200 km in total 
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length in a layer at a height of 250 km. The 
results are exhibited in Fig. 5 where the solid 
lines represent reflections of rays which departed 
with @ positive, and the dotted lines the reflec- 
tions for @ negative. Widened lines indicate cases 
where two or more rays return to make up a 
given multiple reflection. 

A number of complex rays return as even 
higher multiples than those shown. These are the 
ones which depart with @ negative, reflect first at 
the left of O, then at the right, and finally return 
to O after very long delays. They are omitted 
from the figure because their long path length 
probably reduces their intensity to a value which 
we could not detect. 


THE EFFECTS OF ASYMMETRY 


Figure 5 represents the reflection pattern we 
would expect to record if a triangular wave of 
this sort rolled over the transmitter. It differs 
from the most usual type of record chiefly in its 
symmetry. Occasionally, in the middle of the 
night, a pattern of this sort is observed, one being 
shown in rudimentary form at the left of Fig. 6. 
In general, however, the phenomenon is charac- 


9 
8 
4 7 
5 6 
= 
5 
------ 
: 
z 
| 4--------------- - 
-100 ° +100 


DISTANCE IN KILOMETERS 


Fic. 5. Acomputed reflection pattern. The ordinate may 
be expressed in time to make this figure comparable with 
the experimental records. 


teristic of the three or four hours after sunset, and 
at this time the right-hand part of Fig. 5 is 
absent. The most reasonable explanation of this 
asymmetry is that waves are superimposed upon 
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Fic. 6. An example of the variation of multiple reflections 
with time. The portion of the record centered at 2300 
should be compared with Fig. 5. 


a sloping layer. Examination of a number of 
records gives us an average increase in layer 
height (from west towards the east) of about one 
kilometer in a hundred, during the evening. If we 
superimpose the ripple of Fig. 4 upon this slope 
we find that the slope AB is only one-third of 
that between B and A’. The region BA’ is now 
very effective in turning back rays which origi- 
nated at the left of B, but rays which are 
traveling towards the left will, in general, escape. 
This asymmetry is shown in Fig. 7, which was 
drawn up under the same conditions as Fig. 5, 
except for the slope of —0.01 upon which the 
ripple was superimposed. The seventh, eighth, 
and ninth multiples appearing just to the right of 
O are part of a sort of second-order pattern most 
of which is contained in multiples higher than 
those shown. This part of the pattern is sensitive 
to small changes in the geometry of the layer and 
may be expected to be rudimentary in many 
cases. It is probable that, when it is present, it is 
blurred sufficiently to merge it with the primary 
pattern, the combination giving a record of the 
type of Fig. 2. 


THE MAGNITUDE OF RIPPLES IN THE F LAYER 


We have remarked that the amplitudes of 
ripples of this sort cannot be more than one or 
two kilometers, or the height variations would be 
directly measurable. Since this is the case, we 
expect that the phenomenon of the forbidden 
multiple reflections should be observed only when 
the receiver and transmitter are very close 
together. This is because some of the transmitted 
rays must pass laterally beyond the receiver and 
return to it, and this is impossible unless @ and a, 
in Fig. 4, are of the same order of magnitude. For 
transmission over distances comparable to the 
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height of reflection, 6 must be many times larger 
than any angle between the earth and the ionized 
layer could possibly be. This conclusion has been 
checked experimentally. The records reproduced 
in Figs. 1, 2, and 6 were made with the receiver 
within 250 meters of the transmitter, while 
simultaneous records were made at a distance of 
240 kilometers and showed no indication of the 
phenomenon. As a check on the maximum slopes 
in the layer, the receiving equipment was re- 
moved to a location about four kilometers from 
the transmitter and operated there for about ten 
days. The records made at this distance exhibited 
the effects of focusing in the variation of the 
number of multiples, but the echelon formation 
of multiple reflections was present only in 
rudimentary form in one or two cases although 
the records made at 250 meters distance im- 
mediately before and after the removal indicated 
a number of occurrences nearly every day. 

If we were to draw a circle through the points 
A, B, and A’ in Fig. 4 we would find a radius of 
curvature of 2500 km, and this figure may well be 
increased to 4000 km to allow for the defocusing 
effect of the convex surface of the earth. It is 
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Fic. 7. The computed reflection pattern corresponding to 
the layer of Fig. 4 and the records of Figs. 1 and 2. 


easy to calculate the variations in intensity of a 
ray reflected between a cylindrical mirror and a 
plane reflector, as the number of multiple 
reflections increases. The results of such a calcu- 
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Fic. 8. The computed variation in intensity of a ray 
reflecting between a plane and a cylindrical mirror sepa- 
rated by 300 km. 


lation, for a layer height of 300 km and a radius 
of curvature of 4000 km, are shown in Fig. 8. 
Perfect reflection and no attenuation in the 
absorbing layers have been assumed, but the 
effect of attenuation may be simulated by giving 
the ordinates in the figure an arbitrary negative 
slope. Since a truly cylindrical curvature cannot 
be expected, it is probable that the cusps com- 
puted for the concave case are blurred out of 
existence so that the values fall more or less along 
a smooth curve which, for multiples beyond the 
second or third, varies inversely as the multiple 
number rather than as the square of that number. 
The most interesting conclusion from the figure 
is that we may expect the high multiples, such as 
the eighth or tenth, to vary as much as 20 db if 
the reflecting surface becomes alternately concave 
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Fic. 9. A record showing the presence of a fourth 
multiple reflection which is stronger than either the second 


or the third. 
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Fic. 10. The computed variation in relative intensity of 
several multiple reflections for a case where the curvature 
is large. 


and convex. An effect of this magnitude would be 
caused by the passage of a sine wave, such as 
that shown in Fig. 3, having a wave-length of 200 
km and an amplitude of 0.25 km. 

Occasionally .random curvatures of greater 
magnitude occur. The record corresponding to 
one of these is reproduced in Fig. 9 and, although 
made with the receiver sensitivity below normal, 
shows that the fourth multiple was, for a time, 
nearly as strong as the first, and that the second 
is absent while the third is represented by a very 
weak trace. The variations in intensity compare 
favorably with those computed for a cylindrical 
curvature having a radius of 1630 km which are 
exhibited in Fig. 10. 


THE FADING OF LONG DISTANCE SIGNALS 


It is obvious that the parts of this discussion 
dealing with phenomena which are detectable 
only in the special case of vertical incidence 
observations can have little importance in practi- 
cal communication. They have, however, enabled 
us to postulate the frequent occurrence of sig- 
nificant changes in the curvature of the F layer. 
For communication over short distances these 
changes are unimportant, but when the total 
length of the transmission path approaches the 
radius of curvature they become very interesting. 
It is possible that much of the long period fading 
observed on long distance circuits may be a 
result of this effect. This type of fading occurs 
too slowly to be a result of polarization changes 
or interference between a number of rays. It 


seems to us that the hypothesis of varying 
attenuation in the absorbing strata may not 
explain the cases of slow fading observed at night 
and at high frequencies, while the mechanism we 
are suggesting is not, within limits, a function of 
frequency. An example of what may well be a 
focusing effect is shown in Fig. 11, which is a 
record of the field strength of a short wave 
broadcasting station operating on a frequency of 
about 12 megacycles at a distance from Cam- 
bridge of 5200 kilometers. The transmission 
begins at 0115 E.S.T., when the ionization 
density has fallen so low that only one ray path is 
effective. This is indicated by the absence of the 
rapid fading caused by the interference between 
several rays. The signal remains steady and 
pleasant to hear but with a smooth fading cycle 
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Fic. 11. A record of the field strength of a transatlantic 
short wave broadcasting station. Note that the time scale 
reads from right to left. 
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having a period of about ten minutes and a full 
amplitude of nearly 20 db. This lasts until about 
0155 E.S.T. when the ionization density falls 
below the limiting value, the signal strength 
decreasing by 30 db in the next few minutes. It 
appears to us that this slow, cyclic fading without 
change in the general signal level is evidence of 
the focusing effect we have postulated. 


An HyporuHEsis OF RIPPLE FORMATION 


Under the conditions obtaining in the F region 
in the evening of a magnetically quiet day, the 
equivalent height of reflection increases by about 
twelve or fourteen kilometers per megacycle in 
the region near three or four megacycles. This 
can be interpreted to mean that the density of 
ionization increases by about five percent for one 
kilometer increase in height. If, then, irregularities 
in height of this order of magnitude are to be 
explained by assuming that some non-uniform 
source of ionizing energy is operative, the local 


3s 
300 

| 


EASTERN STANDARD TIME 


Fic. 12. The variation of F layer height with time in the 
sunrise period (3.5 megacycles; November, 1938). 


variations in intensity of this source must be of 
the order of ten percent of the intensity which 
would maintain the ionization in equilibrium at a 
value of 10° electrons/cm*. While such variations 
in incident energy are not impossible and are, in 
fact, smaller than variations which are known to 
affect the E region, it is improbable that the 
effect we are discussing is due to this cause. If the 
total ionization in part of the F region were 
changing in such a way as to cause variations in 
the equivalent height of reflection at a fixed 
frequency, the changes would give rise to phe- 
nomena of other types, such as cyclic variations 
in the critical frequency, which have not been 
reported. 


We prefer to believe, therefore, that these 
effects are due to a temporary redistribution of 
the free electrons already existing in the lower 
part of the F region. The simplest mechanism 
which will cause such a redistribution is an 
atmospheric compressional wave. We have noted 
that our phenomenon occurs generally in the 
period immediately following sunset. This per- 
mits us to postulate a disturbance, set up by the 
cooling of the upper atmosphere at sundown, 
which is propagated in the direction opposite to 
that of the advancing shadow. It will be noted 
that this direction of propagation, together with 
the slope in the equivalent layer, is necessary to 
explain the observed asymmetry in our records. 

It is difficult to estimate the velocity of propa- 
gation in the tenuous medium which exists at the 
height of the F region, where the mean free path 
is probably several kilometers and the tempera- 
ture and pressure can only be guessed at, but it is 
obvious that waves can be set up whose length is 
large compared to the mean free path. Such 
waves would cause compressions and rarefactions 
of the free electrons in the F region in essentially 
the same ratio as in the atoms or molecules of gas 
if, as is unquestionably the case, the period of 
half-recombination is long compared to the 
period of the wave, and these changes in electron 
density would appear as changes in the height of 
the equivalent reflector. 

We are not able to deduce the wave-lengths of 
these disturbances accurately, but they seem to 
be of the order of two or three hundred kilometers. 
The records which exhibit a definite cyclic 
character give a mean period of about 20 minutes, 
although it often varies between half and twice 
this figure. The velocity corresponding to these 
estimates is some six or eight hundred kilometers 
per hour or, say, 200 meters/second. This figure 
is of the magnitude of the velocity of sound and 
is, therefore, somewhat interesting. It is barely 
possible that a careful study of this phenomenon, 
involving simultaneous observations at two or 
more locations, might lead to a_ reasonably 
accurate figure for the velocity of sound at the 
300-kilometer level. 


FocusING OF ECHOES AT SUNRISE 


In contemplating the effects of changes in 
curvature of the equivalent F layer reflector it is 
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necessary to consider the curvatures associated 
with sunrise. This region is important because 
there is a large change in effective height in a 
horizontal distance of a few hundred kilometers, 
with correspondingly large curvatures, both con- 
vex and concave. Fig. 12 shows the average 
variation of equivalent height with time for a 
number of magnetically quiet days in November, 
1938. The dotted line indicates the height at 
which the sun’s rays strike the region after 
grazing an “‘opaque”’ sphere of 6600 km radius, 
and is shown here only because it is interesting to 
note that the layer height begins to fall before 
there is any directly incident ultraviolet radia- 
tion, indicating that diffusion is an important 
factor at these heights. 

If we assume that there is no variation in 
electronic density with longitude (as distin- 
guished from time) the heights of Fig. 12 may be 
plotted as a fynction of distance by use of the 
fact that sunrise, in these latitudes, advances at 
some 1200 km hour. On this sort of graph radii of 
curvature may be measured which vary, in this 
average case, from infinity to about 5000 kilo- 
meters. If the curvature is somewhat greater than 
average we should immediately expect an in- 
crease in intensity of the echoes returned from 
the concave region exactly as in the case of the 
random curvatures discussed above. Under 
magnetically quiet conditions and when the fre- 
quency of observation is well below the critical 
frequency, this deduction is confirmed experi- 
mentally. An example of this sort of record is 
exhibited in Fig. 13, where the focusing of the 
fourth and fifth multiple reflections is clearly 
marked. 


THE RETURN OF LoNG DELAY ECHOES 


At this point we should emphasize the distinc- 
tion between the sort of concavity in a layer 
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Fic. 13. A record showing the focusing of the fourth and 
fifth multiple reflections at sunrise (November, 1938). 


which can cause the focusing of a beam of 
incident energy, and the effect of slope in a layer 
which can turn back a ray which has had a 
component of its path in the direction of the 
lower equivalent height. The latter effect is of 
especial importance in the sunrise zone. A ray 
which has been propagated along the dark side 
of the earth will, upon entering the sunrise zone, 
meet the reflecting layer and the earth at an 
angle of incidence which decreases with each 
reflection until it reaches zero, or until the ray 
passes on into a region in the sunlit zone where 
the layer is again concentric with the earth. If the 
angle of incidence is reduced to zero and if the 
horizontal component of the direction of the ray 
is perpendicular to the edge of the shadow, the 
ray will be returned toward the transmitter 
along a path similar to that followed on its way 
to the sunlit zone. This mechanism can be used 
to explain echo paths of long delay, such as some 
which have been ascribed to “G’’ layers. 

In Fig. 14 we have shown the average F layer 
curvature in the sunrise zone deduced from Fig. 
12, with the addition of a sporadic E layer. 
Between these layers we have indicated a ray 
which enters the region from the left and which 
finally meets a reflecting layer perpendicularly at 
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Fic. 14. A cross section of the ionosphere in the sunrise zone, showing the path 
of a ray which may be returned to the neighborhood of the transmitter. 


n « 
le 
1e 
tO 
d 
th 
to 
Ss. 
a- 
he 
th 
is 
is 
ch 
ns 
lly 
as 
of 
he 
on 
of 
of 
to 
rs. 
‘lic 
es, 
ice 
ese 
ers 
ure 
ind 
on, 
or 
the — 
in 
= 


P so that it will return along its outgoing path. 
This construction is intended to illustrate the 
mechanism of reflection of one echo of long delay 
which was observed, using a sweep-frequency 
technique, on November 7, 1938. This echo was 
delayed by almost exactly 0.1 second, corre- 
sponding to a total equivalent path of 30,000 km. 
The delay was not, within limits, a function of 
frequency, indicating that an explanation of its 
existence must be based upon the geometry of the 
reflecting surface rather than upon variations in 
ionic density. At the time of this record, which is 
not clear enough for reproduction here, the great 
circle distance from Cambridge to the point 
where a ray would meet the sunrise line perpen- 
dicularly was somewhat less than 12,000 kilo- 
meters. Thus, allowing for the vertical com- 
ponents of the ray path shown in Fig. 14, the 
reasonable assumption that the average velocity 
of propagation was slightly less than the velocity 
of light would enable us to explain this echo 
completely. The £ layer ionization was assumed 
because the record showed a tendency towards a 
family of reflections differing in their ‘apparent 
height”” by about fifty kilometers, which is 
approximately half the spacing between the two 
layers at the point where reflection occurred 
perpendicularly. A second ray, after following a 
path almost identical with that shown in Fig. 14, 
could be turned back (by making two successive 
reflections at the same angle of incidence) with- 
out having encountered either layer normally but 


0000 0100 0200 EST 


Fic. 15. A record showing, at the right, a group of G 
layer reflections. Note that the base line is below the 
bottom edge of the photograph, the lowest continuous 
trace being the first F layer reflection. Part of the pattern is 
repeated and normally the height scale would be repeated 


also. 
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having traversed the region between the layers 
once more (or less) than the ray shown in the 
figure. Thus a family of reflections having very 
long delays which differ by multiples of the time 
corresponding to the distance between the E and 
F layers could be returned to the transmitter. 
The (sporadic) E region ionization may have 
existed only in the sunrise zone, or may have 
extended a considerable distance along the path 
towards the transmitter. The latter assumption 
would help to explain the fact that the reflections 
were strong enough to be recorded, because 
propagation between the E and F regions should 
be subject to much less attenuation than propa- 
gation between the F region and the earth, where 
the ray would pass repeatedly through low level 
absorbing strata. 

Figure 15 is a reproduction of a record ex- 
hibiting an unusually uniform structure of “‘G 
layer” reflections which appeared at about 0200 
E.S.T. at the end of April, 1934. Here the 
apparent “‘height”’ of the group is about 1500 km. 
Since the use of 60-cycle pulses causes the whole 
reflection pattern to be repeated with a period 
corresponding to 2500 km of effective height, the 
equivalent distance is 1500 plus a multiple of 
2500 km. At this time the edge of the sunlit zone 
was about 3000 km north by east of Cambridge. 
Comparison of this pattern with Fig. 14 will 
indicate that the effective distance of the re- 
flections was about 4000 km, the additional 1000 
km being accounted for chiefly by the vertical 
components of the ray direction. The family of 
reflections referred to above is seen distinctly in 
this record, and their spacing is that to be 
expected if the height of reflection, at the 
point where the rays turned back, was about 
220 km. 

In conclusion, we should refer to a paper! which 
came to our attention after our oral presentation 
of this report. The authors of the paper con- 
cluded an attempt to explain certain com- 
plex echo patterns with the paragraph: 

“However, we have to consider the other 
possibility put forward by one of us, that at 
times the ionized regions may be an undulatory 


1G. R. Toshniwal, B. D. Pant, R. R. Bajpai and B. K. 
Verma, Proc. Nat. Acad. Sci. (India) 6, 161-174 (1936). 
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surface with various ledges here and there. The 
presence of such a structure along with the 
possibility of arrival of echoes from angles other 
than the vertical can easily account for the com- 
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plex echoes. If this view of the production of the 
complex echoes is accepted, we have to search 
for the agency which can cause such a structure 
of the ionosphere.” 
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Nuclear Isomerism in Selenium and Krypton 


ALEXANDER LANGsDORF, JR. AND EmiLio Secret 
Radiation Laboratory, Department of Physics University of California, Berkeley, California 
(Received November 4, 1939) 


A study has been made of the radioactive chain 


The last member is an isomer of the stable krypton of mass 83. An isomeric pair Se” ° * has 
been found and the two genetically related nuclear isomers have been separated chemically. 
Kr*® and Br*® are produced by the fission of uranium and thorium as well as by the bombardment 


of selenium with deuterons and slow neutrons. 


 - the course of an investigation of radioactive 
products of thorium and uranium fission, we 
observed a radioactive krypton growing out of a 
radioactive bromine with a period of 2.4 hours. 
This bromine period coincides with that assigned 
by Snell' to Br*, which he found to be the 
product of the decay of Se* produced by deu- 
teron bombardment of selenium. .He did not, 
however, observe radioactive krypton growing 
out of the bromine. We therefore repeated this 
deuteron bombardment of selenium and were 
able to extract from the products a radioactive 
isotope of krypton having the same period as 
the fission product, and emitting the same type 
of radiation. Since the radioactivities of selenium 
and bromine appeared rather complicated, we 
thought it desirable to go over the different 
points in detail. 

The main results of our investigation are the 
identification of an excited state of the stable 
Kr®, the identification and chemical separation 
of an isomeric pair in Se7® °F *!, and confirmation 
of the identification of Se*. Although discovered 
by Snell, there had been some uncertainty as to 
the true half-life of Se**. 


1A. H, Snell, Phys. Rev. 52, 1007 (1937). 


THe CHAIN Se®—Br®—Kr*8— Kr 


We have produced Se*® by means of slow 
neutrons with the reaction Se® (n,y) Se*® and by 
deuterons with the reaction Se*® (d,p) Se*. The 
latter reaction of course gives much more con- 
centrated samples and has been used for all the 
chemical separations. In order to determine the 
period of Se*,. repeated extractions of Br® are 
necessary to identify its parent in the complex of 
selenium activities. The Br® can be recognized 
easily from the fact that it is the only bromine 
isotope giving rise to an active krypton. 

In order to make repeated extractions of 
bromine, the bombarded selenium was converted 
to selenic acid free of the primary 34-hour 
bromine activity produced by the reaction 
Se® (d,2n) Preliminary experiments with 
the more rapidly prepared selenious acid had 
shown serious difficulties as a result of the 
apparent impossibility of achieving a complete 
extraction of bromine from this material, so that 
the longer chemical preparation of selenic acid 
was found to be necessary. To the selenic acid 
solution we added periodically five milligrams 
each of potassium bromide and iodide, and 
precipitated them with silver sulphate solution. 
Decay curves of the bromine activities thus 
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extracted are shown in Fig. 1, curves I-VII. 
Curves I-—V all indicate a half-life of 140 minutes 
which is assigned to Br*. In each of these curves 
the point marked with a cross represents the 
bromine activity extrapolated back to the time 
of extraction. Curve VIII, drawn through these 
crossed points, indicates that the period of the 
parent Se* is about 30 minutes. . 

The absorption curve of the Br® beta-rays 
gives an end-point of 1.05 Mev by Feather’s 


method.? 


— + Se® decay curve 
© Br® decay curve 


Activity (log scale) 


0.1 


Fic. 1. Decay curves of successive bromine extracts from 
deuteron-bombarded selenium. The crosses represent the 
bromine activity extrapolated back to the time of extrac- 
tion. The curve VIII a the crosses represents the 
decay curve of the parent selenium 83. 


Decay curves of the last two extracts of 
bromine from the active selenic acid are shown 
in curves VI and VII of Fig. 1. They are dis- 
tinctly anomalous in intensity and slope. Al- 
though we have not investigated this anomaly, 
it might be due to a short-lived bromine ac- 
tivity which grows out of a selenium activity 
with a period longer than that of Se*. Since the 
curves of Fig. 1 are compounded of the decay 
curves of Br® and another unknown activity, 
there is some uncertainty in the assignments of 


2 N. Feather, Proc. Camb. Phil. Soc. 34, 599 (1938). 
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Fic. 2. Activity of an emanating bromine sample inside 
an ionization chamber. Triangles represent activity due to 
bromine; crosses represent total activity. The growth of 
activity is all due to accumulation of pton. The dots 
give the activity due to krypton alone after 3.5 hours of 
— The line through the dots gives the period of the 

romine parent of the krypton activity. 


the periods of Se*® and Br® from the data of this 
figure. We have, therefore, made use of the 
simple activity of Kr®* to confirm the above 
assignments of periods to the Br® and Se®. 

Radioactive krypton can be extracted from a 
solution containing Br®* simply by boiling the 
solution. However, rather prolonged boiling com- 
bined with bubbling of a gas stream through the 
solution is required for a quantitative extraction. 
Although we have used this technique, we found 
it more convenient to use emanating solid 
samples when we wanted to make repeated 
extractions or when krypton was wanted with 
but little air admixture. We will describe later 
how to prepare emanating samples. 

The decay curve of krypton is a perfect expo- 
nential; we followed it over a factor of ten 
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thousand, and the half-value period is 113 


minutes. 
The relation of Kr®* to Br® has been shown by 
the following experiment. An emanating sample 


of radioactive bromine was suspended inside an 


ionization chamber which measured the sum of 
the bromine and krypton activities. The air was 
pumped out and replaced by fresh air at equal 
time intervals, so that each exhaustion was an 
extraction of krypton. The decrease of activity 
due to this extraction is, therefore, a direct 
measure of the krypton activity, and the residual 
activity immediately after letting in fresh air is 
bromine activity. In Fig. 2 the crosses represent 
the total activity, the triangles show the bromine 
activity, and the circles represent the krypton 
activity. Since the times of extraction of krypton 
were uniformly spaced, the curve of krypton 
activity is simply the decay curve of its parent 
substance. It will be seen that the period of the 
bromine so obtained agrees with that of the Br® 
found directly. 

As indicated before, we found it desirable to 
check the period of the Se* by the Kr® activity 
because we were not sure that the bromine decay 
curve was simple. This was done by extracting 
bromine periodically from active selenic acid 
with carbon tetrachloride. Half of the extracted 
bromine was converted into an emanating sample 
and half was measured directly. All the emanat- 
ing samples were allowed to accumulate krypton 
for 3.5 hours so that the activity due to krypton 
after this interval of time was proportional to 
the amount of parent bromine initially present. 
Fig. 3 shows the results of measurements of the 
krypton and bromine activities. The initial 
bromine and krypton activities for different 
extractions are in a constant ratio. This shows 
that practically all of the bromine activity in 
our extracts is due to Br® and affords also an 
independent check on the half-life of Se®. 

Krypton 83 emits a particularly soft radiation. 
We have been able to detect it only by intro- 
ducing the gas into an ionization chamber. In 
order to have an estimate of the energy of the 
electrons emitted, the active krypton was ad- 
mitted into a spherical ionization chamber of 
12.6 cm diameter, which had previously been 
evacuated. Air was then admitted and Fig. 4 
shows the ionization as a function of pressure. 


Similar curves have been obtained in x-ray work, 
and a discussion of their shape as a function of 
the absorption coefficient of the electrons is to 
be found in Bothe’s article.* His discussion refers 
to a plane condenser, whereas our experimental 
arrangement corresponded to a spherical con- 
denser, which introduces some uncertainty in the 
determination of the dimension of the vessel to 
be considered as corresponding to the plate 


ACTIVITY LOG SCALE 


Fic. 3. Successive extractions of bromine from selenium. 
ene and crosses same as in Fig. 1. Circles represent 
re 4 activity which grows out of bromine extracts in 


separation of the plane condenser. Nevertheless, 
the absorption coefficient varies so rapidly in 
the region concerned that a considerable error in 
the estimation of the absorption coefficient makes 
only a small error in the estimation of the 
energy. From Fig. 4 the absorption coefficient 


*W. Bothe, Geiger, Scheel, Handbuch der, Physik 
(Berlin, 1933), 23/2, 26. 
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Fic. 4. Ionization due to krypton electrons as a function 
of the air pressure in the spherical ionization chamber of 


diameter 12.6 cm. 


is about 1630 cm?/g, and the practical range is 
about 4.2 mg/cm?. Both these results give an 
energy of about 35 kev for the electrons. 

The natural interpretation of the observations 
is that the electrons are conversion electrons of a 
49-kev gamma-ray. This energy is obtained by 
adding to the observed electron energy the 
binding energy of the K level in krypton. This 
interpretation is further confirmed by the detec- 
tion of the characteristic x-radiation of krypton 
by absorption measurements on a gas sample. 
X-rays were detected with a ‘freon’? chamber. 
We found that arsenic absorbs the x-rays more 
strongly than selenium. Since the energy of the 
characteristic K radiation of krypton lies be- 
tween the absorption edges of arsenic and 
selenium, this experiment shows conclusively 
that the x-rays are characteristic of krypton. 

Since there is no doubt that the parent of this 
radioactive krypton is the bromine isotope of 
atomic weight 83 (i.e., Br®*), we must conclude 
that this krypton is a metastable excited isomer 
of stable Kr*®. 

In this connection, we want to point out that 
Kr® is known to have a nuclear spin of 9/2.‘ 
According to the computations of H. Hebb and 
G. E. Uhlenbeck,* the observed energy and half- 


4H. ergy Zeits. f. Physik 109, 349 (1938). 
5H. Hebb, G. E. Uhlenbeck, Physica 5, 605 (1938). 
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life indicates a change of spin of four, so that the 
excited state probably has a spin of }. 

We think it is within the possibilities of the 
cyclotron to prepare enough excited Kr® to 
observe the hyperfine structure in its spectrum 
and its change when Kr* reverts to the ground 
state. This would permit a direct test of the 
Weizsacker theory of nuclear isomerism. 


SELENIUM 79 oR 81 


The selenium activity is so complex that we 
thought it worth while to look for isomers, as a 
means of furthering the analysis, especially since 
we know that the chemical method of separation 
of nuclear isomers® is applicable in the case of 
tellurium’ which is chemically similar to selenium. 

Radioactive selenic acid and hydrochloric acid 
were mixed so that the final solution was four 


TIME HOURS 


Fic. 5. Successive extractions of the lower isomeric 
state of Se’*. * from the total selenium activity. Dots for 
separated isomer activities; crosses for the initial activity 

ive the period of the parent isomeric state approximately. 
otal fatenesbaonesnes activity (triangles) is on a difierent scale. 


nr Segré, R. S. Halford and G. T. Seaborg, Phys. Rev. 
Livingood and J. W. Kennedy, 


55, 321 am 
7G. T J. 
Phys. Rev 35 55, 794 (1939). 
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normal with respect to hydrochloric acid. The 
solution was then saturated with sulphur dioxide. 
Under such conditions selenic acid is reduced 
very slowly whereas selenious acid is reduced 
rapidly to selenium which can be filtered out of 
the solution. To this selenic acid solution we 
added at intervals of time selenious acid and 
separated the elementary selenium. Fig. 5 shows 
a decay curve of the parent selenium and of the 
extracts. The extracts show a half-life of 19 
minutes, whereas the parent selenium has a 
half-life very close to one hour. 

The assignment of this isomeric pair to Se7® °* *! 
is based upon the observation that it can be 
obtained also by neutron bombardment of 
bromine, which has only two stable isotopes, 
Br7? and Br* in equal amounts. The activity 
of selenium obtained from bromine by the (n,p) 
reaction is much less complex than the activity 
which one obtains by deuteron bombardment of 
selenium. Using this latter method of preparation 
of we could directly measure the half- 
life of the upper isomeric state with good 
accuracy. It is very close to one hour, in agree- 
ment with an observation of Snell.! 

We also checked that this isomeric pair of 
selenium activities does not decay into any 
active product, a result that is to be expected. 

These observations explain the discrepancy 
between the results of Snell,! and Bothe and 
Gentner.* The latter probably observed the 
lower isomeric state of Se*! °* 7* which they ob- 
tained by gamma-ray bombardment; Snell, on 
the other hand, had not distinguished Se* from 
the lower isomer of Se7**". 

By the chemical method we extracted the 20- 
minute lower isomeric state of Se7*: *' produced 
by bombardment of bromine and found thereby 
that we had extracted at least one-half of the 
beta-ray activity of the parent selenium. It 
must, however, be borne in mind that the yield 
of our separation, due to its chemical features, 
may be rather low. It is quite possible, therefore, 
that all of the beta-ray activity of Se7® * is 
due to the lower isomeric state. The absorption 
curve in aluminum of the beta-rays of both 
activities is the same and gives by Feather’s 
method? an upper limit of the energy of about 


*W. Bothe and W. Gentner, Zeits. f. Physik 112, 45 
(1939). 


Se 


Fic. 6. Beta-spectrogram of Se bombarded with deu- 
terons showing K- and L-conversion electrons ejected in 
the isomeric transitions of and 


1.5 Mev. These results indicate that the upper 
state decays to the lower one, which in turn dis- 
integrates to bromine. The conversion electrons 
which correspond to the transition from the 
upper to the lower state have been observed 
as a very soft component in the absorption 
curve of the electrons emitted by the one-hour 
activity and have been detected in a magnetic 
spectrogram. Both K- and L-conversion lines, 
are illustrated in Fig. 6. The energy of the 
K-conversion electrons is about 85 kev, corre- 
sponding to an energy difference between the 
isomeric states of 98 kev. The same value has 
been found by absorption measurements on the 
electrons.® 


CHEMICAL PREPARATIONS 
Selenic acid 


The bombarded selenium was dissolved in 
nitric acid containing some bromide ion, then 
evaporated to dryness to eliminate most of the 
strong primary bromine activities, redissolved 
and boiled with sodium hydroxide, then filtered 
to remove metallic contaminations such as 
copper from the target plate. The solution was 
acidified with nitric acid and treated with an 
excess of normal silver nitrate to precipitate 
silver selenite. This compound was washed, 
suspended in water, and the selenite was oxidized 
to selenate by treatment with bromine, of which 
the excess was removed by boiling. The filtered 
solution contained selenic acid free of primary 
radioactive bromine. 


Separation of bromine 


As an alternative method to the precipitation 
of silver bromide out of selenic acid solution as 
already described, we have added equivalent 


* We thank Mr. A. C. Helmholz for photographing the 
line spectrum of electrons. 
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quantities of sodium bromide and sodium bro- 
mate in order to liberate bromine, which was 
then extracted with carbon tetrachloride. This 
method, while not so efficient as the former 
precipitation technique, has the advantage of 
giving the bromine in a form well suited to the 
preparation of emanating samples. 


Emanating samples 

We prepared silica gel impregnated with silver 
nitrate by immersing the granular silica gel 
(8 to 16 mesh) in 0.1 WN silver nitrate solution 
for about 15 minutes. The gel was then drained 
and dried for several hours on a hot-plate at 
about 70°C. The radioactive bromine was pre- 
pared in the form of free bromine with several 
milligrams of bromine carrier, dissolved in about 
5 cc of carbon tetrachloride (free of any visible 
water layer but not specially dried). Several 
cubic centimeters of the silver nitrate-impreg- 
nated silica gel were added to this solution and 
shaken with it until the bromine color had com- 
pletely disappeared. This usually took about 15 
minutes. The silica gel was drained and dried on 
the hot-plate for about one-half hour. It was 
then ready to be placed in a vessel for collection 
of the emanating krypton. This method of 
preparing an emanating solid has been found to 
be equally applicable in the case of radioactive 
xenons which are produced by decay of iodine 
extracted from thorium and uranium fission 
products. 


Other observations 


The preceding work is not a complete investi- 
gation of the radioactivities produced by bom- 
bardments of selenium. Our work indicates there 
are several points requiring further study. Per- 


haps there is a radioactive bromine with a half- 
life less than one-half hour which grows from a 
radioactive selenium with a half-life considerably 
longer than that of Se*. This unidentified 
activity has already been discussed in connec- 
tion with Fig. 1, curves VI and VII. 

We have also observed an activity in selenium 
bombarded with deuterons with a half-life of 
more than a day. Since this longer-lived activity 
tends to be selectively extracted from the parent 
selenate solution when making the first isomeric 
separation, for Se7®°®, it may be that this 
longer-lived activity is itself the lower isomeric 


state of another as yet unknown selenium. Fig. 5 ° 


shows evidence on which these conclusions are 
based. Curve I, Fig. 5, shows a strong long-lived 
tail in the first isomeric extract. This tail is much 
less pronounced in later extracts. It cannot be 
due to an imperfect isomeric separation, because 
the tail actually shows a longer period than is 
shown by the fotal selenium activity in curve V 
of Fig. 5. However, the possibilities of the pres- 
ence of impurities in the selenium and of the 
formation of arsenic during the bombardment 
have to be borne in mind in appraising the results 
reported in this paragraph. 

Since Br® and Kr*® are observed in the fission 
of thorium and uranium, it is probable that Se* 
is also a product of the fission. It has not, how- 
ever, been observed directly. An extract of 
selenium from neutron irradiated uranium showed 
the presence of several different activities, which 
prevented a simple confirmation of the presence 
of the Se® activity. 

Our thanks are due to Professor E.O. Lawrence 
for his interest in this work. The financial support 
of the Research Corporation and of the Rocke- 
feller Foundation is gratefully acknowledged. 
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A Quantitative Determination of the Neutron Moment in Absolute 
Nuclear Magnetons 


Luis W. ALVAREZ, Radiation Laboratory, Department of Physics, University of California, Berkeley, California 
AND 


F. BLocu, Department of Physics, Stanford University, Palo Alto, California 
(Received October 30, 1939) 


The magnetic resonance method of determinng nuclear 
magnetic moments in molecular beams, recently described 


by Rabi and his collaborators, has been extended to allow © 


the determination of the neutron moment. In place of 
deflection by inhomogeneous magnetic fields, magnetic 
scattering is used to produce and analyze the polarized 
beam of neutrons. Partial depolarization of the neutron 
beam is observed when the Larmor precessional frequency 
of the neutrons in a strong field is in resonance with a weak 
oscillating magnetic field normal to the strong field. A 
knowledge of the frequency and field when the resonance is 
observed, plus the assumption that the neutron spin is }, 
yields the moment directly. The theory of the experiment 


is developed in some detail, and a description of the - 


apparatus is given. A new method of evaluating magnetic 


moments in all experiments using the resonance method is 
described. It is shown that the magnetic moment of any 
nucleus may be determined directly in absolute nuclear 
magnetons merely by a measurement of the ratio of two 
magnetic fields. These two fields are (a) that at which 
resonance occurs in a Rabi type experiment for a certain 
frequency, and (6) that at which protons are accelerated 
in a cyclotron operated on the mth harmonic of that 
frequency. The magnetic moment is then (for J= 4), 
u=H,/nH,. n is an integer and H,/H, may be determined 
by null methods with arbitrary precision. The final result 
of a long series of experiments during which 200 million 
neutrons were counted is that the magnetic moment of the 
neutron, #,=1.93;+0.02 absolute nuclear magnetons. A 
brief discussion of the significance of this result is presented. 


INTRODUCTION 


HE study of hyperfine structure in atomic 

spectra has shown that a large number of 
atomic nuclei possess an angular momentum and 
a magnetic moment. Since, according to the 
theories of Heisenberg and Majorana, protons 
and neutrons are recognized as the elementary 
constituents of nuclear matter, their intrinsic 
properties and particularly their magnetic mo- 
ments have become of considerable interest. The 
fundamental experiments of Stern and _ his 
collaborators! in which they determined the 
magnetic moments of the proton and the 
deuteron by deflections of molecular beams in 
inhomogeneous fields gave the first quantitative 
data of this sort. The approximate values which 
they gave for the two moments,’ 


bp=2.5, (1) 


1R. Frisch and O. Stern, Zeits. f. Physik 85, 4 (1933); 
I. Estermann and O. Stern, Zeits. f. Physik 85, 17 (1933); 
O. Stern, Zeits. f. Physik 89, 665 (1934); I. Estermann and 
O. Stern, Phys. Rev. 45, 761 (1934). , 

? We shall throughout this paper give magnetic moments 
in units of the nuclear magneton eh/4xMc (e and M 
=charge and mass of the proton, c= velocity of light), and 
angular momenta in units h=h/2x. The angular mo- 
mentum of the deuteron is taken to be 1, that of the 

ton and the neutron to be 34; the last is an assumption 
or which a direct experimental proof is still lacking, but 
which is generally accepted. 


0.8, (2) 


suggested that in all likelihood, one would have 
to ascribe to the neutron a magnetic moment of 
the approximate value 


(3) 


The negative sign in formula (3) indicates that 
the relative orientation of their magnetic mo- 
ments with respect to their angular momenta is 
opposite in the case of the neutron to that of the 
proton and the deuteron. 

The technique of molecular beams has been 
greatly developed during the last few years by 
Rabi and his collaborators;* their ingenious 
methods have allowed them to determine the 
magnetic moments of many light nuclei with high 
precision, and to establish the existence of an 
electric quadrupole moment of the deuteron. 
Their values for the magnetic moments of the 
proton and deuteron are 


Mp = 2.785+0.02, (4) 
wa=0.855+0.006. (5) 
They have also demonstrated that both moments 


*I. I, Rabi, Phys. Rev. 51, 652 (1937); I. I. Rabi, J. R. 
Zacharias, S. Millman and P. Kusch, Phys. Rev. 55, 
526 (1939). 
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are positive with respect to the direction of the 
angular momentum. 
An experimental proof that a free neutron 


_ possesses a magnetic moment, and a measure of 
‘its strength, could also be achieved in principle 


by deflection of neutron beams in an inhomo- 
geneous magnetic field. But while the great 
collimation required for this type of experiment 
may easily be obtained with molecular beams, it 


would be almost impossible with the neutron . 


sources available at present. Better suited for the 
purpose is the method of magnetic scattering, 
which was suggested a few years ago by one of 
us.‘ It is based upon the principle that a notice- 
able part of the scattering of slow neutrons can be 
due to the interaction of their magnetic moments 
with that of the extranuclear electrons of the 
scattering atom. In the case of a magnetized 
scatterer this will cause a difference in the 
scattering cross section, dependent upon the 
orientation of the neutron moment with respect 
to the magnetization, and particularly in the case 
of ferrornagnetics, it will cause a partial polariza- 
tion of the transmitted neutron beam. The mag- 
netic scattering of neutrons, and thereby the 
existence of the neutron moment, has been 
proved experimentally by several investigators, 
particularly by Dunning and his collaborators.® 
The magnetic scattering, however, can yield only 
a qualitative determination of the neutron mo- 
ment since the interpretation of the effect is 
largely obscured by features involving the nature 
of the scattering substance. 

Frisch, v. Halban and Koch® were the first to 
attempt to use the polarization of neutrons 
merely as a tool, and to determine the neutron 
moment by a change of the polarization, pro- 
duced by a magnetic field between the polarizer 
and the analyzer. Such a change should indeed 
occur, because of the fact that the moment will 
precess in a magnetic field; by varying the field 
strength, one can reach a point where the time 
spent by the neutrons in the field is comparable 


‘F. Bloch, Phys. Rev. 50, 259 (1936); 51, 994 (1937); 
Ann. Inst. H. Poincaré 8, 63 (1938); J. Schwinger, oa 
Rev. 51, 544 (1937). 

N. Powers, H. Carroll and J. R. 
Rev. 51, 1112 (1937); N. Carroll Bever 
and J. R. Dunning, Rev. 52, 38 (1937) ; "Pp, 
Powers, Phys. Rev. T1988). 

Frisch, H. v. Halban and J. Koch, Phys. Rev. 53, 


*R. 
719 (1938). 
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to the Larmor period. In this way, one could 
obtain at least the order of magnitude of the 
moment. Although these investigators have re- 
ported an effect of the expected type, yielding the 
order of magnitude 2 for the neutron moment, we 
have serious doubts that their results are signifi- 
cant. Their polarizer and analyzer consisted of 
rings of Swedish iron, carrying only their 
remanent magnetism (B= 10,000 gauss), while in 
agreement with Powers’ results, we were never 
able to detect any noticeable polarization effects, 
independent of the kind of iron used, until it was 
magnetized between the poles of a strong electro- 
magnet with an induction well above 20,000 
gauss. Although we cannot deny the possibility 
that, due to unknown reasons, their iron was far 
more effective for polarization at low values of 
the induction than that used by other investi- 
gators, we think it more likely that in view of 
their rather large statistical errors the apparent 
effect was merely the result of fluctuations. 

Although most valuable as a new method of 
approach, the experiment of Frisch, v. Halban 
and Kochrcould in any event give only qualitative 
results. The slow neutrons which one is forced to 
use emerge from paraffin with a complicated and 
none too well-known velocity distribution. The 
time during which they precess in the magnetic 
field will therefore be different for different 
neutrons and vary over a rather large range. 
Since it is that time which together with the field 
of precession determines the value of the moment, 
the latter will be known only approximately. A 
quantitative determination of the neutron mo- 
ment therefore requires an arrangement which 
does not contain such features. 


METHOD 


Sometime ago, we conceived of an experimental 
method which could yield quantitative data of 
this sort. The method was independently pro- 
posed by Gorter and Rabi,’ and most successfully 
used by the latter in his precision determinations 


71. I. Rabi, J. R. Zacharias, S. Millman and P. Kusch, 
Phys. Rev. 53, 318 (1938). During the summer of 1938 a 
first attempt to measure the neutron moment with such 
an arrangement was made at Stanford by Bloch, Bradbury 
and Tatel, using a 200- Ra-Be source of neutrons. 
However, it was. unsuccessful due to instrumental effects 
—— we hard to discover with the small intensity 
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of nuclear moments. Its principle consists in the 
variation of a magnetic field Hy to the point where 
the Larmor precession of the neutrons is in 
resonance with the frequency of an oscillating 
magnetic field. The ratio of the resonance value 
of Ho to the known frequency of the oscillating 
field gives immediately the value of the magnetic 
moment. 

The observation of the resonance point is 
based upon the fact that in its neighborhood 
there will be a finite probability P for a change of 
the orientation of the neutron moment with 
respect to the direction of the field Ho. Let this 
field be oriented in the z direction and let there be 


perpendicular to it, say in the x direction, an - 


oscillating field with amplitude H, and circular 
frequency® w, so that the total field in which a 
neutron is forced to move, is given by its 
components 


H,=H, cos (wt+ 8); H,=0; H,= Hyp. (6) 


The solution of the Schroedinger equation for a 
neutron with angular momentum 3 and magnetic 
moment yu gives the probability that a neutron, 
which at time ‘=0 in such a field had a z com- 
ponent =}4 of its angular momentum, will be 
found at the time =T7 with a value m= —}, in 
the form 


P=sin? [(uHiT/2h) 
X (1+ {24H 


where 


(7) 


AH = H,— = Ho—hw/ 2p (8) 


is the difference between the constant field Hy 
and its value at resonance, 


H,* = hw/ 2p, (9) 


for which the Larmor frequency 2H ou/h is equal 
to the frequency w of the oscillating field. Since 
the time T which the neutrons spend in the 
oscillating field will, for different neutrons, vary 
over a wide range, it will be a good approximation 


* P. N. Powers, reference 5, has employed a modification 
of this method to determine the sign of the neutron 
moment. Instead of encountering an oscillating field 
between polarizer and analyzer, his neutrons were acted 
ee by a field which appeared to them, in consequence 

their velocity, to be rotating with approximately their 
Larmor precession frequency. The sign of the moment 
was found to be negative, in oy with theory. 

* We shall always use circular frequencies although we 


sometimes refer to them merely as ‘‘frequencies.” 
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to substitute for the sin? in the numerator of (7) 
its average value 4. This means that, at reso- 
nance, complete depolarization of an originally 
polarized neutron beam will occur, and leads to 
the simplified formula 


(7a) 


Formulae (7) and (7a) have been derived by 
neglecting in the differential equation under con- 
sideration certain rapidly oscillating terms with 
respect to slowly varying ones; this neglect is 
justified as long as the condition 


(10) 


holds. In our final experiments, Hi/H» was less 
than 2 percent; so we felt justified in applying 
formulae (7) and (7a).* 

The dependence of P upon AH, and particularly 
its maximum value for which the field Hy has its 
resonance value H,* given by Eq. (9), is detected 
by a polarizer-analyzer arrangement, similar to 
that used by Frisch, v. Halban and Koch. We 
shall show that a beam of neutrons passing 
through two successive plates of magnetized 
material will change its transmitted intensity by 
an amount proportional to P. Consider two 
plates F, and F, both magnetized in the +z 
direction and call f;* and f=- (i=1,2) the fraction 
of neutrons with m=+} and m=-—4, respec- 
tively, transmitted by F;. f;* will be different from 
fi since the magnetic scattering will be prefer- 
ential for one of the two values of m, thus giving 
rise to a partial polarization of the neutrons. If 
there does not occur any change of polarization. 
between F; and F», the fraction of the neutrons 
transmitted by both will be given by 


thtr. (11) 


If, on the other hand, there is a probability P of a 
change of polarization between F; and Fs, the 


transmitted fraction will be 


(12) 
One will thus observe a relative change AJ of the 
intensity J of the transmitted neutrons, given by 


AI/I=(f'—f)/f=—aP, (13) 


* The corrections, introduced by this neglect are by no 
means of the order of 2 percent, but very much smaller, 
since they are essentially given by (H:/Ho)*. The mathe- 


matical t of these corrections will soon be published. 
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Fic. 1. Schematic drawing of electrical circuits. 


with 
(14) 


where a is a constant, characteristic of the mag- 
netic and geometric conditions of F; and F; 
only. We see that this constant a can also be 
determined independently of any possible changes 
of polarization. For two identical plates F; and 
F, with common thickness a and with the values 
=e—¢*+42(Ag =0 for unmagnetized ma- 
terial) as they follow from the theory of magnetic 
scattering, one can see that a@ also has the 
significance of the relative increase of the total 
transmitted intensity which one observes by 
magnetizing F; and F, to the stage in which the 
fractions f; have the values assumed in formula 
(14). Thus under ideal conditions, and if P 
assumes its maximum value } at resonance, the 
relative change in intensity will just be } of that 
observed when demagnetizing F, and F; in an 
experiment similar to that of Dunning and his 
collaborators. 

Once the maximum value of (13) has been 
experimentally established, i.e., the value of Ho 


‘has been found for which AH=0, or where 


according to (8) and (9) 


Hy=Ho* =hw/2y, (15) 
one will know the magnetic moment u in the form 
(15a) 


provided that the values of Hy* and the fre- 
quency w of the oscillating field are known. 
One can of course measure the latter by means 
of an ordinary standardized wavemeter and the 
former by the usual method of flipping a coil in 
the magnetic field and comparing the ballistic 
throw of a galvanometer with that produced by 
breaking the current in the primary of a known 
mutual inductance. We have, however, thought 
of an alternative method, which not only can be 
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carried to almost unlimited accuracy, but also 
has the highly desirable feature of yielding the 
moment to be determined directly in units of the 
nuclear magneton. It consists in comparing the 
frequency w=w, and the field H)>*=H, with the 
values w, and H,, where w, is the circular 
frequency of an oscillator driving a cyclotron 


which accelerates protons in a magnetic field H,. 


We have them from (15) 
(2H (16) 


if instead of using absolute units as in (15), we 
express the neutron moment in units of the 
nuclear magneton, i.e., if we write 


(17) 


On the other hand, for the ideal resonance of 
protons in a cyclotron, we have the relation 


w,p=eH,/ Mc. (18) 
Dividing (16) by (18), we now obtain 
(wn/wp)(H,/Hn). (19) 


This means that in order to determine the 
moment in absolute nuclear magnetons one has 
merely to determine the ratio of two frequencies 
and the ratio of two fields instead of determining 
w, and H, absolutely. In addition to this great 
experimental advantage, we see that no knowl- 
edge of the universal constants of physics is 
required. The value of Planck’s constant does 
not enter since the moment is expressed in 
nuclear magnetons and the value of e/Mc is 
eliminated by use of the cyclotron. 

If the radiofrequency generator for the 
cyclotron were a power amplifier driven by a 
harmonic of the “neutron oscillator,”’ the ratio 
w/w, would be simply an integer. In this case 
one would merely have to measure the ratio 
H,/H,, in order to know u,. This measurement 
can be reduced to the determination of the ratio 
of the wound areas (area X turns) of two flip 
coils, which can in turn be reduced to that of two 
currents. If the two coils, located in the same 
position in a current system and both connected 
through the same galvanometer circuit give the 
same ballistic throw, then the ratio of the 
currents for which this condition is reached gives 
the ratio of their wound areas. Once that ratio is 
known, one also has that of any two fields which, 
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using the two coils, give the same ballistic throw 
in the galvanometer. Gehrcke and v. Wogau"® 
have developed this technique, by the use of null 
indicator methods, to a point where the ratio of 
two fields may be determined to five significant 
figures. Although we have determined the ratios 
of the frequencies and the fields by essentially 
these methods the errors, entering from elsewhere 
did not justify the application of all their 
refinements described above. 


DESCRIPTION OF APPARATUS 


Up to the present time, almost all of the 
precision measurements of neutron intensity have 
been made using Rn-Be, or similar sources not 
subject to the erratic fluctuations occurring with 
artificial sources. Although it is possible to keep 
the cyclotron ion beam constant to one or two 
percent for long periods of time, this is not suffi- 
cient for the present work, where an accuracy of 
the order of 0.1 percent is needed. Also, one can- 
not assume that the neutron intensity is propor- 
tional to the beam current, since the distribution 
of the latter over the target, and the proportion 
of stray ions striking the insides of the dees are 
both important in this respect. This rules out the 
possibility of using an integrating beam meter, or 
auxiliary BF; chamber, as a monitor. The 
counting-monitoring circuit which was finally 
adopted for this experiment is shown sche- 
matically in Fig. 1. 

The neutrons were detected in a BF; chamber 
operated at 5000 volts. Its diameter and length 
were 9 cm and 15 cm, respectively, and the gas 
pressure was 50 cm of Hg. The arrangement of 
electrodes is quite similar to that of a chamber 
recently described by Powers.? A conventional 
linear amplifier of 4 stages fed the pulses into the 
modulation-distortion circuit, where all pulses of 
magnitude greater than a certain value were 
amplified to uniform size, and all below this 
height were suppressed. This circuit has been 
described in detail before, in connection with the 
production of monochromatic neutron beams," 
but was not essential for most of this work, since 
modulated beams were not used in the final 
experiment. The uniform pulses from this stage 

1 E. Gehrcke and M. v. Wogau, Verh. d. D. Phys. Ges. 


11, 664 (1909). 
UL. W. Alvarez, Phys. Rev. 54, 609 (1938). 
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were fed to a scale-of-16 of the type described by 
Reich,” which uses hard tubes. After passing 
through the scale, the pulses entered a switching 
circuit which fed them to either of two scales of 4 
plus a Cenco counter. The switching was carried 
out automatically at regular intervals of from 
2-10 seconds, depending upon the adjustment. 
At the same time, the switching circuit turned a 
magnet on and off (from magnetic scattering 
experiments), or keyed the oscillator (for the 
resonance experiment). 

The first circuits designed to carry out the 
switching used relays; these were found to be 
unsatisfactory due to differences in spring ad- 
justments, so an all electronic switching circuit 
was devised. This is shown in Fig. 2. Negative 
pulses were fed to the grids of the two right-hand 
6K7’s, each of which had another 6K7 connected 
to it in the familiar Rossi coincidence circuit. If 
the grids of either of these left-hand 6K7’s were 
negative, the plate current of its partner would be 
cut off, delivering a positive pulse to the corre- 
sponding scale-of-4. But if the grid was positive, 
the voltage pulse at the common plate connection 
would be too small to affect the scale. The circuit 
at the left is merely a device for making one grid 
positive and the other negative, and reversing 
this condition instantaneously at will. The two 
885’s are arranged in a conventional inverter, or 
scale-of-2 circuit. When current flows through 
the upper Thyratron and is blocked in the lower, 
pulses will pass through the switch into the upper 
scale-of-4. If a positive pulse is now applied to 
the grids of the 885’s, the scale-of-2 will invert, 
and the lower scale-of-4 will receive the neutron 
pulses. The positive pulse which flips the scale- 
of-2 originates in a relaxation oscillator of the 
type used in cathode-ray oscillograph sweep 
circuits. It was found that the interval between 


we 


Fic. 2. Electronic switching circuit. 
12H. Reich, Rev. Sci. Inst. 9, 222 (1938). 
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pulses was not constant when R and C were 
large enough to make the time constant 5 
seconds; so the frequency was increased, and the 
pulses from the oscillator were fed through a 
scale of 32 into the switching circuit. Fluctuations 
in the time interval were reduced to a negligible 
amount by this procedure. The length of the 
interval could be varied by means of a variable 
resistance in the oscillator circuit. Small relays in 
the plate circuits of the Thyratrons made it 
possible to control the magnets or R.F. oscillator 
in phase with the counter switching. A lag in 
these operations merely reduced slightly any 
real effect, while a lag in the counter switching 
would introduce a spurious effect of the same 
magnitude. 

On first thought it would seem that the 
statistical errors would be considerably increased 
by the method of switching a scale-of-16 from 
one scale-of-4 to another. But calculation of the 
effect shows that it is of no importance, and our 
experimental fluctuations have always checked 
with those to be expected on the assumption that 
this arrangement is equivalent to two separate 
scales-of-64. We have also made runs to test the 
equality of the counting time under the two 
conditions, by feeding the output of a beat fre- 
quency oscillator into the distortion amplifier. It 
was always found that the two Cenco counters 
recorded equal numbers of counts per cycle of the 
switching. 

We spent several months developing our 
polarizer-analyzer equipment. Our first attempts 
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were along lines similar to those used by Frisch, 
v. Halban and Koch. Flat rings of iron with a 
diameter of 1 meter, a thickness of 1 cm, and a 
width of 10 cm were wound with copper wire in 
toroidal fashion to minimize the stray field. We 
were unable to observe any difference in the 
transmitted neutron intensity, which depended 
upon the magnetization of the iron. This is in 
agreement with the Columbia school which has 
shown that this iron (Armco) has a sharp mag- 
netic effect threshold very near the saturation 
value, which most likely was not reached with 
our low magnetizing currents. 

Since they also report that Swedish iron shows 
neutron magnetic effects at lower fields, we next 
wound bars of Swedish iron with flat copper strip 
through which 120 amperes could be passed. The 
bars were 1.5 meters long, 10 cm wide, and 1.5 cm 
thick. Although the H in the iron was 250 gauss 
(B=20,000), we were still unable to observe a 
difference in neutron transmission. We finally cut 
a piece from one of the Swedish iron bars, and 
magnetized it strongly between the pole pieces of 
a Weiss magnet and observed the magnetic 
effect. With an iron ‘piece 4 cm thick, we obtained 
differences of about 6 percent in transmission. 
With half of the iron in the Weiss magnet, and 
the other half in a similar strong magnet, we 
confirmed Powers’ results on double trans- 
mission, i.e., that the change in intensity was 
about the same for the pieces in one magnet, as 
when they were separated in two, at some dis- 
tance. Many experiments were made to determine 
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Plan of the apparatus. Section taken 85 cm above floor of cyclotron room. 
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the optimum thickness of iron, and the mag- 
netizing currents needed to obtain the maximum 
effects. 

It was then possible to consider the construction 
of the precessing field, and the oscillator system. 
A plan of the apparatus is given in Fig. 3. 
Neutrons from the Be+D reaction were slowed 
to thermal velocities in the simple howitzer, and 
diffused down a cadmium-lined tube through the 
water tank to the polarizer. After passing through 
the precessing and flipping field, and the analyzer 
magnet, they were recorded in the BF; chamber. 
The large blocks of paraffin which defined the 
beam made the fast neutron background quite 
small. The iron plates reduced the thermal 
neutron intensity to about 5 percent, but the fast 
neutron background, as measured by covering 
one of the plates with Cd, was still only 25 
percent of the total counting rate. This small 
background made it unnecessary to use the 
backgroundless, modulated beam, as was origi- 
nally planned. The somewhat increased magnetic 
effect at low neutron velocities was not enough 
to offset the decrease in intensity which attended 
the use of the modulation. In the earlier stages of 
the experiment, and actually in one of the 
resonance experiments, 120-cycle modulation was 
used successfully, but it was abandoned in favor 
of the greater intensity available without 
modulation. 

The oscillator was of the conventional Hartley 
type, using an Eimac 1507, self-excited. The 
frequency was found to remain constant to better 
than 0.05 percent for days. The solenoid which 
produced the oscillating magnetic field was wound 
with a 9-cm square cross section and 15 cm 
height. The conductor was flat copper strip 
0.500.010” wound so that adjacent turns 
almost touched. Neutrons passed through the 
strip on entering and leaving the region of the 
oscillating field H,. The solenoid was connected 
in parallel with a section of the tank coil chosen 
to give 11 amperes magnetizing current as meas- 
ured on an R.F. ammeter in the solenoid circuit 
(H,:~10 gauss). The oscillator was keyed by 
means of a relay in the primaries of the plate 
supply transformers. 

The steady precessing field Ho was at first 
produced by a pair of Helmholtz coils plus a pair 
of cylindrical bars of iron to concentrate the flux. 
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This field was quite inhomogeneous, and was used 
for the first two resonance dips only. For the 
later measurements, we used a pair of large coils 
which had energized the magnet of the original 
i-mv cyclotron of Lawrence and Livingston. 
These magnet coils were loaned to us by Pro- 
fessor R. B. Brode, who had been using them in 
connection with a cosmic-ray cloud chamber. An 
arrangement of iron was found which produced a 
field of about 600 gauss +1 gauss throughout the 
volume occupied by the R.F. solenoid. The fluc- 
tuations in the laboratory d.c. supply were too 
great to allow its use in the coil circuit, so we 
utilized the constant current supply of the 
cyclotron magnet. This automatic stabilizer keeps 
the. voltage fluctuations across the cyclotron 
magnet coils to less than 0.1 percent and gradu- 
ally increases the value of the stabilized voltage 
to compensate for the heating of the coils. The 
magnet coils were then connected in parallel with 
the cyclotron coils through a variable resistance. 
The current was watched constantly, and the 
resistance was varied slightly when necessary, to 
allow for the differential heating of the two sets of 
coils. 

A pair of large grounded copper plates were 
placed against the polepieces, so that the magnet 
coils were shielded from the R.F. fields of the 
solenoid. If this shield was not grounded, radio- 
frequency currents were induced in the d.c. lines 
and fed back into the stabilizer. A partial 
rectification of these currents would cause the 
stabilizer to change the magnet current, and 
therefore throw the cyclotron out of resonance. 
The resulting drop in the neutron intensity would 
be in phase with the switching cycle, and would 
therefore produce an apparent effect of the same 
sign and magnitude of that to be expected. 
Another cause of apparent effects was the inter- 
action of the R.F. oscillator with the relaxation 
oscillator in the timing circuit. When the scale- 
of-32 was not properly grounded, it would count 
an extra pulse whenever the oscillator was turned 
on. This gave a 3 percent apparent effect of the 
correct sign and magnitude, due to the change in 
switching period. Several other spurious effects 
were encountered at various times during the 
experiment; it is due to the large intensities of 
neutrons available that we were able to trace 
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down and eliminate them after comparatively 
short counting times. 

The test for reliability of the apparatus in 
which we placed most faith was the following: A 
series of observations was taken at 20ya of 
deuterons, and it was desired to know whether 
the effect, e.g., difference in neutron counts with 
and without the oscillator on, was real. The 
current was then increased to 80 amperes, and 
the analyzer was covered with a cadmium sheet. 
-The counting rate in all parts of the amplifier, 
from the BF; chamber to the final mechanical 
counters remained the same as before, and all 
relays, magnetic and electric fields were operating 
in the same manner. The only difference was that 
the counted neutrons did not for the most part 
pass through the polarizer and analyzer, and that 
part of high velocity which did could not exhibit 
any magnetic effects. Any difference between the 
relative counting rates under these two condi- 
tions was ascribed to magnetic effects of the 
thermal neutrons. It was also possible in this way 
to eliminate causes of systematic effects. The 
closing of relays would occasionally jar the BF; 
chamber, and the resulting noise would introduce 
extra counts in synchronism with the timing 
cycle. Several effects of this nature were dis- 
covered and eliminated after long blank runs. It 
was also important at these high counting rates 
to guard against saturation in the BF; chamber. 
The intensity was always adjusted so that the 
counting rate was still a linear function of the 
deuteron current. 

These sources of trouble are mentioned to 
emphasize the great desirability of powerful 
neutron sources in experiments of this sort. In the 
final set of experiments, in which all these 
spurious effects were eliminated and on which our 
value of the neutron moment is based, only about 
4 million neutrons were counted. This number of 
counts could be taken with Rn-Be sources of 
moderate strength in a reasonable length of time, 
but it would be almost impossible to track down 
obscure systematic effects in the manner outlined 
above. We counted more than 200 million 
neutrons in test runs before we felt sure enough of 
the apparatus to have confidence in the results 
outlined below. Perhaps the most convincing 
proof of the reality of the effect was that we 
varied Ho and w over a factor of almost 10, and 
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each time the frequency was increased, the magnetic 
field setting for the center of the new dip was found 
within the accuracy, to be expected at the value 
calculated from the previously determined value 
of o* / ®. 


MEASUREMENTS 


In the first experiments with the small 
Helmholtz coils, the magnetic field was measured 
with the aid of a fluxmeter-search coil combi- 
nation. Although somewhat crude, the measure- 
ments did show that Ho*/w was constant for two 
frequencies differing by about 20 percent, and 
that the effect could therefore be attributed to 
the magnetic moment of the neutron, as outlined 
above. The shallowness of the resonance dip 
observed in these first experiments (AJ/J=0.6 
percent) could be attributed to the inhomogeneity 
of the field Ho. 

When the large magnet coils were installed, it 
was clear that a more accurate method of 
measuring the field would be necessary. A brass 
flip coil was therefore constructed in such a way 
that it could be alternately lowered into a 
standard position within the R.F. solenoid, or 
withdrawn to allow the neutrons to pass through. 
Stops were placed so that the coil turned through 
exactly 180°, as determined by optical means. 
The coil was turned by spring tension, after 
being cocked by hand. The leads from the coil 
were connected in series with a ballistic galva- 
nometer of very long period, a Hibbert magnetic 
standard, and a search coil at the center of a 
standard solenoid. A knowledge of the constant 
of the standard solenoid plus the wound area of 
the flip coil and the search coil allowed the 
strength of the field to be ascertained. The 
current through the standard solenoid was meas- 
ured by means of a standard resistance, standard 
cell, and potentiometer. 

Before and after each counting period, the 
relation between the current through the coils, 
and the magnetic field was determined. This con- 
sisted in noting the ratio of the throws of the 
ballistic galvanometer d/s, where d was the 
deflection due to the turning of the flip coil, and s 
was the throw of the magnetic standard. Effects 
due to any nonlinearity of the galvanometer were 
minimized by adjusting this ratio close to one by 
rewinding the standard. d/s was determined as a 


| 

| 

} 

| ce 
tc 
it 
th 
| 
| m 
ce 


NEUTRON MOMENT 


14.60 


Fic. 4. Proton resonance curves. The potential drop in 
millivolts across the magnet current shunt is plotted as 
abscissa against the electromotive force in millivolts of the 
acne ad used in the probe on which the proton 
current fell. 


function of the magnetizing current, and the 
neutron intensity was also measured as a function 
of this current. The resonant field was then found 
from a comparison of these two curves. 

The frequency of the oscillator was measured 
on a General Radio precision wavemeter. A 
knowledge of the field and frequency yielded the 


magnetic moment of the neutron by substitution 
in Eq. (15a). 

The most trustworthy value of the moment is, 
however, given in absolute nuclear magnetons by 
the comparison method, outlined above. Al- 
though we ultimately trusted more this method it 
was gratifying to obtain well within the experi- 
mental error the same value for u, by the two 
entirely different methods of calibration. For this 
purpose it was necessary to construct a new flip 
coil which when turned through 180° in the 
proton resonance field of the 60-inch cyclotron 
would give a d/s almost equal to that obtained in 
the measurement of Ho*. From a knowledge of 
the approximate values of the two magnetic 
fields involved, and the wound area of the original 
coil, it was a simple matter to construct such a 
coil. This coil was wound on a hard rubber form 
to ascertain the absence of magnetic impurities in 
it or in the brass used for the first coil. 

It was then necessary to determine the ratio of 
the wound areas of the two flip coils. This 
measurement was carried out in the following 
manner: The two coils were mounted near the 
center of a pair of carefully constructed Helmholtz 
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coils, which were designed for accurate cloud- 
chamber beta-ray curvature measurements. The 
constant of the coils did not enter the calculation ; 
it was merely assumed that the field was directly 
proportional to the current through the coils. 
This current was measured in terms of the 
potential drop across a standard resistance in the 
circuit, on a Wolff potentiometer standardized 
with a Weston cell. The small coil was flipped in 
the strongest field available, and the current was 
then reduced until the large coil gave the same 
ballistic throw as the small one in the higher 
field. (The effect of the earth’s field was sub- 
tracted from all readings.) The ratio of the wound 
areas of the two coils was then given by the ratio 
of the currents as determined on the potentiome- 
ter. The two coils were connected in series with 
the galvanometer; so the total resistance of the 
circuit remained constant. The current ratio was 
found to be constant for many values of the 
fields. To test for magnetic impurities in the 
search coil, the coil was connected in series 
opposition, and the large coil was turned through 
such an angle that when a current was broken in 
the Helmholtz coils no kick was observed on the 
galvanometer. This condition prevailed for a 
wide range of magnetizing currents. A very 
sensitive galvanometer was used in these meas- 
urements; so this zero method gave us confidence 
that no appreciable ferromagnetic contamination 
existed in the flip coils, one of brass and the 
other of hard rubber. 

The new coil, ballistic galvanometer, and 
magnetic standard were then carried to the 
Crocker Radiation Laboratory, where the mag- 
netic calibration was made in terms of the 60-inch 
cyclotron field. The rod carrying the flip coil 
projected into the cyclotron chamber through a 
window so that the coil could rotate in the gap at 
a point about 25 cm from the center of the 
cyclotron. Tests with a large search coil con- 
nected to a fluxmeter showed the field to be 
constant throughout this central region to better 
than 0.1 percent. The magnetization curve of the 
cyclotron magnet was then measured in the 
region around the known proton field. The 
measurements consisted in a determination of d/s 
against the potential drop across the magnet 
current shunt, as measured on the Wolff 
potentiometer. 
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Fic. 5. Neutron resonance dip. The magnet current in 
arbitrary used is plotted against the fractional change 
(AI/I) of the intensity of the neutron beam. 


The cyclotron was then evacuated and hydro- 
gen let into the chamber. The heating effect of 
the proton current to a probe placed at the point 
occupied by the flip coil was then measured as a 
function of the potential drop across the magnet 
current shunt. A comparison of these two curves 
gave the d/s value for the acceleration of protons 
at the frequency in use at the time. The effects of 


hysteresis on the resonance curve were found to 


be negligible, by observing it when the magnet 
current had been increased from zero to the proton 
value, and also when it had been carried through 
the cycle—zero, deuteron-field, proton-field. 

The proton resonance curve, Fig. 4, was meas- 
sured at two values of the dee voltage, and the 
center of gravity of the curve was found to move 
toward higher fields as the voltage was reduced. 
This was shown not to be caused by a change in 
the cyclotron frequency, by measurements taken 
during each run. The increase in width of the 
curves with voltage is due to the fact that at 
higher voltage, the protons have to circulate 
fewer times. The shift of the peaks is due mainly 
to the focusing properties of the cyclotron. Since 
the magnetic field was uniform out to the probe, 
the focusing was entirely electrostatic, and this is 
known to be dependent upon the phase of the 
protons with respect to that of the oscillator." 


Phys. Rev. 53, 408 (1938). 
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Professor E. McMillan and Mr. R. Wilson have 
made a study of these effects in this laboratory, 
and they agree that the proper method of 
choosing the magnetic field, H,, would be to runa 
large series of resonance curves with decreasing 
dee voltage, and extrapolate the resonance field 
to zero dee voltage. Changes in the apparatus 
made this impossible, so they have suggested that 
we choose a point midway between the maximum 
of the curve taken at lower dee voltage and its 
high field intercept. Since we have seen both 
experimentally and theoretically that the ex- 
trapolated resonance curve will lie somewhere 
between these two points, separated in field by 1 
percent, the maximum error introduced in this 
way could be only 0.5 percent. 

The ratio of the cyclotron frequency to the 
neutron frequency was determined with the aid 
of a Super-Pro short wave radio receiver. The 
second harmonic of the cyclotron oscillator was 
tuned in on this receiver in the same frequency 
range as the sixth to tenth harmonics of the 
neutron oscillator. The agreement between the 
various values of the frequency of the neutron 
oscillator as determined from these measurements 
attests to the linearity of the frequency scale of 
the receiver. The receiver was thus used merely 
as an uncalibrated, linear wavemeter. 


RESULTS 


Since the width of the resonance dip in gauss is 
constant for all values of the precessing field 
(Eq. (7a)), it is obvious that it is desirable to 


TABLE I. Data for neutron resonance curve (Fig. 5). Each 
reading is based on 120,000 counts. Magnetic field i 
neutron resonance, 622 gauss; frequency of oscillator, 1543 
kilocycles; magnetic moment of the neutron, 1.94 nuclear 
magnetons; d/s for coil 1 in neutron field, 0.866; d/s for coil 
2 in proton field, 0.844. 


CURRENT PERCENTAGE EFFECT MEAN EFFECT 
1.67 —0.8, +1.0 +0.1 +0.3 
1.72 +0.6, —0.9, —0.8 —0.4 +0.23 
1.76 —0.5, —0.9, —2.1, —1.1, —1.1 +0.16 

—1.2, —0.8 
1.78 —2.0, —2.2, —1.6, —1.5, —1.4, —1.52+0.13 
—0.5, —1.7, —1.2, —1.6, —1.5 
1.79 —1.7, —1.7 —-1.7 +0.3 
1.80 —1.4, —0.5, —0.7, —1.5, —0.71+0.16 
+0.1, —0.3 
81 —1.5 —-1.5 +04 
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perform exploratory experiments at low pre- 
cessing fields, to find an approximate value of 
the magnetic moment. After the relatively wide 
dip has been found at low fields, it is possible to 
increase the precision of the determination by 
raising the field and the frequency in the same 
ratio. Our first measurements were made at 260 
and 315 kc, and merely served to show that the 
value of the moment was close to 2 nuclear 
magnetons. When the large homogeneous field 
was installed, we repeated the measurement at 
315 kc; after observing the dip again at 810 kc, 
we increased the frequency to 1843 kc, where the 
field calibrations were made. Since the resonance 
dips at all the lower frequencies were consistent 
with that at 1843 kc, we will not describe those 
measurements in detail. The curve at 1843 kc 
is shown in Fig. 5. 

The experimental data from which this curve 
is plotted are given in Tables I, II and III. It 
will be noticed that two of the points were not 
plotted in the curve. These are the values at 1.79 
_ and 1.81 amperes, which have very large probable 
errors, due to the small numbers of neutrons 
counted. These two points have no appreciable 
effect upon the position of the minimum, and 
although they lie far off the curve, their positions 
are quite consistent with their probable errors. 
The curve was obtained on two separate days a 
week apart, and since no systematic difference 
could be noted between the two sets of data, they 
have been plotted together. The value of the 
magnetic moment is then given from Eq. (19), 


n= (1.856/9.685) X (10.60 X0.844/0.886) 
= 1.935 nuclear magnetons. 


The probable errors entering from various sources 
are: position of dip, 0.5 percent; galvanometer 
readings, 0.2 percent ; coil calibration, 0.2 percent; 
frequency calibration, 0.2 percent ; and cyclotron 
calibration, 0.6 percent. The final result may 
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for equal galvanometer throws. 


RatTIo 


10.61 
10.58 
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TaBLe III. Frequency calibration: reading, R, on radio 
receiver for nth harmonic. 


NEUTRON OSCILLATOR CYCLOTRON OSCILLATOR 


n R 
2 19.37 


R/n 
9.685 


All values in megacycles. 


then be expressed as 
bn = 1.935+0.02 
absolute nuclear magnetons. 


DISCUSSION 
The now rather accurately known values 


Mp=2.785+0.02 ua=—1.935+0.02 
ba = 0.855 0.006 


of the magnetic moments of proton, neutron and 
deuteron are of considerable interest for nuclear 
theory. The fact alone that yz, differs from unity 
and xu, differs from zero indicates that, unlike the 
electron, these particles are not sufficiently 
described by the relativistic wave equation of 
Dirac and that other causes underly their mag- 
netic properties. 

Whatever these causes may turn out to be one 
has to notice that there holds to well within the 
experimental error the simple empirical relation 


(20) 


This relation is far from being obvious and it 
would in fact seem rather surprising if it were 
rigorously satisfied. To explain it in simple terms 
one would have to make both the following 
assumptions: 

(a) The fundamental state of the deuteron is a 
3§ state so that there are no contributions to yu, 
arising from orbital motion of the particles. 

(6) The moments and are “additive,” 
i.e., their intrinsic values are not changed by the 
interaction of the proton and the neutron, 
forming a deuteron. 

The first assumption has been disproved by the 
recent discovery" that the deuteron possesses a 
finite electric quadrupole moment which is 


“J. M. B. Kellogg, I. I. Rabi, N. F. Ramsey and J. R. 
Phys. Rev. 55, 318 (1939). 
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incompatible with the symmetry character of a 
pure *S state. The second cannot be discarded on 
an experimental basis but it ceases to be plausible 
if one admits the possibility, that ultimately the 
same causes may underly both the magnetic 
properties and the mutual binding forces of the 
proton and the neutron. 

It is conceivable that the departure from any 
one of the two assumptions (a) and (6) would 
separately cause a considerable deviation from 
(20) but that for unknown reasons both together 
cancel each other very closely. Until reliable 
estimates of these deviations can be obtained 
we consider it, however, more likely that neither 
of them amounts to more than a few percents. 
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The methods of contact potential measurement de- 
scribed in previous reports of this series have been extended 
to zinc, for which the published work function data are 
discordant. Measurement was by the retarding potential 
method with variable anode, in tubes sealed from the 
pumps and gettered with vaporized barium. Barium was 
the reference metal and the measured surfaces of both 
barium and zinc were prepared by fractional distillation 
followed by revaporization of the middle fractions. In 
observations on three tubes and fifteen pairs of Ba-Zn 
surfaces the observed contact P.D., mean value 1.76 v, was 
reproducible and constant to +0.02 v. Combination of this 
value with the work function of barium, now well estab- 
lished at 2.52 ev, assigns the value 4.28+0.02 ev to the 
work function of zinc at room temperature. The initial 


HIS paper reports further results in a 

program of study of the external work 
functions of the pure metals as determined by 
measurement of their contact differences of 
potential with respect to a suitable reference 
metal, barium. Previous work under this pro- 
gram! has shown, (1) that the work functions of 
barium surfaces are reproducible to at least 


1P. A. Anderson, Phys. Rev. 47, 958 (1935); 49, 320 
(1936) ; 54, 753 (1938). 


measurement on a surface was taken within 8 seconds of its 
deposition. As a consequence of this small time lapse it 
becomes possible to answer the fundamental but neglected 
question of whether a work function measured under 
optimum vacuum conditions is characteristic of the clean 
metal or of a gas-contaminated surface in equilibrium with 
the residual gas. The rate of deposition of oxygen on zinc is 
estimated from the behavior of a tungsten “‘test surface” 
in a comparable vacuum. Since the regimes of measurement 
and of equilibrium gas film formation overlap and no change 
of work function with aging is detectable, it is concluded 
that the observed work function is characteristic of clean 


' zinc. From the point of view of the adsorption processes 


involved, the measurements show that the mean time of 
sojourn of oxygen is much less on zinc than on tungsten. 


+0.01 ev when the surfaces are formed by 
fractional distillation in a sealed-off tube, and 
(2) that measurements of work functions by the 
contact potential method are now capable of the 
precision and accuracy characteristic of the best 
photoelectric work. Since the contact potential 
method, unlike the thermionic, is applicable to 
all the metals and, unlike the photoelectric, 
yields measurements of uniformly high sensi- 
tivity, it appears to offer the most promising 


| 


CONTACT DIFFERENCE OF POTENTIAL 


road to a systematic knowledge of the work- 
function characteristics of noncomposite metal 
surfaces. 

For the purpose of testing our method of 
contact potential measurement, previous studies 
have been made upon metals for which the 
photoelectric data are especially reliable. The 
method is now being extended to metals for 
which the photoelectric data are contradictory, 
or for which no work function studies have yet 
been made under adequate experimental con- 
ditions. Zinc is one of the numerous common 
metals for which these data are discordant. 
Published values? range from 3.08 to 4.24 ev and 
the agreement has scarcely been improved by 
the most recent work. DeVoe,’ using the photo- 
electric method, has obtained 4.24 ev for the 
Fowler threshold while Oatley,* employing a 
contact potential method, has reported a work 
function of 3.40 ev. Differences in the crystal 
structure of the zinc surfaces, determined by 
variations in the methods employed by different 
investigators for preparing their specimens, are 
no doubt partly responsible for the divergent 
values reported, but structural differences cannot 
account for work function deviations greater 
than a few decivolts. As has been pointed out in 
previous papers of this series, it is desirable that 
a standardized technique be selected for the 
preparation of each metal which is subjected to 
study, that the work function of this standard 
surface be determined under the most rigorous 
vacuum conditions, and that studies of the 
structural dependence of the work function be 
then referred to this standard surface. It has 
been suggested, also, that in work on most of 
the common metals, standard surfaces of repro- 
ducible structure can be best prepared by 
fractional distillation of the metal in.a gettered 
vacuum, with deposition on a given substrate at 
given substrate temperature. This choice of 
technique, borne out by several extended series 
of measurements in this laboratory, has received 
further support from photoelectric measurements 
made at Northwestern by Cashman and his 
associates, who have developed a similar tech- 


* For a bibliography of the earlier work see A. L. Hughes 
and L, A. DuBridge, Photoelectric Phenomena (McGraw- 
Hill, 1932), p. 76. 

C.F. DeVoe, Phys. Rev. 50, 481 (1936). 

*C. W. Oatley, Proc. Roy. Soc. A155, 218 (1936). 
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nique. Recent results on magnesium films pre- 
pared by this standardized method and measured 
in independent photoelectric and contact-poten- 
tial studies have given work function values 
which agree to 0.01 ev.® It was the purpose of 
the present work to extend this method to the 


Fic. 1. The tube. First and second vaporizers, A and B, 
are mounted on single 3-lead press, one of 60-mil leads 
common. A (rotated 90° in sketch) comprises 20 x 20 x 5- 
mm pan of 1-mil Mo foil, with flat helical heater of 20-mil 
W wire. B is single 25X30-mm sheet of foil; Mo foil for 
barium, W foil for Zn (see text). Target assembly TT is 
of glass; rotatable on 60-mil W lead by magnetic traction 
on glass-enclosed armature. Contact C is 10-mil Pt wire 
embedded in face of target. Vapor shield S, with porthole 
P, is of glass with glass-enclosed armature. Electron gun 
E of 4-mil Ta foil; emitter 1-mil W foil 1.5 mm wide. 
A second vaporizing bulb, connected at V, is identical in 
design to the bulb shown. 


5R,. J. Cashman, Phys. Rev. 54, 971 (1938); P. A. 
Anderson, Phys. Rev. 54, 753 (1938). 
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determination of the contact difference of poten- 
tial between barium and zinc. Both metals were 
fractionally and multiply distilled and condensed 
on glass at room temperature. 


EXPERIMENTAL PROCEDURE. TUBE DESIGN 


The general features of the method and the 
electrical circuits have been described in pre- 
vious reports. Details peculiar to the present 
measurements may be described in connection 
with the design of the tube, Fig. 1. The main 
improvement over our previous tubes lies in 
the construction of the first and second vapor- 
izers and in their thorough shielding to eliminate 
any possibility of contamination of either metal 
by the other during the vaporizations. The 
first vaporizer A, providing for rapid fusion and 
controlled vaporization of a large charge (5-10 g) 
of metal by radiation heating, was an adaptation 
of a well-known method. The second vaporizer B 
is of new design which provides (1) a condensing 
surface of large area with consequent high 
efficiency in terms of the amount of metal 
finally utilized to that available in the middle 
fractions, and (2) very exact control of the rate 
at which vapor is delivered by the second 
vaporizer to the final condensing surface. This 
vaporizer was made by converting a 25 X30-mm 
sheet of 1-mil tungsten or molybdenum foil into 
a grid by parallel cuts with a razor blade, as 
indicated at B in the plan sketch of Fig. 1. 
The grid was connected to its 60-mil tantalum 
current leads with tantalum rivets and held 
taut by a simple spider of 10-mil tungsten wires 
embedded in a Nonex bead and hooked through 
small holes in the grid,- as indicated in the 
elevation sketch of Fig. 1. Heating current for 
the grid was supplied by a step-down trans- 
former with continuously variable secondary 
voltage (Variac 80-A). The tubular shield S, 
mounted on a simple glass-to-glass thrust bearing 
as indicated, could be spun by a small permanent 
magnet and its porthole P brought into or out of 
register with the tube connecting vaporizing and 
measuring bulbs. During the preliminary fusions 
and vaporizations of barium and zine both 
shields were kept closed; a shield was opened 
only for the deposition of a film on the target T 
and both shields were never opened simul- 


taneously. In addition to their primary functions, 
the shields furnished an invaluable means of 
observing the progress of the vaporizations after 
the inner walls of the vaporizing bulbs had 
become coated with opaque layers of distilled 
metal. Slight rotation of a shield brought its 
porthole opposite an uncoated section .of the 
outer bulb. Burning the heater of the first 
vaporizer at low current gave adequate illumina- 
tion and approximate regulation of the rate of 
vaporization could be obtained by observation 
of the rate of deposition on the clear section of 
bulb thus exposed. 

_As in the tubes previously described, the 
measured films were laid down on a target which 
could be moved into position to receive any 
given fraction of either metal and then into the 
measuring position before the electron gun. The 
targets TT were, in this tube, duplicate Pyrex 
disks 18 mm in diameter, carefully mounted so 
that one reproduced the position of the other 
before the gun. A modification of the outgassing 
procedure applied in our previous work was 
necessitated by the relatively high volatility of 
zinc. The preliminary outgassing at 500—550°C 
was followed, after filling with dried air and 
introduction of the zinc and barium charges, by 
baking at 150°C rather than at 350°. Before 
sealing the tube from the pumps both barium 
and zinc were melted down two or three times to 
remove absorbed gases and these fusions were 
repeated after seal-off until approximately two- 
thirds of the original charges had been vaporized. 
The second vaporizers were flashed after each 
fusion to remove the condensates. The next 
fractions, amounting to an estimated 10-20 
percent of the total charges, were utilized for 
the measurements and were distilled slowly to 
the second vaporizers. One of the targets was 
then coated with a thick mirror-like film of zinc 
by slow distillation from the second vaporizer 
and the potential setting for this surface deter- 
mined. A succession of zinc films was then 
formed on top of the original zinc surface and 
each measured to determine the reproducibility 
and constancy of the potential setting for zinc. 
In certain experiments designed to test the 
possible effects of gas contamination, and dis- 
cussed below, the films were deposited almost 
instantaneously by raising the voltage applied 
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to the second vaporizer. The zinc film was then 
brought before the second barium vaporizer and 
a succession of barium films laid down and 
measured. Thick films of zinc and barium were 
then formed alternately on this same target and 
the contact potential difference Ba—Zn thus 
checked repeatedly. The whole series of measure- 
ments was then repeated with the second target, 
starting with a thick film of barium in immediate 
contact with the glass target. This procedure was 
designed to show whether or not the work func- 
tion of either metal was influenced by changing 
the substrate from glass to the other metal of 
the pair. No such influence was, however, 
detected. 


MATERIALS 


The barium employed in these measurements 
was taken from the same stock used in our earlier 
work and previously described. The zinc was the 
“spectroscopically pure’’ metal prepared by the 
New Jersey Zinc Company and kindly furnished 
to me by Mr. H. M. Cyr, with a spectroscopic 
analysis showing total impurities under 0.001 
percent. 


THE CONTACT DIFFERENCE OF POTENTIAL 
BETWEEN BARIUM AND ZINC 


For all of the fifteen pairs of barium-zinc 
surfaces measured the contact differences of 
potential, mean value 1.76 v, were reproducible 
to less than +0.02 v. As indicated above, a 
measurement for a given “pair’’ of surfaces 
included in some instances the determination of 
potential settings for as many as five individual 
zinc surfaces and as many barium surfaces. 
Potential settings for the individual films of 
such a series were reproducible to +0.01 v for 
both zinc and barium and the observed contact 
P.D. between aging surfaces remained constant 
to within +0.02 v for periods of measurement 
extending from 8 seconds after deposition to 
several days. The measurements were the most 
satisfactory yet obtained for any pair of metals 
studied during the course of this program, in 
which barium has consistently shown a repro- 
ducibility limited only by the constancy of the 
measuring circuit but the second metal, e.g., 
silver, magnesium, has been less reproducible. 
The entire series of measurements on barium- 


zinc included observations with three tubes, one 
of which was measured after the electrical circuits 
were rebuilt to improve their manipulation. 
Results for the three tubes were indistinguishable 
within the limits of reproducibility for any one 
tube. 


THE EXTERNAL WorRK FUNCTION OF ZINC 


The work function of barium films prepared 
by the standard technique described above is 
now known with a high degree of probable 
accuracy to be 2.52 ev at room temperature. 
Combination of this value with our result of 
1.76+0.02 v for the contact difference of poten- 
tial between barium and zinc assigns a value of 
4.28+0.02 ev to the external work function of 
zinc surfaces prepared by thermal vaporization 
and condensation on glass at room temperature. 
This result is to be compared with such of the 
earlier photoelectric and contact potential in- 
vestigations of zinc as have been made under 
modern vacuum conditions. In Dillon’s photo- 
electric study of polycrystalline and monocrystal- 
line zinc® the massive specimens were prepared in 
vacuum outside of the photoelectric tube, trans- 
ferred quickly to the tube and subjected to 
prolonged pumping and irradiation with ultra- 
violet light. The irradiation produced progressive 
changes in the initial threshold and, finally, 
equilibrium values of 3.32 and 3.57 ev, respec- 
tively, for polycrystalline zinc and Zn (0001). 
Dillon tentatively attributed these values to 
gas-free surfaces but recognized the possibility 
that they might have been characteristic of gas- 
coated surfaces stable under the conditions of 
the experiment. In the later work of DeVoe,’ 
also done in Mendenhall’s laboratory, the photo- 
electrically measured surfaces were prepared by 
distillation in the measuring tube itself under 
especially rigorous experimental conditions. The 
observed thresholds were highly reproducible and 
and were constant over long periods of time. 
The Fowler threshold was found to be 4.24 ev at 
a residual pressure of 1.5X10-* mm. This is the 
only published value in good agreement with the 
present results. 

Since the present measurements were taken in 
sealed-off tubes it was not possible to test 
DeVoe’s interesting observations: (1) that the 


* J. H. Dillon, Phys. Rev. 38, 408 (1931). 
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work function of zinc remains unchanged as the 
residual air pressure is increased from 1.5 
to 3X10-’ mm, but then drops rapidly to 
3.5-3.6 ev as the pressure is increased to 10-* 
mm, and (2) that the work function of zinc is 
insensitive to pure nitrogen up to at least 10-* 
mm, the highest pressure tested. These observa- 
tions are, however, borne out by the measure- 
ments of Oatley,‘ who obtained a work function 
of 3.4 ev at a residual pressure of 10-* mm. An 
extreme sensitivity of the work function to 
residual air pressure at pressures in the neighbor- 
hood of 10-* mm, such as is indicated by DeVoe’s 
curve, would account also for the work function 
fluctuations observed by Oatley. 


EFFECTS OF GAS CONTAMINATION AT LOWEST 
RESIDUAL PRESSURES 


Is the work function given by the present 
measurements, 4.28 ev, characteristic of clean 
zinc or must this value be attributed to a con- 
taminated surface in equilibrium with residual 
gas? In work function determinations by both 
the photoelectric and contact potential methods 
it has been customary to utilize the most 
rigorous vacuum and outgassing techniques 
possible and then, provided only that the ob- 
served values are constant and reproducible, to 
assign these observed values to clean metal 
surfaces. But contaminating gas films must be 
built up to some degree even in the highest 
attainable vacuum and, after equilibrium be- 
tween adsorbed film and bulk gas has been 
established, reproducibility and constancy of 
work function is to be expected for a contami- 
nated as well as for a clean surface. The funda- 
mental importance of this difficulty is generally 
recognized but the problem which it presents 
has received little explicit attention. 

The experimental conditions which appear to 
be prerequisite to an attack upon the problem 
are! (1) the metal surface must be so prepared 


as to be initially free from gas, and (2) the first 


work function measurement must be made before 
a gas film, sufficiently concentrated to produce a 
measurable change in work function, has had 
time to build up. If these conditions are fulfilled 
and if further measurements show no drift in 
work function, we may conclude that the equi- 


PAUL A. ANDERSON 


librium gas film is either too dilute to affect the 
work function or is composed of gases to which 
the work function is insensitive. While the first 
of these conditions can hardly be realized in its 
ideal form it is probably approached when the 
metal surface is prepared by multiple distillation 
in a gettered vacuum and the final vaporization 
is not too slow. If the maximum time lapse 
allowable under condition (2) could be calculated 
for any given metal-gas system it would then be 
possible to set up definite criteria for the meas- 
urement of any given metal. Such a calculation, 
involving time rates of formation of adsorbed 
gas films, evidently can be made only by a 
kinetic method of the Langmuir type and the 
minimum data required will not, in general, be 
available. In certain cases however, as shown 
below, enough information can be obtained from 
work function measurements themselves to en- 
able us to get around this difficulty. 

In the present measurements no drift in work 
function could be detected for freshly deposited 
zinc surfaces within a series of measurements 
beginning 8 seconds and extending to several 
days after formation of the surface. To obtain 
an estimate of the amount and probable effects 
of the gas deposited in the 8-second interval we 
compare the behavior of zinc with that of a 
suitably chosen “‘test surface’’ which is initially 
freed from gas and, like zinc, measured in a tube 
which has been gettered by vaporized barium, 
i.e., under a residual gas of comparable pressure 
and composition. The test surface is so selected 
that (1) it is sensitive to the same component of 
the residual gas as is zinc, and (2) its adsorption 
mechanism for this active gas is already known. 
The behavior of the test surface then furnishes 
an index to the rate at which the active gas, of 
otherwise unknown concentration, is supplied 
to the zinc surface. Since DeVoe’s gas tests 
show that oxygen is the active contaminant for 
zinc, tungsten suggests itself as an ideal test 
surface. The reflection coefficient of tungsten for 
oxygen is known to be extremely small and the 
work function of tungsten is extraordinarily 
sensitive to traces of adsorbed oxygen. Now it 
has been found in numerous published’ and 
unpublished measurements in this laboratory 
that tungsten, cleaned by flashing to 2800°K 


7 Reference 1; first citation, 
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and then allowed to cool, invariably shows a 
drift of work function, but that, in well gettered 
tubes, this drift does not become measurable 
until several minutes after flashing. In tubes 
immersed in liquid air the drift may be postponed 
for as much as 20 minutes. The drift is toward 
higher work function and continues for several 
hours before equilibrium is attained. In the 
earliest stages of adsorption m, the number of 
adsorbed molecules per unit surface of adsorber, 
is too small to affect the time rate of increase in 


‘mand we may write 


dn=(1—a) 


where a is the reflection coefficient and yv is the 
number of oxygen molecules which strike unit 
area in unit time. The effect of the inert com- 
ponents of the residual gas upon a presumably 
will be of the same order of magnitude for both 
tungsten and zinc and in the gettered vacuum v 
will be approximately the same for both surfaces. 
In terms of this equation our observations on 
tungsten show that dn/dt is so small that n 
attains the extremely small value necessary for 
detection only after a time interval of the order 
of minutes. Turning to zinc we may conclude, 
then, that even if all the oxygen molecules which 
strike its surface in 8 seconds were to remain 
adsorbed their number would be insufficient to 
produce any significant change in work function 
within this time interval. From the absence of a 
time drift of the work function of zinc, and its 


occurrence for tungsten, we may infer, further- 
more, that the mean time of sojourn, 7, of an 
adsorbed oxygen molecule must be much smaller 
on zinc than on tungsten. For it follows from 
purely statistical considerations which are inde- 
pendent of adsorption mechanisms, that for any 
equilibrium film of adsorbed gas 


r=n/v. 


Even when » is small, 7 is so large for oxygen on 
tungsten that at equilibrium » approaches the 
total number of adsorbing positions available 
per unit area. The fact that zinc, although 
sensitive to oxygen, shows no drift at these small 
values of v can be accounted for only if m and r 
are so small that the equilibrium film when 
finally completed is too dilute to produce a 
measurable change in the work function. We 
have, then, reasonably convincing evidence that 
the work function of a zinc surface prepared and 
measured in a barium gettered vacuum is charac- 
teristic of the uncontaminated metal. It may be 
mentioned, finally, that from the point of view 
of residual gas effects the contact potential 
method possesses a unique advantage over 
emission methods for determining work func- 
tions. Since it measures the average work 
function of a surface, the contact potential is 
relatively insensitive to slight contaminations 
which, localized on emitting areas, may have a 
relatively serious effect upon the emission charac- 
teristics of the surface. 
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The Water Vapor Molecule 


Byron T. DARLING AND Davip M. DENNISON 
University of Michigan, Ann Arbor, Michigan 
(Received November 9, 1939) 


The problem of the vibration rotation spectrum of water 
vapor is treated by means of the theory of semi-rigid poly- 
atomic molecules developed by Wilson and Howard. The 
potential energy is expanded as a power series in the normal 
coordinates and involves three zeroth-order constants, six 
first-order and six second-order constants. The positions of 
the band centers are calculated and found to depend upon 
ten quantities, X;, Xix, and y which are functions of the 
potential constants. A new feature of the treatment is the 
recognition of a resonance interaction between certain of 
the overtone bands which arises from the near equality of 
v, and »3. Eighteen band centers are known experimentally. 
These serve to determine the X;, Xix, y and furnish eight 
self-consistency checks which are very adequately satisfied. 
There exists no discrepancy between the Raman and infra- 
red spectra as reported earlier. In order to obtain the geo- 
metric displacements corresponding to each normal co- 
ordinate it is necessary to examine the spectrum of D,O. 
This not only furnishes the required information but also 


allows two independent checks upon the theory both of 
which turn out to be in nearly perfect accord. The inter- 
action between vibration and rotation is considered and 
the effective moments of inertia are calculated. These are 
functions of the normal frequencies and of the first-order 
potential constants. It is shown that A=I¢—I4—Tp de- 
pends only upon the normal frequencies and hence may be 
computed at once. A comparison between the observed and 
predicted A yields a very satisfactory agreement. The 
analysis of the rotational structure made by Mecke is 
supplemented by taking account of the rotational stretch- 
ing. The resulting molecular constants fix the valence angle 
to be 104°31’ and the O-H distance to be 0.9580A. From 
the effective moments of inertia the first-order potential 
constants may be evaluated and these, together with Xy 
determine the second-order potential constants. It is now 
possible to compute the interaction constant y and a 
comparison with the observed y again results most satis- 
factorily. 


1. INTRODUCTION 


HE vibration rotation spectrum of water 
vapor is a very extensive one and has been 
the subject of a great many investigations. The 
first strong line of the rotation series must occur 
in the far infra-red at about 500u and the 
spectrum extends with very few gaps through the 
near infra-red, where the fundamental vibration 
frequencies lie, up to about the middle of the 
visible spectrum. It has been accurately mapped 
from 1354 to 5700A. The far infra-red lines, which 
correspond to changes in the rotational energy 
only, have been measured with great precision by 
Randall, Dennison, Ginsburg and Weber' and 
lead to a determination of the rotational energy 
levels of the molecule. 

It is well known that the water molecule, which 
has the form of an isosceles triangle, possesses 
three normal modes of vibration. These funda- 
mental oscillations are shown in Fig. 1. Each of 
them should be active both in the infra-red and in 
the Raman spectrum. In point of fact, however, 
only v2 and vs have been observed? in the infra- 


1H. M. Randall, D. M. Dennison, N. Ginsburg and L. 


R. Weber, Phys. Rev. 52, 160 (1937). 
?E. K. Plyler and W. W. Sleator, Phys. Rev. 37, 1493 


(1931). 


red. They lie near 1595 and 3755 cm™, re- 
spectively. These bands are both extremely in- 
tense and the third fundamental 1 is probably 
weaker and is masked by the absorption of 13. 
The principal feature of the Raman spectrum? 


‘is a strong line at about 3650 cm which can be 


ascribed to »:. In addition to the three funda- 
mentals, some fifteen harmonic and overtone 
bands have been measured. These extend through 
the near infra-red into the vidible region of the 
spectrum where they are observed as atmospheric 
absorption lines in the solar spectrum. The ac- 
curacy with which these latter bands have been 
mapped is very high indeed although not es- 
sentially better than that attained in measuring 
the far infra-red rotation lines. 

The fine structure of each band, whether 
fundamental or overtone, is quite extensive since 
the moments of inertia of the molecule are small 
and the rotational energies large. The structure 
appears to be very irregular, as must be expected 
since no two of the moments of inertia are equal 
and the molecule belongs to the class of asym- 


3D. H. Rank, K. D. Larsen and E. R. Bordner, J. Chem. 
fies) 2, 464 (1934); D. Bender, Phys. Rev. 47, 252 
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metric rotators. Mecke‘ and his co-workers have 
made a careful study of the rotational structure 
of all the known bands and have succeeded in 
analyzing them in the light of the theory of the 
asymmetric rotator. The results of their work 
may be stated as follows. Neglecting a rotational 
stretching effect which will be discussed later, the 
rotational levels of the molecule in any particular 
vibrational state may be organized by means of 
three moments of inertia. These moments, which 
will be called effective moments of inertia; vary 
slightly with the vibrational state in question. 
Mecke showed that they might be expressed as 
constants plus small corrections which are linear 
in the vibrational quantum numbers. By this 
means he extrapolated to the moments of inertia 
of the nonvibrating molecule and obtained the 
molecular dimensions. He found the valence 
angle to be 104° 36’ and the O-H distance to be 
0.9558A. A characteristic feature of the effective 
moments of inertia is that they do not satisfy the 
relationship J4+Js=I¢ which must be expected 
to hold for any coplanar molecule. In a very 
important contribution to the theory of semi- 
rigid rotators Wilson and Howard® have shown 
that the rotational energy levels of a molecule 
should exhibit just the features observed by 
Mecke, namely that they may be correlated by 
means of effective moments of inertia which will 
vary with the vibrational state. 

An analysis of the positions of the centers of 
the fundamental and overtone bands has been 
made by Bonner® who used the conventional 
expression for the vibrational energy, 


W=Wot LX nit 22x 
= 


where ,, m2 and m; are the vibrational quantum 
numbers. The nine constants X;, Xu, may be 
determined from the eighteen observed bands 
and thus allow nine self-consistency checks. 
These checks are moderately good with one 
important exception. The fundamental frequency 
v1, measured as a Raman line, is predicted to lie 
from forty to fifty wave numbers below its 
observed position. 

*R. Mecke, W. Baumann and K. Freudenberg, Zeits. f. 


Physik 81, 313, 445, 465 (1933). 
asa Wilson and J. B. Howard, J. Chem. Phys. 4, 260 


*L. G. Bonner, Phys. Rev. 46, 458 (1934). 


? 
V, = 3650 cm" 


? 


e- 


Vs = 3755¢em" 


Fic. 1. Fundamental oscillations of the water molecule. 


The solution of this difficulty appears to lie in 
the fact that the two fundamentals »; and v3 have 
nearly the same value, thus paving the way for 
resonance effects. These effects are not of first 
order, however, since »; and v3 have different 
symmetry characters and consequently may not 
interact directly. It is well known that the 
vibrational energy levels of the nonlinear, sym- 
metrical molecule YX_2 belong to either of two 
symmetry classes. One of these is even upon 
reflection in a plane perpendicular to the X—X 
line and through the Y atom, and it is charac- 
terized by m; being equal to an even integer. The 
other class is odd upon reflection and occurs 
when ns is an odd integer. In designating a level 
it will be convenient to employ the symbol’ 
(m3,%1,M2). It is thus obvious that the levels 
which may interact due to the near equality of 
vy, and v3 will be of the type (m3,%1,%2) and 
(n3s+2,n,:—2,m2). The hypothesis that a reso- 
nance of this sort constitutes an essential feature 
of the analysis receives immediate confirmation 
when we examine the overtone bands which 
occur with greatest intensity. A glance at Table I 
shows that the strong overtone bands usually 
appear in pairs, thus (120), (300); (121), (301); 
(310), (130), etc. The members of these pairs 
comprise just those levels which could interact 
by resonance. An analysis of the positions of the 


7 The order (m3,%,m) corresponds with the notation used 
by Mecke which it seems wise to retain here. 
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TABLE I. Observed and calculated values of W/he. 


(nsnin2) W/hc(obs.) W/he(calc.) 
(001) 1595.0 _ 
(002) 3151.0 _ 
(010) 3650 _ 
(100) 3755.41 
(101) 5330.8 5330.6 
tt 10) 7251.0 7250.4 
111) 8807.05 8805.52 
(120) 10613.12 _ 
(300) 11032.36 _ 
Col 12151.22 12148.46 
301) 12565.01 12567.74 
(130) 13830.92 
(310) 14318.77 _ 
oH 15347.91 15346.27 
311) 15832.47 15834.12 
(132) 16821.61 16822.70 
(140) 16899.01 16894.25 
(320) 17495.48 17482.60 


energy levels will require in addition to the nine 
X; and X% of the usual theory, one further 
constant y which measures the effect of the 
resonance. It will be shown that, upon de- 
termining the ten constants with ten of the 
observed levels, the remaining eight levels are 
satisfactorily predicted by the theory. 

The general program which we hope to carry 
out in the present paper is as follows: The 
Hamiltonian function will be set down and the 
potential developed as a power series in the 
normal coordinates. The zeroth-order potential 
contains three constants, namely, the normal 
frequencies. The first-order potential consists of 
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2. THE HAMILTONIAN 


cubic combinations of the coordinates and con- 
tains six constants. The second-order potential 
which involves quartic combinations possesses 
also six arbitrary constants. A solution is made 
which yields the energy values of the system and 
is a function both of the vibrational and of the 
rotational quantum numbers. The band centers 
are obtained by setting the rotational numbers 
equal to zero and in this way the ten constants 
Xi, Xu, and y which were previously determined 
from the observed band centers are found as 
functions of the fifteen potential constants. To 
proceed further we turn to the rotational energy 
of the system and calculate the effective moments 
of inertia of the molecule. These are found to 
depend linearly upon the vibrational quantum 
numbers in agreement with Mecke’s observations, 
They are functions of the six cubic constants but 
not of the quartic constants. The effective mo- 
ments of inertia as determined experimentally 
furnish us with nine data, the coefficients of 
1, Ne and nz for each of the three moments, in 
addition to the equilibrium values of the mo- 
ments. These nine data fix the six cubic constants 
and allow three independent checks of the theory, 
which are made by comparing the observed and 
predicted values of A=Ic¢—I,—Izp. The agree- 
ment with the observations is most satisfactory. 
We now possess sixteen data, the six cubic 
constants from the effective moments of inertia 
and the ten X;, X% and y with which to calculate 
the fifteen potential constants. This implies the 
existence of one further test of the self-consistency 
of the theory and upon calculation it proves to be 
adequately fulfilled. 


In recent years a number of investigators have attacked the problem of finding the Hamiltonian 


of a semi-rigid, rotating polyatomic molecule. Of these, the formulation which seems to be most 
convenient in application is that by Wilson and Howard.5 A second derivation of this Hamiltonian by 
a somewhat shorter and more direct method has been made by one of us (B. T. D.) and will be 
published shortly. The result of this calculation is to obtain a Hamiltonian which differs only slightly 
and almost unessentially from the Wilson and Howard function, but which is Hermitian while it may 


be shown that theirs is not. 


Our Hamiltonian, when expressed in the Wilson and Howard notation, is: 


H= *(Ps— +3 + V. 
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The P. are components of the total angular momentum and p, are the components of the local 
angular momentum in the rotating coordinate system defined by Wilson and Howard. ses are the 
cofactors of a determinant yu which contains the moments of inertia. These quantities depend upon the 
normal coordinates and will be calculated in detail later. p, are the momenta conjugate to the normal 
coordinates gx. 

In the case of the water molecule, some simplifications occur. It will be shown that, of the local 
angular momenta, the only nonvanishing component will be perpendicular to the plane of the 
molecule. We denote it by ps. The potential function V we develop as a power series in the normal 
coordinates where the zeroth-order, first-order and second-order portions will contain terms which 
are, respectively, quadratic, cubic and quartic in the normal coordinates. Since the molecule is 
assumed to be symmetrical with respect to a plane perpendicular to the line joining the equilibrium 
positions of the hydrogen atoms and passing through the oxygen, the potential must obviously be an 
even function of gs, where gs3 describes the antisymmetrical motion v3 (see Fig. 1). The Hamiltonian, 
when written out, thus becomes: 


3 3 


+43 + +A3gs”) +3 (2191? +4292? +459193° 


The constants dj, Az, As, a6, constitute the fifteen arbitrary potential constants which 
were mentioned in the Introduction.® . 


3. THE VIBRATIONAL ENERGY 


In order to obtain an expression for the band centers we have merely to set the total angular 
momenta components P, equal to zero and examine the remainder of the Hamiltonian. It will prove 
convenient at this point to use, instead of the normal coordinates g;, the dimensionless variables x;. 
The x; are the arguments of the Hermitian orthogonal functions describing the normal vibrations and 
are related to g; as follows: x;=2 [wic/h }'gi. The w; appearing here and throughout the remainder of 
the paper are the normal frequencies expressed in waves per cm; c is the velocity of light. The po- 
tential energy in these coordinates becomes: 


V/he=4(w1x1? + + wars”) + + + 
+ (Bix14+ Bux + 3? + 


The relationship between the a;, 8; and the aj, 5; of the earlier formula may be easily obtained when 
needed. 

The solution of this part of the problem which describes the vibrating but nonrotating molecule 
may be readily obtained by the usual perturbation methods. One finds 


*It should be pointed out that the quartic part of the potential may contain three additional terms, namely ¢:@292*, 

gras? and q:*q2. These terms, since they are odd in at least one coordinate, will not contribute to the-energy. They would, 

ever, have had a slight influence in transforming from the Hamiltonian for H,O to that for D,O, since the normal 
coordinates of the latter molecule are not exactly equal to those of the former. 
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_ The constants X; and Xx, are the following functions of the potential constants: 


15a2? ay? a4 


as” a;* as? a6 


a3" a3” a, a4 


X12= K.E.P.+84— 


3aja5, as? 


we 2(2wstwi) 2(2w3—w;) 


ae? ae” 
we 2(2wst+we) 2(2w3—we) 


A calculation similar to ours was first made by Bonner.® Our results for that part of the X x. which 
depends upon the potential constants is in agreement with his when allowance is made for the 
difference in notation. 

The symbol K.E.P. which appears in the above formulas for the X represents the contribution 
from the kinetic energy part of the Hamiltonian. An examination of the Hamiltonian Eq. (1), shows 
that the K.E.P. may come from a number of terms, many of which turn out to be quite negligible. 
Thus the difference between $>-u'p.u-*p,u? and the zeroth-order expression 35° p,2 yields nothing to 
the order of approximation to which we are working. Similarly }u'psussu*psu' may be replaced by 
4ussps*. This latter quantity which involves the single nonvanishing internal angular momentum 3, 
will in general contribute both to X13 and X23. A calculation shows’ that for the water molecule the 
K.E.P. of X13 is so small that it may be set equal to zero but that the K.E.P. of Xe3 is +26.74 cm™. 

The final step in discussing the formula for the positions of the band centers is to consider the 
influence of the resonance interaction between levels of the type (m3m:m2) and (m3+2,m,—2,m2). The 
matrix elements in the Hamiltonian connecting these states are all off diagonal and they are all of 
second order. Thus the calculation of the X x, will be unaffected, since these depend upon second-order 
elements on the diagonal but only upon the first-order off-diagonal elements. A second way of stating 
this is that none of the terms in X%, contain the resonace denominator (w;—ws3). The method of de- 
ducing the final energy levels is now to diagonalize separately those portions of the Hamiltonian 


® The calculation is easily made when p; is explicitly determined as will be done in the next section. us; may be given its 
equilibrium value 1/J.°. The physical reason why X; is so little affected is that in water vapor the motions of the hydrogen 
atoms in the vibration »; are very nearly along the legs of the isosceles aay Thus the compound motion of »; and » 
together almost no angular momentum. For this reason the K. E. P. of X23 is just equal to the corresponding 
contribution to X23 in CO, which has heen'calculated by Weinberg and Eckart, J. Chem. Phys. 5, 517 (1937) and found to 


be 


3 15a;? a;* as” a3” 
2 4 w, 2we 8(2witwe) 8(2w1—we) 
| 3 
X 2= K.E.P.+-82: -—— — 
2 4we 8(2wetwi) 8(2w2—w1) 
as? 
2 2w1 2we 8(2w3—wi) 8(2w3t+we) 8(2w3—we) 
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matrix which contain the elements in question. The off-diagonal terms of these portions of the 
Hamiltonian consist of an interaction b. A short calculation shows that 


where 


as? 


306 


y=38s——-+ 


(3) 


It is seen that y depends only upon the potential constants. There exists a contribution from the 
kinetic energy part of the Hamiltonian but this turns out to be negligible in the case of the water 


molecule. 


As an example, let us consider the interacting levels (120) and (300). They will be given as the roots 


of the determinant 


Xi t+4X 024 2X 12 — W/he 


Vv 3y 


The roots may be expressed as W/hce+}[50?+6,2}! 
where W/hc is the mean value of the levels, in 
this case, is the 
separation of the levels without interaction, 
namely — 2X, 4+2X2—8X 114+ 4X and 6; 
is the separation caused by the interaction alone, 
in the present example 2,/3y. The highest fre- 
quency bands observed are (140) and (320). 
These are evidently members of a triplet whose 
third member is (500) and in this case the 
determinant which fixes the energy levels is of 
third order. 

The observed positions of the band centers are 
collected in the Table I. It will be noticed that 
the numbers given for the bands lying in the near 
infra-red differ slightly from the values usually 
appearing in the literature. This arises from the 
fact that we have corrected the wave-lengths to 
vacuum.!° Several -measurements* have been 
made of the frequency »; which appears in the 
Raman spectrum and the values which have been 
reported range from 3648 to 3654 
This discrepancy probably has its origin partly in 
the broadness of the Raman line (about 5 cm~'). 
We have adopted the even value 3650 cm-'. In 
order to calculate the ten constants X;, X % and 
y, we need ten data and these have been selected 
as indicated in the third column of Table I. In 
the case of (121), (301) and (131), (311) the 
individual levels were not used but only their 
sums and this is denoted by bracketing their 
calculated values. 


1 M. Rusch, Ann. d. Physik 70, 373 (1923). 


3X1 +9X11— W/he 


Since there are eighteen bands which have been 
observed, there will occur eight checks upon the 
theory. The agreement between the measured 
and calculated values is very satisfactory and is 
definitely much better than that obtained by 
Bonner. We feel that the small differences which 
do exist are due to the fact that we are ap- 
proximating the potential function by a power 
series development which is broken off after the 
quartic terms. This view is strengthened by the 
circumstance that the agreement appears to 
become poorer for the higher levels. The nu- 
merical values of the constants in waves per cm 
are as follows: 


X,=3693.89 Xi,:=—43.89 20.02 
X3= 1614.5 19.5 Xi3= — 155.06 
X3=3801.78 X33=—46.37 X23=— 19.81 

ly[= 74.46 w,=3825.32 

we = 1653.91 

W./hce=4631.25 w3= 3935.59. 


The zero-point energy Wp» has been calculated 
on the basis that the energy constant may be 
developed as a power series in n;+}, thus 


W/he=Ywi(ni+4) FLEX td) 


We are aware that the conventional second-order 
perturbation theory does not quite give this 
result as has been pointed out by King.'' We 
believe that this is due to the fact that the 
perturbation method involves a development 


"G, W. King, J. Chem. Phys. 5, 414 (1937). ‘ 
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which breaks off after a certain number of steps. 
A similar situation occurs in the theory of the 
diatomic molecule where the perturbation calcu- 
lation yields an energy which is not exactly a 
power series in (+4). In this case, however, it is 
possible to obtain exact solutions for a number of 
particular potential forms (i.e., the Morse po- 
tential) and in every instance the energy is 
expressed as a function of n+}. 

It will prove necessary for the rotation vibra- 
tion interaction which is to be treated in the next 
section that we know the actual geometric 
displacements corresponding to each of the 
normal coordinates. The following geometric” 
displacement coordinates will be used. Let us 
consider the Cartesian axes x and y whose origin 
lies at the center of gravity of the system and 
whose x axis is parallel to the line joining the 
equilibrium positions of the two X atoms. We 
now choose the x component of the change 
(increase) of distance between X atoms to be the 
coordinate g. Let y be the y component of the 
displacement of the Y atom with respect to the 
center of gravity of the X atoms. The final 
coordinate x is illustrated in Fig. 2 and is equal to 
x:+60, where x; is the x component of the dis- 
placement of the Y atom against the center of 
gravity of the X atoms, @ is the angle between the 
line joining the X atoms and the «x axis, while } 
is the equilibrium distance between the Y atom 
and the line joining the X atoms. From the 
conservation of angular momentum it is easily 
shown that x=(1+,b?/2ma?)x;. Here uw is the 


1 A geometric coordinate is one which defines the con- 
— of the molecule without reference to the masses 

the atoms. The constants appearing in the potential 
function are thus the same for any isotopic molecule. 
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reduced mass 2mM/(2m+M) and 2a is the 
equilibrium distance between the X atoms. 

The kinetic energy of the vibrating system may 
be calculated in the usual manner and is 


The general potential function will contain only — 


four constants (since it must be even in the 
coordinate x) and may be written 


2V =ax?+ by*+cq?+ 2dyq. 


The normal frequencies are given in terms of 
the roots (A;=47*w,*c*) of a determinant and may 


be expressed as follows: 


AitA2=b/u+2c/m, 
As=a(1+yb?/2ma*) 


The first minors of the determinant yield the 
relationships between the displacement coordi- 
nates and the normal coordinates qs. 


51191 + 41292, 
= 62191 + 
x= 63393, 
where 
(4) 


where the index i= 1, 2. 
[1/u+6?/2ma? }}. 
The following relations exist between the 6; 


1, 5217+ 520 =1/y, 

521511 4+522512=0, 


It is clear that the three normal frequencies of 
H.O will not serve to determine these relation- 
ships, since there are four potential constants 
involved. The normal frequencies of D,O will 
supply the remaining datum. 

The fundamental frequencies of D2O have been 
measured® although with somewhat less accu- 


% The frequency appears in spectrum and 
was reported b DH. Rank, . Larsen and E. R. 
Bordner, J. 2, 464 (i530). The bands and 
were observed by E. F. Barker and W. W. Sleator, J. 
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racy than the corresponding work for H,O and 
have the values »;= 2666, ve=1178.7, v3=2789 
cm-!. These frequencies are not the normal 
frequencies w, but, as may be seen from Eq. (2), 
are connected with them by the formulas, 
— 2X11 12 — 3X1, etc. We shall treat 
the problem of obtaining the X%. by a method 
- similar to that used by Bonner. In the diatomic 
molecule the coefficient of n? varies as the square 
of the frequency when changing from one 
isotopic molecule to another. We shall accordingly 
assume that varies as w? and X% as It 
will turn out that all three normal frequencies 
change by about the same factor and conse- 
quently we feel that this method, while not 
rigorously accurate, will lead to no serious errors. 
A preliminary solution gives the variation of 
each w. We find in this way that w;= 2757.86 and 
w2=1210.17.These numbers, which are derived 
from the observed fundamentals, may be tested 
by noticing that the ratio of wiwe for H,O to 
wwe for DzO must be equal to the square root of 
the ratio of mu for D.O and H,O, respectively. A 
calculation based on H=1.00813, D=2.01473, 
and O= 16.0000 shows that this relation is almost 
perfectly fulfilled; the observed ratio being 
1.89566, while the predicted ratio is 1.89539. Such 
an excellent agreement implies that our method 
of obtaining the X x for D2O has been a valid one. 
For purposes of computation w; is increased by 
0.20 cm! and we by 0.08 cm~ in order that the 
theoretical ratio of w,;w2 for HxO and D.O shall be 
exactly satisfied. 

The third fundamental will be treated slightly 
differently : w; of DxO may be obtained by using 
Bonner’s values for the molecular dimensions, 
a=0.75625A and b=0.58450A, together with the 
ratio of masses. We find w3= 2883.79. From this, 


Chem. Phys. 3, 660 (1935). Of these two frequencies v3 
is the less accurately known and was determined as the 
center of the zero branch associated with the band. We 
believe that the published value of 2784 cm™ is subject 
to the following corrections. In H,O, the center of the 
zero branch of v3 lies at about 3741 cm whereas the true 
band center as revealed by the rotational analysis is at 
3755.4 cm™'. This difference of 14.4 cm™ is caused by a 
type of convergence. Now in the case of the diatomic 
molecule the convergence factor varies as the cube of 
the normal —_ ay | and consequently we believe that 
the correction for the frequency v3 of DO should be 
approximately 5.7 cm. A reduction to vacuum wave- 
lengths yields the final result, »»=2789 cm. Fortunately 
the frequency v3 is not involved in the calculation of the 
potential constants, although it does furnish a test of the 
self-consistency of the theory. 
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together with the X x we calculate that v3= 2787.7, 
which is in very satisfactory accord with the 
observed number 2789. 

Collecting these results, we have for D,O: 


2758.06, Xi.=—22.81, Xi2=—10.56, 
we = 1210.25, Xo2= — 10.44, X13= — 81.92, 
w3= 2883.79, Xs33=—24.90, —10.62. 


The zero-point energy of D,O is then 3385.74 

The values of the potential constants and of 
the 6; for H:O now follow readily. 


a=10.672X10' dynes/em c=3.1344, 

b=7.1810, d= 3.1588, 
= +1.14678m-, 5:2 = +0.82758m-, 
+0.43910m-!, 522= —0.60845m-}, 
0.92828m-}, 


where m is the mass of a hydrogen atom, which 
we will take to be 1.673410- g. 

The sign of d is not determined directly but if 
we assume that the force field must be nearly of 
the valence type it follows that d is positive. 
We wish to remark that the accuracy with which 
the potential constants are known is certainly 
not as high as would be indicated by the number 
of significant figures given above. In making 
calculations, however, it is wise to make sure 
that no inconsistencies arise from rounding off 
numbers prematurely. 


4. THE VIBRATION ROTATION INTERACTION 


The main objective of the present section is 
to calculate the effective moments of inertia of 
the molecule. In a very recent paper Shaffer 
and Nielsen“ have also undertaken to solve this 
problem for the nonlinear symmetrical molecule 
YX2. They begin, as we do, with the Wilson and 
Howard Hamiltonian and carry through their 
computation along lines which in a general way 
parallel our work. Their final results for the 
moments of inertia are partly in agreement with 
ours but contain certain differences which appear 
to us to be vitally important in a discussion of 
the problem. Thus we find that the quantity 
A=I¢—I,4—T[sz is independent of the coefficients 
of the cubic terms in the potential, a result which 


“ W. H. Shaffer and H. H. Nielsen, Phys. Rev. 56, 188 
(1939). 
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is not attained by Shaffer and Nielsen.* We feel 1 
justified in giving a brief account of our calcula- et Oe His=ue3=0. 


tions, particularly since they appear to be very ; 
simple and straightforward. The Z, are obtained in the following manner. 


The moments and production of inertia in the The only component of the internal angular 
rotating coordinate system defined by Wilson momenta p, which does not vanish is ps3 and has 
and Howard may be readily expressed as func- the value 


tions of the displacement coordinates x, y and q = — — 
of the last section. Letting 1/f=(1+yb?/2ma’), ba= uf (xy — ty) — (ufb/2a) (xq — 29). 


we find: 


A =pb?+ py”, 


This may be expressed as a function of the normal 
coordinates by means of Eq. (3). 


bs=Ki(q3g1 — 9193) + 9293), 


C=A+B, where 
D= 2ubfxt+ Ki 0611/22), 
E=F=0. 


The sas which appear in the Hamiltonian Jt may readily be shown from the relations con- 
function are ; necting the dag that K,°+K,.?=1. A substitution 
B A D yields the result that for H,O,!° K,?=0.000019 

and K,?=0.999981. The quantity Z:, is the 


AB—D? AB-—D? AB—D* _ coefficient of in the formula for ps3. Thus 


Z3= — Kiqi—Koq2. 


Z:=Kigs, 


The was may now be expressed as functions of the normal coordinates. It is convenient to develop 
them as power series in the g; and to retain no terms beyond the second. 


Mir = 1/pb? —2(52191 + 52292) /ub*+ 3 + 3 (52191 + 
= 1/2ma?— (51191 + 51292) + 


M33 = uf?533"/2ma? — 


ma 
— f*u?5 334K 1K 


f533q3/ma2b — 3 f533q3(52191 + 522g2) /4ma*b? — 3 + 51292) /4ma*d. 


A feature of these expressions which will later turn out to be important is that the linear part of 
1/33 is equal to the linear part of (1/u11+1/y22). This arises directly from the fact that A+B=C. 

The Hamiltonian (Eq. (1)) may now be re-examined and each term evaluated to the desired order 
of approximation. The method which is employed is to diagonalize the Hamiltonian with respect to 
the vibrational part of the matrix elements, The resulting expression may then—in agreement with 


Wilson and Howard's predictions—be put in the form 
energy. 


The J4, Ig and I¢ are the effective amounts of inertia and are functions of the vibrational quantum 


* Note added in proof.—A re-examination of the article by Shaffer and Nielsen leads us to believe that the differ- 
ences may well be due to several typographical errors and that our results are essentially in agreement. 

46 The physical meaning of this very unequal division is interesting. K, and Kz represent the contributions to the 
internal angular momentum from the combined motion », v3 and v2, v3, respectively. In H,O the motion of the 
hydrogen atoms (whose amplitude is so much larger than that of the oxygen) is nearly along the legs of the isosceles 
triangle in both the motions »; and v3 but nearly perpendicular to the legs in v2. For this reason the combined 
motion »;, v3 will possess practically no internal angular momentum as compared with that of v2, vs. 
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numbers. It must be stated that the above equation is true only to second-order approximation ; 
higher orders will bring in quartic combinations of the P.. 

The first part }"yue.P.? contains, in addition to the constant pae?=1/J.°, terms which are linear 
in the g, and terms which are quadratic. These latter are immediately diagonalized, the q;* yielding 
h(nit+3)/42°wic while gig, gives no contribution. The linear part of uaa may be treated in either of 

two ways which lead to the same results. The first method recognizes that the presence of a linear 
term means that the centrifugal forces of rotation have displaced the equilibrium positions of the 
atoms. We may shift to the new origin by balancing the centrifugal force against the Hookes’ law 
force. This process will eliminate the linear part of wae and will bring in contributions from the cubic 
part of the potential. The second method is that of the standard second-order perturbation theory ; 
the terms linear in the g; and those which are cubic in the g; combine to produce an addition to the 
energy which may be readily calculated and which is equal to that obtained above. 

The second quantity in the Hamiltonian, }y12(P1P2+P:2P;), yields nothing in this order of approxi- 
mation. It may likewise be shown that the final term 


—}(psusstussps) Ps 


can be replaced by the simpler expression —y33°p3P3. 

Obviously —33°p3P contains no elements which are diagonal in the vibrational quantum numbers 
and therefore no first-order contribution to the energy. In second order it does, however, give an 
important term which may be calculated from the usual form 


(uss°)?P (WS—W,’). 


It should be mentioned that this term does contain a matrix element connecting the resonating 
states (m3m\m2) and (m3+2,n,—2,m2) but that its magnitude for H,O is negligibly small and need not 


be considered. 

The methods which have just been outlined allow us to obtain the effective yoo. The effective 
moments of inertia are found by inverting the wea and are given by the following expressions. We 
present them as functions of the potential constants rather than Let s*=h/*b’c and 
P=h/4r°a'c. 


I4=14°(1 «1! + (m1 +4) — 5214/1!) 
— f?533?/8ma*ws3+ 52105/w1!+ J, 
+351 101/013 + (m1 +}) (m2+}) 
— |, 
(ws? — w1?) (mi +4) + [hK — we") ](m2+3) 
—[(hK — w1*)) + (AK — w2”)) 


These formulas differ from those of Nielsen and Shaffer only in the parts containing the cubic 
constants a: --as. Our results, in contrast to theirs, show that A=I¢—I,4—Tz, is independent of the 
cubic constants. The reason why this comes about is easily understood. The cubic coefficients in the 
effective moments of inertia arise only from the terms in wag which are linear in the g;. These terms, 
as has been pointed out, satisfy the relation that 1/u33;=1/u11+1/u22 and consequently they can 
contribute nothing to the A. 

The quantity A plays an important role in the theory of the water vapor spectrum since it may be 
calculated from a knowledge of the normal frequencies and the shape of the molecule only and hence 
compared directly with the experimental data. Explicitly, 


A=(0.1644(m2+4) —0.0122(m3+43) ]X10-*. (5) 
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The coefficient of (w,+ 4) in the expression for A is practically zero since the numerical value of K,2 


is so small. 
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It will prove convenient to substitute the numerical values for the normal frequencies of the 
water molecule into the formulas for the effective moments of inertia in order to facilitate a correlation 


with the experimental data. We thus obtain 


(14 —I4°) X10 = (—0.0151 — 0.0001 127a1+0.000183 (21+4) + (—0.0672+0.0005493a2 
— 0.000037 (m2+3) + (—0.0149 — 0.000037 6a;+0.000183 
(In X10 = (— 0.0290 — 0.0002136a1 —0.0001807as) (1 +4) + (—0.0349 — 0.000542 
— 0.00007 (m2+4) + (—0.0280 — 0.00007 12a; —0.0001807.a5) (3+4). 


In the case of the water vapor molecule in which there exists a resonance between the states 


(m3,%1,M2) and (m3+2,m,;—2,mn2) there will occur a mixing of the effective moments of inertia for these 
states. If J,’ and J,? represent the moments of inertia as calculated from Eq. (6) and if I.’ and J,!! 
are the actual moments of the resonating states, then it may readily be shown, that 


I." = Tq (6+ 59) /26+1.2(6— 5o)/26, 
TM =I +59) 


where in these expressions 59 stands for the separation which would exist between the vibrational 
levels if there were no resonance interaction and 6=[6,?+6,"]! is the actual separation between 


the levels. 


5. CORRELATION WITH EXPERIMENT 


The final problem of this investigation will be 
to compare the experimental fine structure of 
the water bands with the theoretical formulas for 
the effective moments of inertia and thereby to 
obtain the molecular potential constants. Mecke* 
has analyzed the fine structure lines of each 
band and has shown that the rotational energy 
levels are just those which are predicted by the 
theory of the asymmetrical rotator. He thus 
calculates the effective moments of inertia for 
every particular vibrational state. We believe 
that his values may be somewhat improved by 
taking into account the rotational stretching of 
the molecule in the manner proposed by Wilson.'* 
It was stated in the last section that the applica- 
tion of a higher order perturbation to the 
Hamiltonian function would yield terms which 
are quartic in the components of total angular 
momenta P,. The coefficients are functions of 


TABLE II. Calculated values of —5W. 


20 


J 21 1, lo 
—6W(cm™) 0.331 0.032 0.045 0.015 0.006 0.047 — 


16 E. B. Wilson, J. Chem. Phys. 4, 526 (1936). 


(7) 


the normal frequencies only and may be easily 
evaluated. The final results of our calculation as 
given in Table II for a limited number of rota- 
tional states show that each level will be de- 
creased by what amounts to a small correction.” 

The effective moments of inertia are then 
obtained from the relations 


h/2n*cI 4 =(21—1_1) exp +0.28, 
h/2n*cI p= (20—1o)exp +0.03, 
Liexp+0.03 


and are shown in the second, fourth and sixth 
columns of Table III in units of 10-* g cm’. 
(We take h=6.624 X10-*7.)* 

The observed values of A=I¢—IJ4—TIzp are 
given in the seventh column, next to which 
appear the A as calculated from Eq. (5). The 
agreement is very satisfactory and is best for 
those levels which lie in regions most favorable 
for experimental determination. Thus the rota- 
~ 47 The work of Randall, Dennison, Ginsburg and Weber 
indicates that for the high levels the rotational stretching 
may be very large (for the state 11,, it is about 280 cm™); 
in fact, so large that it could not be simply computed. 

* Note added in ee he physical constants appearing 
in this paper are taken from a private communication 
kindly furnished by Professor R. T. Birge. 
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tional energy levels of the ground state (000) are 
undoubtedly known with the greatest accuracy. 
The strong bands (120) to (310) which appear 
as atmospheric lines in the solar spectrum have 
also been measured with high precision. The 
weaker solar bands, as well as the near infra-red 
bands, have not been as exactly determined. We 
regard the agreement between Aobs. and Acalc. as 
being a real and substantial test of the theory 
and as constituting one of the more important 
points of this investigation. 

The theory also predicts that the effective 
moments of inertia should be linear functions of 
the vibrational quantum numbers. A study of 
the data reveals that this is essentially the case, 
although small discrepancies exist which are 
usually of the order of the experimental errors. 
We have made the following treatment: The 
seven most accurately observed bands (000) and 
(120) to (310) are selected and their moments 
of inertia are ‘fitted to a linear formula by 
means of a least-square solution which is subject 
to two supplementary conditions; first that 
Ic¢°’=Ia°+Ip° and second that the A is just 
equal to the calculated A. Of course here, as in 
every case involving the resonating levels, it is 
necessary to take account of the mixing of the 
effective moments of inertia caused by the 
resonance. This is easily done by means of Eqs. 
(7) and we arrive at the result, 


I4 = 1.0229+ (0.0213) (m1 +3) 
— (0.1010) (m2+4) + (0.0486) 


Ig X10" = 1.9207 + (0.0398) (m1 +4) 
— (0.0249) 


TeX 10* = 2.9436 + (0.0611) (m1 +4) 
+ (0.0385) (m2+3)+ (0.0441) (ns+3). 


The values for the effective moments of inertia 
as calculated from the formula just written are 
given in the third, fifth and seventh columns of 
Table III. The agreement with the observed 
values is satisfactorily good and, as in the case 
of the A, is best for those bands which are most 
accurately known experimentally, thus leading 
to the belief that the discrepancies may be due 
to observational errors. 

The equilibrium values for the moments of 


TABLE III. Observed and calculated values of the effective 
moments of inertia and of A. 


|I 4 ops. I4 pops. Ip ops. catc.|A ons. A calc. 
(000) | 1.007 1.007 | 1.929 1.932 | 3.013 3.016 (0.077 0.076 
(001) | 0.904 0.906 | 1.899 1.907 | 3.107 3.055 | .304 241 
(100) | 1.057 1.056 | 1.928 1.939 | 3. 3.059 | .063 .064 
(101) | 0.988 0.955 | 1.923 1.915 | 3.133 3.009 | .222 .228 
(110) | 1.000 1.076 | 1.995 1.979 | 3.117 3.120 | .122 064 
(111) | 0.977 0.976 | 1.896 1.954 | 3.200 3.160 | .327 .228 
(120) | 1.110 1.105 | 2.014 2.012 | 3.185 3.178 | .061 061 
(300) | 1.141 1.147 | 1.969 1.962 | 3.155 3.151 | .045 042 
(121) | 1.000 1.003 | 1.978 1.987 | 3.240 3.217 | .262 226 
(301) | 1.046 1.046 | 1.940 1.947 | 3.201 3.190 | .215 207 
(130) | 1.133 1.136 | 2.040 2.040 | 3.244 3.233 | .071 057 
(310) | 1.161 1.158 | 2.006 2.014 | 3.217 3.219 | .050 047 
(131) | 1.030 1.034 | 2.003 2.015 | 3.287 3.272 | .254 .221 
(132) | 0.961 0.935 | 1.998 1.990 | 3.352 3.310 | .393 
(140) | 1.107 1.224 | 2.056 2.064 | 3.283 3.286 | .120 —.004 
(320) | 1.183 1.231 | 2.006 2.055 | 3.365 3.282 | .176 —.006 


inertia give the O—H distance to be 0.9580A 
and the valence angle to be 104° 31’; figures 
which differ only slightly from those obtained by 
Mecke* and Bonner.* A comparison of these 
formulas with Eq. (6) allows us to determine 
the six cubic constants. We find 


o=—322cm—, aa=+1, —909, 


ay=+216, +1060. 


—47, 


Knowing the cubic constants together with the 
numerical values of the six Xx, it is of course 
possible to obtain the six quartic constants. 


Bs=+35, Bs=+212, 
B2= +2, Bs=—116, Be= —108. 


It is interesting to notice that the perturbing 
terms in the potential which connect the levels 
m, and m3, namely x,x;* and x,°x;?, appear to be 
unusually large. This no doubt accounts both 
for the resonance between the levels (m3m2), 
(ms+2,m,—2,m2) and for the fact that they 
appear to be the most intense overtones in the 
spectrum. 

Since we now possess all twelve of the anhar- 
monic potential constants, together with the 
normal frequencies, we may calculate the inter- 
action constant y by means of Eq. (3). We find 
the value Yeaie= —71.5 cm~ which may be com- 
pared with the experimental figure of | Yous! 
=74.46 cm—!. The agreement, within about 4 
percent, is remarkably good and argues strongly 
for the essential correctness of our interpretation 
of the water vapor spectrum. 
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About eighty lines belonging to the triplet systems of the magnesium-like ions, Cl VI, K VIII, 
Ca IX, and Sc X, have been identified. Term values for each ion are reported with respect to 


N two previous papers from this laboratory 

studies of the sodium-like spectra of Cl' and 
K, Ca, and Sc? were reported. It is the purpose 
of this paper to present data on the magnesium- 
like spectra of these same elements. 

The spectra of the first five ions of the mag- 
nesium sequence have been known for some 
time.*: * Bowen and Millikan’ reported the wave- 
lengths of the pp’ group in Cl VI without 


TABLE I. Classified lines of Cl VI. 


Int A(VAC.) CLASSIFICATION 
3* 736.762 135,729 | 3s3p*P,—3p* *P, 
3* 733.891 136,260 
3 3 
4* | 730.311 | 136,928 { = 
3* 727.537 137,450 
3* 724.129 138,097 
3 555.580 179,992 | 3s3p *P.—3s3d 8D, 
20 555.485 180,023 
-2 552.053 181,142 3D, 
10 551.992 181,162 3D, 
5 550.355 181,701 3Py 3D, 
8 399.995 250,003 | 3s3d*D,—3s4f *F, 
7 399.957 250,027 3D, 3F; 
5 399.938 250,039 3D, 3F, 
25 325.161 307,540 | *P.—3s4s 8S, 
20 323.936 308,703 
15 323.356 309,257 3S, 
3 243.883 410,033 | 3s3p *P2.—3s4d 8D, 
12 243.854 410,082 3P, 3D; 
2 243.208 411,170 ap, 3D, 
8 243.194 411,194 3D, 
5 242.885 411,718 2D, 
3 195.227 | 512,224 | 3s3p*P.—3s5d 5D; 
2 194.796 513,358 3D, 


* Intensity as reported by Bowen and Millikan. 


* Now at the University of Chicago. 
1L. W. Phillips, Phys. Rev. 53, 248 (1938). 
(1939) G. Kruger and L. W. Phillips, Phys. Rev. 55, 352 
*R. F. Bacher and S. Goudsmit, Atomic Energy States 
(McGraw-Hill Book Co., 1932). 
‘I. S. Bowen, Phys. Rev. 39, 8 (1932). 
51. S. Bowen and R. A. Millikan, Phys. Rev. 25, 597 


(1925). 
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attempting to make J value assignments. More 
recently Edlén® has reported on the magnesium- 
like spectra of the ions Ti XI to Co XVI. The 
data reported here serve to bridge the gap in 
the sequence from Cl VI to Sc X with the excep- 
tion of A VII. A preliminary attempt in this 
laboratory at the excitation of the intermediate 
stages of ionization of argon has shown some 
indications of success, and it is hoped that these 
data will be available soon. 

Classified lines in the spectra of Cl, K, Ca, and 
Sc are listed in Tables I-IV. In all cases the 
position of the multiplets was first predicted by 
means of a constant-first-difference or constant- 
second-difference graph, and the separations were 
checked by the (Av)! law. Since the 3s3p *P 
term in Cl VI has been identified by its combina- 


TABLE II. Classified lines of K VIII. 


INT A(VAC.) v(cm™) CLASSIFICATION 
5 572.790 174,584 | 
4 569.479 175,599 3Py 
6 565.112 176,956 
12 564.462 177,160 3P, 
6 561.593 178,065 
2 557.029 179,524 3P, 
2 422.642 236,607 3s3p *P2—3s3d *D; 
10 422.511 236,680 3P, 3D; 
1 418.571 238,908 3D, 
7 418.443 238,981 3P, 3D, 
3 416.628 240,022 3D, 
15 230.745 433,379 3s3d *D;—3s4f *F, 
12 230.706 433,452 3D. 3F; 
10 230.678 433,504 3D, 3F, 
15 199.921 500,198 3s3p *P2—3s4s 4S; 
10 198.975 502,576 ap, 3$, 
5 198.536 503,688 
1 156.566 638,708 3s3p *P.—3s4d *D, 
3 156.543 638,802 3P, 3D, 
9 156.507 638,949 3P, 3D; 
4 155.985 641,087 3P, sD, 
7 155.961 641,185 3P, 3D, 
5 155.715 642,199 sD, 


6 B. Edlén, Zeits. f. Physik 103, 536 (1936). 
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TABLE III. Classified lines of Ca IX. 


K, Ca AND Sc 


TABLE IV. Classified lines of Sc X. 


A(vac.) v(cm~') CLASSIFICATION 


A(VAC.) v(cm™) CLASSIFICATION 


Zz 
a4 


515.576 
512.060 
507.105 
506.163" 
503.261 
497.996 


378.550 
378.382 
378.089 
373.979 
373.802 
371.899 


184.433 
184.377 
184.338 


163.229 
162.371 
161.979 


129.970 
129.934 
129.445 
129.425 
129.197 


100.958 
100.645 


193,958 3s3p *P2—3p? *P, 
195,290 3P, 
197,198 3P, 
197,565 3P, 
198,704 
200,805 5P, 


264,166 3s3p *P.—3s3d *D, 
264,283 5P, 8D, 
264,488 
267,395 
267,521 
268,890 


542,202 3s3d *F, 
542,366 3D, °F; 
542,483 sD, 


612,635 | 3s3p *P,—3s4s 3S; 
615,872 3S; 
617,364 3S; 


769,410 3s3p *P2—3s4d *D, 
769,620 5P, 3D; 
772,530 sD, 
772,650 8D, 
774,010 5Po sD, 


990,510 | 3s3p *P,—3s5d 
993,590 


* Identified by Edlen. 
Bl—Blend with another line. 


tions with several other terms, it is now possible 
to make correct J value assignments in the pp’ 
group. Since this group of lines is out of the 
range of the Illinois spectrograph, the intensities 
reported are those given by Bowen and Millikan.*® 

Whitford’ has classified six lines taken from 
Ekefors’* K data as the pp’ group in K VIII. 
The six lines he chose are the correct ones, but 
the J value assignments, and hence his term 
table, appear to be incorrect. The classification 
given in Table II is verified by other combina- 
tions with 3s3p *P term. 

Four of the lines in the pp’ group in Ca IX are 
given in Ekefors’® Ca data and were correctly 
identified by Edlén.* Two of the lines in the 
3s3p *P —3s3d*D multiplet in Ca IX appear to 
be more intense than one would expect them 
to be. Inspection of the plate indicates that in 
reality these two lines are both blends with 
other calcium lines belonging to some other 
stage of ionization, and it has been so indicated 
in Table III. 


7A. E. Whitford, Phys. Rev. 46, 793 (1934). 
*E. Ekefors, Zeits. f. Physik 71, 53 (1931). 


343.100 
342.877 
342.509 
338.115 
337.888 
335.815 


151.001 
150.941 
150.902 


135.921 
135.131 
134.766 


291,460 3s3p *P2—3s3d *D, 
291,650 ‘Dy 
291,963 
295,757 8D, 
295,956 sD, 
297,783 sD, 


662,247 3s3d *F, 
662,511 
662,682 3D, 


735,721 3s3p *P2—3s4s 
740,020 3S) 
742,025 4S, 


TABLE V. Term values. 


CIVI K VIII Ca IX 


0 0 
553 1,495 
1,718 4,735 


136,713 196,785 
137,449 198,698 
138,650 202,300 


181,698 268,890 
181,713 269,017 
181,741 269,223 


3s4s *S, | 309,257 617,367 


3s4d*D, | 411,721 774,017 
8D, | 411,749 774,145 
8D; | 411,800 774,355 


3s4f °F; 811,368 
°F; | 431,742 811,388 
°F, 811,420 


3s5d*Dz | 513,911 995,085 
3D; | 513,942 995,245 


297,777 
297,969 
298,280 


742,033 


Term values for deep-lying states of the 
triplet system are given in Table V for all four 
ions. It has not been possible to make positive 
identification of any lines belonging to the 
singlet system, although the transition 3s? 'S 
—3s3p'P has been tentatively identified in Cl, 
K, and Ca. In Cl and K the separations of the 
3s4f *F term calculated from observed lines turn 
out to be smaller than (or of the same order of 
magnitude as) the error of measurement. The *F 
level in these two ions has accordingly been 
reported as unresolved. 

In conclusion, the authors wish to express their 
appreciation to Professor P. G. Kruger for his 
interest in the problem. 
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In this paper, calculations are presented for the de- 
termination of the magnetic energy, assuming dipole-dipole 
interaction, of a macroscopic crystal possessing para- 
magnetic ions situated on a definite lattice structure. 
Parallel and antiparallel configurations are considered as 
well as specimens of different shapes and lattice structures 
of various types. It is shown that the configuration of 
lowest energy cannot be strictly defined but is dependent 
upon the shape of the specimen. This configuration is an 
antiparallel one for a sphere but changes to a parallel one 


for a long slender specimen. The transition between the two 
configurations appears to take place at a certain critical 
value of the demagnetization factor. For specimens above 
this value the theory accounts qualitatively for the low 
values of the magnetization found experimentally and for 
the apparent constancy of the Curie point with changes of 
shape. Below the critical value the theory implies a pure 
ferromagnetic state, which is not in agreement with 
experiment. Probable reasons for the discrepancy are given. 


N attempting to give a theoretical explanation 
of the anomalous experimental results of 
Kurti and Simon! on paramagnetic salts at low 
temperatures, the author found it desirable to 
first make a definite determination of the configu- 
ration of lowest energy. This paper is devoted to 
that end. The model used is a three-dimensional 
lattice of freely suspended magnetic dipoles, each 
of which is assumed to possess two orientations, 
either parallel or antiparallel to some given 
direction. 

To compute the configurational energy, the 
exact method of attack would be to calculate by 
wave-mechanical means the energy levels of the 
assembly as a whole, considering the crystal as 
one large molecule. Owing, however, to the 
complexity of the wave equation when inter- 
action between particles is considered, this 
method is prohibitive. Therefore we resort to a 
purely classical determination of the magnetic 
energy of the given three-dimensional con- 
figurations, assuming different possible arrange- 
ments of the paramagnetic ions. From spectro- 
scopic analysis, as carried out by Beevers and 
Lipson,” the structure of such paramagnetic salts 
as iron alum is now well known. The unit cell is 
about 12 angstroms in length and the metal 
atoms are placed on a face-centered structure. It 
is due to these comparatively large distances 


*The major part of this paper was written while the 
author was at Cambridge Deruniy, England. 
1 Kurti, Laine and Simon, Comptes rendus 204, 675, 754 


(1937). 
A. Beevers and H. Lipson, Proc. Roy. Soc. 148, 
664 (1938). 


apart that the exchange forces, which fall off 
exponentially with the distance, probably do not 
matter; whereas the magnetic forces, varying 
only as the inverse third power of the distance, 
become of prime importance. 

The dependence of the energy upon the mag- 
netic forces rather than the exchange forces 
between atoms makes the calculations much 
more complex as the magnetic interaction be- 
tween two magnets, considering only dipole 
moments, depends on their orientation relative to 
the connecting line as well as on their distance 
apart. The classical expression for the interaction 
energy is 


(ui ri) (ur (1) 


where yu; and y; are the strength of the magnetic 
dipoles and rj; is the distance between the 
two. If the dipoles have the same direction as 
well as strength, the expression simplifies to 
u’r*(1 —3 cos? 0) where @ is the angle between the 
dipole direction and the line joining the two 
dipoles. 

For a single row of dipoles the lowest energy 
is obtained when they are all parallel to one 
another and directed along the line joining them. 
For a double row or a two-dimensional set-up, 
this is no longer true. Here a direct evaluation of 
Eq. (1), summed over all dipoles, shows that the 
antiparallel alignment—one string pointing to 
the right, the next nearest string to the left, and 
so on gives the deepest energy. 

For the three-dimensional solid, with magnetic 
dipoles situated at all lattice points and all 
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pointing in the same direction, the magnetic 
energy is given by 
pe 212 —m?— n?2 


21? —m*—n?* 


where N is the number of ions per cc and /, m, and 
n are integers running from — © to . This is an 
infinite summation of the type 1/r* in three 
dimensions. It cannot be definitely evaluated as 
the series is logarithmically divergent. This is 
equivalent to saying that the value of the 
interaction energy will depend on the shape. For 
a spherical or cubic specimen, symmetry require- 
ments show directly that W will be zero. For a 
specimen of any other shape it becomes necessary 
to resort to the Lorentz approximation. Lorentz 
considers a small sphere surrounding the dipole in 
question and assumes that the summation over 
all dipoles can be expressed as a summation 
within the sphere and an integration outside. The 
summation contributes nothing to the field at 
the center, provided the dipoles are arranged in 
symmetrical fashion. There will, however be a 
contribution for unsymmetrical arrangements. 
This point has been considered by Sauer and 
Temperly,* who show that the internal field may 
be taken as proportional to the degree of order of 
the dipole arrangement. The value of the inte- 
gration over all the dipoles outside the small 
sphere is determined by considering the effect of 
the surface charge produced by these outside 
magnets. On the surface of the small sphere itself 
there exists a charge density which contributes to 
the field at the center a factor of (427/3)I, where J 
is the intensity of magnetization per unit volume; 
the additional field due to the surface charge on 
the external boundary of the specimen is brought 
in by what is commonly referred to as the 
demagnetization factor. For a sphere the demag- 
netization factor of — (42/3) exactly annuls the 
local field factor of + (42/3) so that the energy is 
given correctly by the Lorentz approximation. 
For parallel alignment of the elementary 
dipoles the magnetic interaction energy, using 


of | ‘. Sauer and H. N. V. Temperly, Proc. Roy. Soc. (in 
print). 
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the Lorentz approximation, can be obtained for 
any external shape of specimen provided the 
appropriate demagnetization factor is given or 
can be determined. But for other alignment 
configurations such as the antiparallel one, the 
Lorentz method fails and it becomes necessary, in 
order to get a reasonable approximation to the 
correct magnetic energy, to take a direct sum 
result over a large number of dipoles. All direct 
sum results hereinafter given will be determined 
by taking into account the effect of all dipoles 
within a small cube surrounding the dipole in 
question and containing approximately 125 
atoms. Careful attention must be paid to alge- 
braic signs in carrying out this summation as 
antiparallel dipoles give an interaction energy of 
opposite sign to that of parallel dipoles. For 
purposes of direct summation the equation for W 
can conveniently be written as 


W= (3) 


Here 7 is the right-angle distance from the center 
dipole to whichever string is being considered. 
Thus if we were determining the energy for an 
antiparallel arrangement in a simple cubic 
crystal, y would take the values 1, 4/5, --+ for 
the antiparallel dipoles, ./2, 2, 1/8, --- for the 
parallel dipoles, while W has the sign given by 
Eq. (3) for the parallel rows and changes sign for 
the antiparallel. 

Although the paramagnetic ions are located on 
a face-centered lattice, it appears worth while for 
reference purposes to include the energy calcu- 
lations for the simple cubic structure and the 
body-centered structure as well. In the case of 
the body- and face-centered structures, it is 
possible that a deeper energy might be obtained 
by assuming the magnets directed along the body 
or face diagonal rather than along the axis of the 
unit cell. Such indeed proves to be the case. We 
also distinguish between two different possible 
types of antiparallel alignment which we shall 
refer to as “A’’ and “B.” In arrangement “A” 
the rows closest to one another are alternately 
directed. For example, in the body-centered 
cubic lattice the corner dipoles would all point 
forward and the body-centered dipoles would 
point backward. Ini rangement “B,” the nearest 
dipoles situated on a plane perpendicular to the 
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Fic. 1. Showing sections taken at right angles to 
direction of dipoles. 


axis are alternately directed; in the example 
given above, this means that the corner dipoles 
are themselves arranged in antiparallel fashion 
and similarly for the body-centered ones. Any 
possible ambiguity arising from these definitions 
will be compensated by illustrating, wherever 
necessary, the difference between the two ar- 
rangements. The magnetic interaction energy for 
each of three alignments, parallel, antiparallel 
“A,” and antiparallel ‘‘B’’ will now be obtained 


for different shaped specimens. 


INTERACTION ENERGY FOR A SPHERICAL 
SPECIMEN 


(1) Parallel Alignment: The magnetic energy 
for a simple cubic lattice structure is found to be 
zero whether the determination is made directly 
from Eq. (2) or by means of the Lorentz method. 
This result still holds good for the body-centered 
and face-centered lattices, as they consist, re- 
spectively, of two and four interlocked cubic 
lattices. 

(2) Antiparallel ‘‘A’’ and “‘B’’: It is necessary 
here to consider each one of the three lattice 
structures separately. 

Simple Cubic Structure : Evaluation of Eq. (3), 
taking into account the interaction of the 
nearest 125 dipoles, leads to an energy value of 
W=—2.7N°y’. 

Body-centered Structure: We consider two 
different possibilities of dipole direction, along 
edge of unit cell or along body diagonal. For each 
of these chosen directions, the magnetic energy is 
calculated for both types of antiparallel align- 
ment. The two configurations ‘‘A’”’ and “B”’ are 
shown in Fig. 1, the circles indicating parallel 
dipoles and the crosses antiparallel ones. The 
interaction energy for the (001) direction and 
with arrangement “A’”’ is zero, while with 
arrangement “‘B,”’ the result is —1.35. For the 
diagonal direction, W is —1.75 for “‘A”’ and “‘B.” 


Face-centered Structure: Fig. 1 may again be 
used for noting the essential difference between 
the two types of alternating alignment. The 
approximate values of the energy are 2.2 for 
arrangement ‘‘A’’ and dipole direction along cube 
edge, and —1.1 for the same direction with “B.” 
With face-diagonal direction of the magnets, W 
decreases to —1.1 for “A’’ and —1.8 for “B.”’ 


INTERACTION ENERGY FOR AN ELLIPSOIDAL 
SPECIMEN 


(1) Parallel Alignment: For this case the energy 
may be expressed as 


(4) 


where by D.F. is meant the demagnetization 
factor. The values obtained will differ from the 
energies in the case of the sphere by having the 
proper demagnetizing factor, depending on the 
axis ratio of the ellipsoid, inserted in place of the 
4/3 of the sphere. For example, a specimen with 
an axis ratio of r=3 possesses a demagnetization 
factor of 1.37 and gives rise to an interaction 
energy of 1.9. As the axis ratio approaches 
infinity, D.F. approaches zero and W is then 
given by in 
units of N2u?. The results are the same for all 
three types of crystal lattice structure. 

(2) Alternating Arrangements: The energy 
values as given above for a spherical specimen 
are true here as well. This is so because in the 
alternating arrangement of dipoles there exists 
no charge density on the surface and hence the 
shape of the boundary does not matter. 

For simplicity of reference and comparison, the 
foregoing results, together with some other 
energy values as yet unmentioned, have been 
assembled in Table I. Some of these magnetic 
energy values have already been given by Van 
Vleck,‘ but as no details of the calculation and 
also no explicit values for the possible diagonal 
directions appeared in his paper, it was thought 
desirable to collect all the values in one Table. In 
addition to the energies for the parallel con- 
figuration, the antiparallel ‘“‘A,”’ and the anti- 
parallel “B,’’ Table I includes the magnetic 
energy for a rather special configuration, pictured 
in Fig. 2. This arrangement of dipoles might be 


4J. H. Van Vleck, J. Chem. Phys. 5, 556 (1937). 
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termed a nonmagnetic, disordered state of the 
crystal as compared to the parallel arrangement, 
which is a magnetic, ordered state and the 
antiparallel arrangement, which is also ordered 
but nonmagnetic. This particular disordered state 
is only one of a number of possible disordered 
arrays, yet it is of interest to compare its energy 
values with those of the ordered arrangements. 

Examination of Table I allows us to draw the 
following conclusions: 

(1) The parallel alignment gives an energy 
independent of the type of lattice structure but 
dependent upon the shape of the specimen. That 
is, for a given shape, the energy is the same 
whether the dipoles are arranged on a body- 
centered, face-centered, or simple cubic lattice 
yet a change in specimen shape from a sphere 
to a long slender specimen changes the energy 
from the value zero to —2.1. 

(2) Conversely, the antiparallel configurational 
energy is dependent upon the type of lattice 
structure but independent of the shape. For 
example, the energy changes in value from —2.7 
to 1.1 in going from a simple cubic lattice to a 
face-centered one, while a change in the external 
boundary leaves the values unaltered. Comparing 
the two types of antiparallel alignment, we see 
that the lower energy is in general associated 
with type ‘“B” rather than “A.” 

(3) Dipoles set along the direction of the face 
diagonal for the face-centered lattice and along 
the direction of the body-diagonal for the body- 
centered lattice, lead to deeper energy values 
than when directed along the cube axis. This is 
what we should expect for changing direction 
from cube axis to body or face diagonal reduces 
the distance between adjacent dipoles of a 
parallel string and hence causes an increase in the 
magnetic energy. 

(4) For parallel alignment dipoles tend to set 
themselves in such a direction that the demag- 
netizing factor is a minimum. In other words, in 
an ellipsoidal specimen, the dipoles would prefer 
to align along the long axis rather than the 
short. 
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TABLE I. Energy constants. All direct sum values (anti- 
parallel ones) given to within 0.05. a—simple cubic, dipole 
direction along (001); b—body-centered, dipole direction 
along (001) ; c—body-centered, dipole direction along (111); 
d—face-centered, dipole direction along (001); e—face- 
centered, dipole direction along (011). 


ENERGY 
CONSTANT 


LATTICE 
STRUCTURE 


ARRANGEMENTS 


Parallel a, b,c, d,e 


ean 


Disordered 


(5) The actual configuration of minimum 
energy is not strictly defined but is dependent 
upon the shape. For a sphere or an ellipse with 
low axis ratio, the antiparallel arrangement is 
definitely the lower. For a long slender specimen 
or an ellipse with high axis ratio (r>6 for the 
face-centered lattice) the parallel configuration 
becomes the more stable. 

(6) A comparison of the values from Table I 
for the parallel and the disordered configurations 
indicates that a spherical specimen, in the 
absence of any external field, has equal energies 
in the completely ordered state and in the 
completely disordered state. This, perhaps, may 
not be an expected relation; yet when it is 
recalled that for a spherical specimen the 
Lorentz and the demagnetizing terms cancel, 
thus causing the magnetic interaction energy to 
vanish, the correlation appears more natural. 
The recorded energy value for the disordered 
array was calculated only for a simple cubic 
lattice but it appears likely that the same result 
would be obtained for the other lattice types. 

The three-dimensional model studied in this 
paper appears to eliminate the chance of a 
spherical specimen becoming spontaneously mag- 
netized as a result of pure dipole-dipole coupling. 
The ferromagnetic state would be possible, 
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though, for a long slender specimen. Now it is 
well known that the usual Ising model, to which 
our model roughly corresponds, will never give 
ferromagnetism. This result, however, holds only 
in the case of a linear chain. Peierls’ has gener- 
alized Ising’s model to include the two- and three- 
dimensional cases and he finds that ferromag- 
netism is possible. It is not correct though to 
assume, in our model, that even a long slender 
specimen will necessarily be ferromagnetic. For 
even if the lowest energy level of the crystal does 
correspond to the parallel state, there may exist 
other nonparallel states of only slightly higher 
energy values; and the greater number of 
complexions of these nonparallel states might 
easily be sufficient to overcome the energy 
advantage of the parallel level. 

By taking the antiparallel configuration as the 
stable one at low temperatures, and by using the 
energy values given in Table I, it is possible to 
determine the Curie points. This is done in a 
paper by Sauer and Temperly* who show that the 
Curie point, for all specimens in which the 
antiparallel state is the stable one, i.e., for 
ellipsoidal shapes with r less than 6, the critical 
temperature is independent of the shape. This 
appears to be in confirmation with the experi- 
mental results as: Dr. Simon has reported that 
for a sphere as well as long samples, there is very 
little change in the Curie point. 

In view of our results, the Curie point for a 
sphere is now to be interpreted as the tempera- 
ture at which a transition takes place between a 
disordered state and the ordered antiparallel 
state. Heretofore, it has usually been taken as an 
indication of a conversion from the disordered 
state to the ordered parallel one ; because of this 
interpretation, one expected the measured mag- 
netization to be in the neighborhood of its 
saturation value. All experiments, however, show 
it to be very much smaller. This has always been 
considered one of the very puzzling features of 
the experimental results but the mystery is 
readily removed when we recognize that the 


* R. Peierls, Proc. Camb. Phil. Soc. 32, 471 (1936). 


ordered state below the Curie point is not 
magnetic but essentially nonmagnetic. 

For an ellipsoidal specimen with r>6 the 
reason for a low value of the magnetization is not 
so apparent, as here the Curie point should 
correspond as usual to the transition to the 
parallel saturated state, provided, of course, that 


our simple treatment is adequate. Before definite — 


conclusions can be drawn there is a great need for 
more experimental evidence, especially in the 
region of r >6. The low value of the magnetization 
is inferred, at present but from a single experi- 
mental value on iron alum with an r ratio of 6, 
Still there does appear to be a definite discrepancy 
between theory and experiment. The probable 
reason for this discrepancy is the too definite 


-postulation of the existence of the parallel state as 


the stable one on the simple evidence that its 
energy value is slightly less than that of the 
antiparallel state. In a rigorous treatment the 
effect of the greater statistical weight of the 
nonparallel states must needs be considered. 

It requires very little change to convert the 
energy constants of Table I to conform to the 


case of electrical dipole-dipole interaction instead 


of magnetic. Nothing in the electrical case 
corresponds to the demagnetization factor in the 
magnetic, hence the D.F. term of Eq. (4) must 
vanish. This converts all the different energy 
constants for parallel alignment to the single 
value —2.1. The antiparallel configurational 
energies remain the same. A comparison now 
leads to the conclusion that, for the simple cubic 
lattice, the antiparallel arrangement still gives 
the deeper energy. For the face-centered and 
body-centered lattices, however, the electric 
dipole coupling leads to different results than the 
magnetic. It indicates that, of the various 
arrangements considered, the parallel state is 
always the most stable regardless of the shape of 
the specimen. This does not imply electric satu- 
ration for the same reasons as given above for the 
magnetic case. 

I wish to express my gratitude to Professor 
Fowler for suggesting this problem and for his 
generous assistance. I am also indebted to Mr. 
A. H. Wilson and Professor Van Vleck for helpful 


-discussions. 
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Within the framework of the general theory of relativity, it is proposed to introduce at each 
point of space-time a Euclidean metric tensor yy, in addition to the usual Riemannian metric 
tensor g,». In this way one imparts tensor character to quantities which in the usual form of the 
theory do not have it. For example, one can obtain a gravitational energy-momentum density 
tensor in place of the usual pseudo-tensor. Furthermore one can impose four additional co- 
variant conditions on the gravitational field and thus restrict the form of the solution for the 
field corresponding to a given physical situation. 


§1 
N the general theory of relativity, Einstein’ 
makes use of Riemannian geometry, charac- 
terized by the existence of an invariant interval 
or line element ds, given by 


ds? = g,,dx*dx’, (1) 


where g,, is the metric tensor and (x!,x?,x3,x*) 
are the coordinates in an arbitrary coordinate 
system. The metric tensor is determined by the 
law of gravitation: 


Gy = —8nT,,, (2) 


where 7,, is the energy-momentum density 
tensor, and 
Gy» = Ry (3) 
with R,, the contracted Riemann-Christoffel 
tensor. 
The tensor G,, satisfies the identity 


G,%. =0, (4) 


where (;) denotes covariant differentiation based 
on gy (g-differentiation). This permits setting 
G,» proportional to 7,, which is taken to satisfy 
the equation (expressing the laws of motion of 
matter or energy) 


=0. (5) 
§2 


It is now proposed to introduce at every 
point of space-time, along with the metric tensor 
£., a second metric tensor y,, corresponding to 
flat space, i.e., for which the Riemann-Christoffel 
tensor vanishes identically everywhere. This may 
be interpreted in various ways. For instance, we 


1A. Einstein, Ann. d. Physik 49, 769 (1916). 


may suppose that we map the Riemannian space 
with metric g,, on a flat space with metric y,, and 
assign the same values to the coordinates of 
corresponding points. Or we may say that the 
two metrics side by side represent a comparison 
of the given space with the space one would 
have if the gravitational field were removed. It 
should be emphasized that in introducing the 
Yu», We are not postulating here any new proper- 
ties of the space. 

In this way one can define a Euclidean line 
element analogous to (1) 


= y,x*dx’. (6) 


Moreover one can now define covariant differ- 
entiation based on y,, (y-differentiation), which 
will be denoted by (,). Since the Riemann- 
Christoffel tensor formed from y,, has been 
assumed to vanish, it follows that one can 
interchange the order of y-differentiation, so 
that the latter obeys all the rules of ordinary 
differentiation, except that the y-derivatives of 
Yu» vanish. Indeed it is always possible to choose 
a coordinate system (special relativity) such 
that the components y, are constant and 
y-differentiation reduces to ordinary ; but it may 
not always be convenient to choose such a system. 
Let us next consider A*,, defined by 


where {} is the Christoffel 3-index symbol 
used in g-differentiation and I,, that occurring 
in y-differentiation. This can be shown to be a 
tensor and is found to be given by 


Me = + Bran — Sura) (8) 
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If in the usual expression? defining R,,, one 
substitutes (7), it is found that, as a consequence 
of the Euclidean character of y,,, one can write 


This relation is interesting in that R,, is seen 
to be obtainable from the tensor A*,, by tensor 
operations, i.e., R,, is a tensor function of gy». 


§3 


Comparing (9) with the usual expression for 
Ry», one sees that {} has been replaced by 
A*,, and ordinary differentiation by y-differentia- 
tion. This type of correspondence turns out to be 
quite general. It appears that one can rewrite all 
the quantities occurring in relativity theory so 
that {)} is replaced by A*,,, an ordinary deriva- 
tive by a y-derivative (in particular, dg,,/dx" by 
Sur, «) and (—g)! by 


x=(g/y)', (10) 


where g and y are the determinants of g,, and Y,», 
respectively. Since g,»,. is a tensor and « a scalar, 
it follows that some quantities, on being re- 
written in this way, take on a tensor character, 
although they do not have it in the usual form 
of the theory. 

To illustrate this point, the following example 
may be useful to compare with the analogous 
calculation usually made.* One can show that 
the gravitational Eqs. (2) are connected with a 
variational principle 


Le(—7)'dr) 


f (11) 


where 
dr=dx'dx*dx*dx', 


L=g* ar A*apA*,») (12) 


and only the components of g,, are varied (with 
the variations vanishing on the boundaries of 
the region of integration). This function L is a 
scalar. From it one can derive a gravitational 


’ 2E.g., L. P. Eisenhart, Riemannian Geometry (Princeton, 
1926), 21. 

* Reference 1, p. 804, and H. Weyl, Raum Zeit Materie 
(Berlin, 1923), p. 272. 
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energy-momentum density tensor (analogous to 
the pseudo-tensor of Einstein). One writes (Eq. 


(4)) 


‘ G,* +G*,A%sa =0, ( 13) 

or, since 
(14) 
(xG%,) — gag, =0. (15) 


Throughout this paper, unless otherwise speci- 
fied, indices are raised and lowered by means 


of gy». 
Now (11) can be shown to be equivalent to 


the relation 
= 0(Lx) (16) 


and from this one obtains (using integration by 
parts) 


{Ln — [0(Lk) (17) 
which is to be substituted in (15). If we now 
define as the gravitational energy-momentum 
density tensor, 


1 
t,*=——{ — [0(Lx) (18) 
16r 


we can write (15) or (4) as 


[x(T+1.) (19) 

Or, taking 
(20) 

we can write 
x= 0. (21) 


By a suitable coordinate transformation, one can 
put the left-hand member into the form of an 
ordinary divergence. The fact that ¢’,’ is a tensor 
serves to remove some objections that had been 


raised in the past.‘ 
§4 


It is to be seen that, in general, the quantities 
in the usual form of the relativity theory can be 
considered as arising from the corresponding 
ones here, through the special choice of y,», 


Yu» =constant. (22) 


*H. Bauer, Physik. Zeits. 19, 163 (1918). 
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However, if (22) is to hold for all coordinate 
systems, Y,» cannot be a tensor. Consequently 
the quantities which depend on y,, lose their 
tensor character. Thus the gravitational tensor 
t, of (18) goes over into the pseudo-tensor in 
the Einstein theory. On the other hand, R,, is 
independent of y,, and hence remains a tensor. 
Obviously, if one wishes to maintain the tensor 
character of quantities depending on +,, one must 
give up (22) and allow y,, to transform as a 
tensor under coordinate transformations. 


§5 


Although g,, and y,, have been considered as 
existing side by side, so far nothing has been 
said about any relation between them. Actually 
one would expect to have some relation between 
them, for it seems reasonable that, if the gravita- 
tional field is made to vanish, g,, should go over 
into Yu. Such_a relation can be obtained by 
imposing four additional covariant conditions on 
the field. This can be done because of the set 
of identities (4) existing among the field equa- 
tions. In working with the linear approximation 
of the gravitational equations, Einstein‘ similarly 
added four conditions on the field for the purpose 
of eliminating (or reducing) apparent fields 
arising from infinitesimal coordinate transforma- 
tions. It has hitherto not been possible to set 
up corresponding (covariant) conditions for the 
exact equations. In the present form of the 
theory this can be done, and it is these equations 
that provide the relation between. g,, and y,». 
They should serve to remove, or at least restrict, 
the ambiguity in the form of a solution arising 
from the possibility of coordinate transforma- 
tions which go over into the identical trans- 
formation in the absence of a field. 


For example, the static spherically symmetric solution 
of the gravitational equations in free space can be written 
in a number of forms depending on the choice of the radial 
variable r, it being asserted that one cannot know which 
choice corresponds to the variable r in flat space. The 
additional conditions would serve to single out one of the 
various possibilities. 


One reasonable set of conditions can be ob- 
tained from the following considerations: 
One can rewrite the expression (9) for R,, in 


‘A. Einstein, Berl. Ber., p. 688 (1916). 


the form 
= 387 — piv — Sip — 
—A% pA” age — A% as*, (23) 
Sy = Maat = fay Ky (24) 


where 


and an asterisk with an index indicates raising 
or lowering it with g,,. 
This suggests taking as additional conditions 


S,=0. (25) 


In consequence, the right-hand member of (23) 
is simplified, only the first term now containing 
second derivatives. It can be verified that the 
conditions (25) in first approximation are the 
same as those imposed by Einstein in the case of 
the linear equations. 

The static spherically symmetric solution of 
the field equations in free space 


(26) 
together with (25), is found to be 


r+m 
ds? = ————dr® — 


rT—m 


r—m . 
+——di?, (27) 
r+m 


with the usual notation, if one takes 
do? = (28) 


Another possible set of conditions can be obtained based 
upon taking 
«=i, ° (29) 
This serves to simplify many of the formulas. It should be 
noted that this differs very essentially from the condition 
g=-—1, often used, in that it is not affected by coordinate 
transformations. Three further conditions are required. One 
might take (as being equivalent to three conditions) 


Srp = 0. (30) 


In this case if one takes for de® the expression (28), 
the static spherically symmetric solution has the Schwarz- 
schild® form: 


ds? = —(1—2m/r)“dr* 
—1*(d#+sin* (31) 


However this solution leads to the difficulty pointed out by 
Schrédinger’ (a not very serious difficulty, however*) that 


® K. Schwarzschild, Berl. Ber., p. 189 (1916). 

7 E. Schrédinger, Physik. Zeits. 19, 4 (1918). 

8 A. Einstein, Physik. Zeits. 19, 115 (1918); Berl. Ber., 
p. 448 (1918). 


to 
13) 
14) 
15) 
ci- 
ins 

to 
16) 
by 
17) 


150 N. ROSEN 


the components of t,” all vanish for r>2m. If one wishes to 
avoid this,one must give up (29) and (30). 


§6 


Enough has been given to show that there are 
advantages from the formal point of view in 
introducing the Euclidean y,, into the general 
relativity theory. It imparts tensor character to 
quantities which otherwise do not have it, and 
allows additional conditions to be imposed on 


the field so as to restrict the form of the solution 
for a given physical situation. 

In conclusion, it is necessary to point out that, 
having once introduced y,, into the theory, one 
has a great number of new tensors and scalars at 
one’s disposal. One can, therefore, set up other 
field equations than (2). It is possible that some 
of these may be more satisfactory for the descrip- 
tion of nature. Further investigation is here 
required. 
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The possibility is considered of interpreting the formalism of the general theory of relativity 
in terms of flat space, the fundamental tensor g,, being regarded as describing the gravitational 
field but having no direct connection with geometry. The resulting theory in general leads 
to the same predictions as the Einstein theory, but there are cases where the predictions 
differ. The present theory may explain the principal results obtained by D. C. Miller in his 
“ether-drift’’ experiments. The implications of the theory for cosmology are briefly touched 


upon. 


§1 

N a previous paper' (hereafter referred to as I) 

it was shown that it is useful to introduce 
into the general theory of relativity the concept 
of the existence at each point of space-time of a 
Euclidean metric tensor y,, in addition to the 
usual Riemannian metric tensor g,,. From the 
standpoint of the general theory of relativity, one 
must look upon y,, as a fiction introduced for 
mathematical convenience. However, the ques- 
tion arises whether it may not be possible to 
adopt a different point of view, one in which 
the metric tensor y,, is given a real physical 
significance as describing the geometrical proper- 
ties of space, which is therefore taken to be flat, 
whereas the tensor g, is to be regarded as 
describing the gravitational field.” 

It has been pointed out in I that the introduc- 
tion of y,, leads to the possibility of other laws 


1N. Rosen, Phys. Rev. 57, 146 (1940). 

2 In some respects this resembles the theory of gravita- 
tion proposed by Nordstrém (cf. report by M. v. Laue, 
Jahrbuch f. Rad. u. El. 14, 263 (1917)). It will be seen that 
there are important differences, however. 


than those adopted in general relativity. In the 
present paper, however, no attempt will be made 
to change the laws of the latter, since they form 
a self-consistent system and have proved to be 
quite satisfactory for the description: of large 
scale phenomena, at any rate. 


§2 


As far as the field equations are concerned, it 
is immaterial what interpretation one gives to 
the variables involved. This is not the case with 
the equations of motion. Let us therefore consider 
the law of motion for a particle in the field. 
The latter is given in the general theory of 
relativity by the equation of the geodesic’ 


(1) 
ds* \|aBj)ds ds 


where ds is the line element, defined in terms of 
the tensor g,,, (I (1)). Let us now introduce as 
independent variable the Euclidean line element 


3 A. Einstein, Ann. d. Physik 49, 769 (1916). — 
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do, defined in terms of the tensor y,», (I (6)), 
instead of ds. After canceling out a factor do/ds, 


the equation can be written 


d sdo dx* dx* dx*® do 
do\ds doe a B) do do ds 


This can be rewritten 
dx* dx® dx* (3) 
do da da da 
d*s sds 
where r=— / —= (4) 
and Tag and A*gg are given by I, (7) and (8). 


If we multiply (3) by y,,.dx’/de, the left-hand 
member vanishes as a consequence of the defini- 
tion of do. Hence it is necessary that 


dx* dx® dx’ 5) 

do do da 

which gives the relation between s and co. If this 
expression for \ is substituted into (3), one gets 
an equation from which s has been completely 
eliminated. 

We can get the equation in a more familiar 
form if we go back to (2), multiply it through by 
a constant mo, and write 


where m will be considered as the “proper mass”’ 
of the particle. Then the equation can be 
written as 


d dx* dx* dx® dx* dx® 
—— —= —mA*.s—- —-. (7) 
da da do 


On the left we have the absolute derivative 
with respect to o of 
dx* dx* (8) 
}, 
do ds 
which is the energy-momentum vector. Hence the 
right-hand member is to be interpreted as the 
gravitational force. We see from (7) that, 
although inertial and gravitational mass are 
equal, they are not necessarily constant. 
If we define the scalar 


dx* dx’ 


(9) 


do 


then (6) can be written 
m=mo/w (10) 


as the law of dependence of proper mass on the 
gravitational field. 

Conversely, if we start from (7), multiply by 
£udx"/do and make use of I (7), we get 


3mdu*/do = —w*dm/de, (11) 


which is equivalent to (10). 

Thus we see that, from the standpoint of flat 
space, one can put the equations of motion of a 
particle into the Newtonian form (generalized to 
four dimensions) if we assign to the particle a 
variable mass. In a coordinate system in which 
the particle is at rest the dependence of mass on 
the field is given by 


m= (12) 


In a weak gravitational field this can be written 
in first approximation‘ 


m=m)(1—¢), (13) 


where ¢ is the gravitational potential. 

In the case of a light ray we have the addi- 
tional condition (in the presence of a gravita- 
tional field) 

w=0. (14) 


We therefore let mo vanish also, but take m, as 
given by (10), to be finite and constant. This 
allows one to cancel out the m’s in (7). That (14) 
is consistent with (7) is to be seen from (11). 


§3 ‘ 

To see the connection between Eq. (7) and 
the gravitational field equations, (I (3)), let us 
suppose we have a stream of particles each of 
which moves according to (7). Let mo be the 
proper particle density and let us assume that all 
quantities that we may need are smooth functions 
of the coordinates. We write (7) in the form 


dx* dx* 


o do 


and multiply through by mo. On the left we have 


dx* dx* dx* 
( mne— No =). 
da do 


‘Since it is shown, iia 3, that for a weak field 
£44 =1+2¢ in Galilean coordinates. 
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The scalar (mo(dx*/do)), . can be evaluated in a 
Galilean coordinate system. Thus 


dx* dt dx* dx* 
= —) = (— (16) 
do dt/ « dt] 


where ¢ is the time coordinate and 
(17) 


is the particle density in a reference frame with 
respect to which the particles at the point in 
question have a velocity v, the velocity of light 
being taken as unity. Hence 


8) 


dx* On A(nvz) A(nvy) A(nv,) 
( + + 
da/ dat Ox oy dz 


and this vanishes (equation of continuity) if the 
number of particles is conserved, which we shall 
assume. If we set as the proper mass density 


(19) 
(14) becomes 
( dx* dx* dx* dx® (20) 
da/ a do do 


Taking as the energy-momentum density 
tensor 
dx* dx” (21) 
w= ’ 
do do 
where the prime has the same significance as in 
I (20), one readily verifies that I (5) and hence 
I (21) are satisfied. In this way we see that the 
equations of motion and the field equations are 
consistent with one another. 


§4 


The preceding calculations show that the 
standpoint of flat space at any rate does not 
lead to internal contradictions. Let us now com- 
pare the theory based on this standpoint with 
the Einstein general relativity theory. 

First of all, insofar as the three crucial tests® 
are concerned the two theories lead to the same 
predictions. Since the equations of motion have 
been taken the same, except for the independent 
variables, it is clear that relations among the 


5 A. S. Eddington, The Mathematical Theory of Relativity 


(Cambridge, 1924). 
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coordinates, giving the motion of a particle or 
the path of a light ray, will also be the same in 
both theories. But more than that, the gravita- 
tional red shift will be the same in the two 
theories, although this is not quite as obvious. 

The gravitational red shift in the general 
relativity theory is usually derived from the 
assumption that the interval ds connected with 
the period of oscillation of an atom dt remains 
unchanged if the atom is put into a gravitational 
field, together with the relation 


ds? = 


to which I (1) reduces if the space coordinates 
of the atom are fixed. It is clear, however, that 
this change in the period of oscillation must also 
be a consequence of the equations of motion of 
the particles constituting the atom. These equa- 
tions being the same in the two theories, the 
predicted red shift will also be the same. (It 
should be mentioned that the equations of 
motion referred to will involve the electromag- 
netic field and may be of a quantum-theoretical 
nature.) 

One way of describing the situation is to say 
that all the formal relations involving ds are 
taken over into the present theory but ds is not 
regarded as the geometrical interval. Putting the 
matter in this way makes the comparison of the 
two theories almost trivial. Thus one sees, for 
example, that the principle of equivalence is 
valid here as long as we deal with coordinates 
and not with intervals. 

Indeed it appears at first sight as if the differ- 
ences between the theories were only questions 
of conventions. 


For example, if a small measuring rod is moved from 
one place to another in a gravitational field, the interval 
ds between its end points will remain constant, but the 
corresponding de may change.* One may say that the 
length of the rod does or does not change, depending on the 
point of view. As far as the changes in the coordinate- 
differences of the end points are concerned, however, no 
disagreement exists. Incidentally, there is some analogy 
here to the change in length of a measuring rod with 
temperature. If all substances had the same temperature 
coefficients of expansion, one could assume that the length 
of the measuring rod is independent of the temperature and 
no contradictions would arise. This corresponds to the 
point of view of the Einstein theory. 


® Reference 3, p. 818. 
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A similar example involving time is to be found in the 
case of the gravitational red shift. 


Actually it turns out that there are other, 
more important, differences between the two 
theories. They are based upon the fact that, 
after all, the fundamental ideas of the two 
theories are different. The present theory involves 
less relativity. This may be seen from the follow- 
ing example: 

Let us take for the moment the flat-space 
point of view. Suppose that in a given frame of 
reference, we have a static gravitational field 
described by g,,». Then the velocity of light will, 
in general, be different from that in field-free 
space. The space containing the field can there- 
fore be looked upon as a medium with an index 
of refraction different from unity. From the 
standpoint of special relativity one can expect 
that such a medium will exert a “drag” upon 
light, as in the case of the well-known experiment 
of Fizeau. If we go over to a second frame of 
reference moving uniformly with respect to the 
first (Lorentz transformation), it will be possible, 
by measurements carried out on light in this 
system, to determine its motion relative to the 
other system. Thus the static gravitational field 
determines an “absolute” frame of reference. 
This will be true whenever the velocity of light 


differs from its special-relativity value, even in. 


the limiting case when the components g,, are 
constant in a part of space, so that there is no 
gravitational force (but the gravitational poten- 
tial differs from zero). 

Such a situation may account for the principal 
results of the ‘‘ether-drift’’ experiments of D. C. 
Miller.’ A discussion of this question will be 
presented in a separate paper, since it involves 
some considerations which would lead us too far 
afield here. 

If we consider the same situation from the 
point of view of the Einstein theory, we arrive 
at entirely different conclusions. We cannot 
carry out a Lorentz transformation in the 
presence of a gravitational field but must first 
get rid of the latter by going over to a suitably 
accelerated frame of reference. Having obtained 
a coordinate system in which the components g,, 


7D. C. Miller, Rev. Mod. Phys. 5, 203 (1933). 
SI am indebted to Professor Einstein for a helpful 
discussion of this question. 


are constant (at least for a small portion of 
space), we can always make a further coordinate 
transformation so as to give them their special 
relativity values. But then no light-drag or 
“ether-drift”’ can be observed. 


§5 

The question of viewpoint being discussed 
here has a bearing on the application of rela- 
tivity theory to cosmology. Thus, if one takes 
seriously the interpretation in terms of flat 
space, one must give up the A-term in the 
cosmological field equations,® since this prevents 
the g,, from going over into the y,, at large 
distances from matter. Of course, any solution 
of the equations corresponding to an “expanding 
universe’’ must be interpreted as representing 
the expansion of a distribution of matter and 
radiation in a flat space. There appears to be no 
necessity for imposing the condition of homo- 
geneity, as is usually done. 


§6 

In the Einstein general relativity theory 
gravitation is explained in terms of geometry. 
In the theory suggested here based on flat space, 
this geometrization of gravitation has been given 
up. Perhaps this may be regarded by some as a 
step backward. It should be noted, however, that 
the geometrization referred to has never been ex- 
tended satisfactorily to other branches of physics, 
so that gravitation is treated differently from 
other phenomena. It is therefore not unreason- 
able to wonder whether it may not be better to 
give up the geometrical approach to gravitation 
for the sake of obtaining a more uniform treat- 
ment for all the various fields of force that are 
to be found in nature. 

It might be remarked that the theory, as 
presented here, was obtained by a reinterpreta- 
tion of the Einstein general relativity theory. 
This amounts to using the formalism of the 
latter as a means for obtaining a consistent set 
of equations. Such a procedure is not entirely 
satisfactory. It ought to be possible to build up 
the theory independently from suitable postu- 
lates. This will be attempted in the near future. 


* A. Einstein, Berl. Ber. p. 142 (1917). As a matter of 
fact the existence of the nebular red-shift weakens the 
chief argument advanced there for introducing the \-term. 
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Note on Ether-Drift Experiments 


N. RosEN 
Massachusetts Institute of Technology, Cambridge, Massachusetts 
(Received October 30, 1939) 


It is shown that it is possible to account for the principal results obtained by Miller in his 
“ether-drift” experiments, on the basis of a flat-space theory of gravitation, as due to the 
motion of the earth in the gravitational field of all the matter and energy in the universe. 


OME time ago, D. C. Miller presented a 
report! summarizing the large number of 
“ether-drift” experiments which he had carried 
out to determine the absolute motion of the earth 
through space. He concluded from his experi- 
ments that the earth is moving through space 
with a velocity of several hundred km/sec., but 
that the effect of this motion on the velocity of 
light (as measured on the earth) corresponds to a 
reduced velocity of about 10 km/sec., so that 
there is a drag-coefficient which he finds to have 
a value of about 0.051. 

Unfortunately these results have not received 
satisfactory confirmation from other observers. 
Furthermore, Miller finds that the projection of 
the earth’s translatory motion on the surface of 
the earth in the course of the day takes on 
directions which are unsymmetrical with respect 


to the north and south line. This seems difficult to. 


explain on the basis of any simple theory. It may 
conceivably be an indication of some extraneous 
effect which accounts for his other results as 
well. However, the experiments appear to have 
been carried out very painstakingly and the 
results he obtained on the whole are remarkably 
consistent. Undoubtedly his experiments deserve 
_ serious consideration. Yet thus far no satisfactory 
explanation has been proposed. 

In a recent paper by the writer,? the idea was 
considered of interpreting the formalism of the 
general theory of relativity in terms of flat space 
(special relativity). There, mention was made of 
the possibility of explaining the main result 
obtained by Miller on the basis of this idea. In 
the present note this question will be treated in 
some detail. 

Let us suppose that, in a frame of reference S, 
we have a weak static gravitational field de- 


1D. C. Miller, Rev. Mod. Phys. 5, 203 _— 
2 N. Rosen, Phys. Rev. 57, 149 (1940). 


scribed by the gravitational potential (¢< 0). 
Let us take a Galilean coordinate system in S 
with coordinates (x,y,z,t) so that the Euclidean 
line element do is given by 


do? = (1) 


the velocity of light in field-free space being 
taken as unity. Then the Riemannian line ele- 
ment of general relativity is given in first 
approximation by? 


ds*= (2) 


The velocity of light u is determined by the 


condition 
ds=0 (3) 


and one finds that it is independent of direction 
and is given in first approximation by 


u=1+2¢. (4) 


We see that the space can be looked upon as a 
medium having an index of refraction n, where 


n=(1+29)=1-—2¢. (5) 


From the law of addition of velocities in the 


special theory of relativity it follows‘ that a 
medium of index of refraction m has a drag 


coefficient 
k=1—1/n’. (6) 


This means in our case that, if there is an 
observer in a frame of reference S’ moving with a 
uniform velocity v with respect to S, he will find 
that the velocity of light in the direction of 
motion will be 


uy’ =u—ke (7) 
and that in the opposite direction 

ue’ =u+kv, (8) 
Vol II, p. 280. 
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where, in first approximation, 
k= —Ag. (9) 


One can verify (7), (8), and (9), neglecting 
terms of the order of v?, by carrying out a Lorentz 
transformation from coordinates (x,y,z,t) in S to 
(x’,y’,2’,t’) in S’. This leaves (1) unchanged in 
form when expressed in terms of the primed 
coordinates. From (2) and (3) then follow the 
preceding equations. 

Thus we see that a static gravitational field 
provides a frame of reference with respect to 
which the uniform motion of an observer can be 
determined by experiments on light performed 
in his own system. It should be emphasized that 
this does not require the presence of a gravi- 
tational force; all that is necessary is that the 
potential ¢ be different from zero. | 

From the value 0.051 for k, as obtained by 
Miller, one has 


o=—0.013 (10) 


in the dimensionless units of general relativity 
theory. The question naturally arises: where does 
a potential of so large a magnitude come from? 
The obvious answer seems to be that this 
potential is due to the distribution of matter 
throughout the universe. Since we do not have 
enough information about the latter to carry out 
an accurate calculation, we shall restrict ourselves 
to a rough estimate of the order of magnitude. 
The gravitational potential at the center of a 
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sphere of density p and radius R is given by 
g= —2rxpR?/c, (11) 


where « is the gravitational constant and c the 
velocity of light in field-free space. Let us take for 
p the order of magnitude of the estimated mean 
density of matter in space: 


(12) 


From (11) one can calculate the value of R which 
will give the value (10) for ¢. It is found to be 


10° light-years. (13) 


To obtain the same value of ¢ for a point not at 
the center of the sphere, one must take a larger 
value for R, of course. The value given in (13) is 
not unreasonable as an “effective radius’’ of the 
universe. 

Thus we see that it is possible to account for 
Miller’s observation of an ‘‘ether-drift,”’ together 
with the value he finds for the drag-coefficient, 
on the basis of a flat-space theory of gravitation, 
as being due to the motion of the earth through 
the gravitational field of all the matter (and 
energy) in the universe. This field, the resultant 
of a large number of small contributions, is nearly 
static and serves to single out a frame of reference 
for determining “absolute” motion. The asym- 
metry referred to above remains unexplained. 

It seems highly desirable to have more experi- 
mental work done on this question. 


5 E. Hubblé, Astrophys. J. 79, 8 (1934). 


I 


in 
1) 
| 
e- 
st 
2) | 
he 
3) 
on 
| 
4) 
a 
re 
5) | 
he 
a 
7 | 
6) 
an 
la 
nd 
of 
8) | 


JANUARY 15, 1940 


PHYSICAL REVIEW 


VOLUME 57 


LETTERS TO THE EDITOR 


Prompt publication of brief reports of important discoveries in physics may be secured by 
addressing them to this department. Closing dates for this department are, for the first issue of the 
month, the eighteenth of the preceding month, for the second issue, the third of the month. Because of 
the late closing dates for the section no proof can be shown to authors. The Board of Editors does 
not hold itself responsible for the opinions expressed by the correspondents. 


Communications should not in general exceed 600 words in length. 


The Propagation of Electron Waves in Ionic Single Crystals 


The interaction between conducting electrons and the 
lattice structure of insulating crystals can be studied by 
the effects preceding and accompanying electric break- 
down.' Especially striking is the dependence of the 
breakdown direction on the crystallographic orientation of 
the sample. A plate of rocksalt, for instance, cut parallel 
to the cleavage planes (100) and exposed to an intense 
homogeneous field perpendicular to it, does not break down 
first in the direction [100] of highest stress but shows a 
breakdown path parallel to the face diagonal [110] or, 
if overvoltage is applied, in the [111] direction. This 
“direction effect’’ can be explained by the potential 
structure of the crystal (Fig. 1): A surplus electron, in 
traveling through the lattice, encounters the lowest 
potential humps in the [110] direction and the smallest 
number of lattice points in the [111] direction, while the 
cube edge [100] is unfavorable in both respects. It might 
be expected that this directional effect is paralleled by a 
_ “magnitude effect.”’ If the breakdown field strength for a 
plate cut parallel to (100) is Emax, it should be Emax/2! for 
a plate orientated parallel to the dodecahedral plane (110) 
because the face diagonal now points into the field direc- 
tion. Accordingly a plate cut parallel to the octahedral 
plane (111) should break down at Emax-34/2 if the break- 
down field strength in [110] direction is decisive, that is, 


Emax( 111] _ 
Emax[110] 


Preliminary experiments carried out some time ago? did 
not fully confirm this expectation. While (111) plates 


10 
Fic. 1. Potential structure of 8 
the crystal. 
+ + + 
A 
+ 4 - + - ™ 


showed the anticipated lower value of dielectric strength, 
the (110) plates withstood about the same stress as crystal 
plates cleaved in [100] direction. Both results were of 
dubious value because of the difficulty of faultlessly 
preparing very thin and polished plates of NaCl in direc- 
tions differing from the natural cleavage plane. 

In the meantime the question as to the existence of the 
magnitude effect has become of new importance. Following 
a note of Zener,? a wave-mechanical theory of the break- 
down process has been carried through by Franz‘ based 
on the assumption that at very high field strengths the 
distinction disappears between conductor and insulator. 
According to the simple zone theory, the insulator is 
characterized by a completely filled lower band and an 
empty higher zone separated from the first one by about 
three volts. The forbidden interspace is created by Bragg 
reflections of the electron waves on the lattice planes, 
resulting in standing waves and thus preventing the 
acceleration of electrons. This mechanism becomes in- 
efficient under very high field strengths because the energy 


difference between the filled band and the conducting 


levels can be overcome by accelerating the electrons 
through relatively few lattice planes, which are not 
sufficient for setting up strong reflected waves. In con- 
sequence, electrons leak through; the insulator becomes a 
conductor. If the number of planes traversed in the 
different lattice directions per unit length is the decisive 
factor for the intensity of reflection preventing conduction, 
a magnitude effect should be expected: 


Emax(111] 
Emaxf110] 


as above, but furthermore 


Emax(111] _ 
Emax[100] 


in contradiction to the direction effect. 

Since every theory of electric breakdown has to know 
how far directional effects enter into the value of the field 
strength measured, the authors undertook a careful study 
of the existence of the magnitude effect in NaCl. The 
crystal plates were prepared by a new and simple method 
allowing one to secure faultless samples for every crystallo- 
graphic direction in shortest time. After cutting thick 
samples in the right orientation out of synthetic single 
crystals, the material was ground down on glass plates 
under water to about 0.1 mm thickness. Wiped off with 
acetone, the crystal got a high polish, and heating in an 
electric furnace finished the treatment. Crystal plates cut 
parallel to the octahedral plane (111) showed some 
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tendency to crack, resulting in lower breakdown values. 
But after eliminating every disturbing effect, the measure- 
ments demonstrated without a doubt that there does not 
exist a magnitude effect in rocksalt, at least not at room 
temperature. In every direction tested, the breakdown 
strength is 1.65 million volts/em with an error of about 
+5 percent. This error is not so much due to the measuring 
technique, which would permit about two percent, as to 
inhomogeneousness of the material and roughness of the 
electrodes, resulting in locally higher stresses. Very slow 
raising of the voltage seems to have an equalizing effect, 
as already pointed out elsewhere;® in consequence the 
average breakdown field strength here reported lies about 
10 percent higher than our earlier value. 

The nonexistence of the magnitude effect disproves the 
wave-mechanical theory as given above. It is at all events 
hard to believe that the electrons of the Cl- and Na* are 
not locally bound but move freely inside of their filled 
band through the whole crystal. But if no magnitude 
effect exists, the ionization theory of breakdown proposed 
by one of us and calculated in different ways by H. 
Frohlich,* R. J. Seeger and E. Teller,? and W. Franz‘ also 
seems to come into difficulties. Our result shows that 
apparently the electrons do not feel the difference of the 
lattice directions until they have been accelerated to some 
extent. But this has to be expected, as Dr. Teller kindly 
advises us, because of the wave-length of the electrons. 
The breakdown strength of 1.65 million volts/cm corre- 
sponds to about 4.7 X 10? volt potential difference between 
the ions in [100] direction. The electrons start therefore 
with energies of the order of 10~' volt corresponding to a 
wave-length of about 410-7 cm, that is, one electron 
covers about 15 ions in the [100] direction. This spread 
of the electron waves over a large lattice area makes them 
disregard the structure until they have shrunk by accelera- 
tion to the appropriate size. 

A. von HIPPEL 


J. W. Davisson 
Electrical Engineering Department, 
Massachusetts Institute of Tec ‘ 
Cambridge, Massachusetts, 
December 18, 1939. 


1A. von Hippel, Z App. Phys. 8, 815 (1937). 

2A. von Hippel, Zeits. f. Physik 75, 161 (1932). 

3C. Zener, Proc. Roy. Soc. London 145, 523 (1934). 

4W. Franz, Zeits. f. Physik 113, 607 (1939). 

5 R. C. Buehl and A. von Hippel, Phys. Rev. 56, = (1939). 
Frohlich, Proc. Roy. Soc. London 160, 230 (19. 37). 

7A. J. Seeger and E. Teller, Phys. Rev. 54, 515 (1938). 


The Nonexistence of Transuranic Elements 


It is not clear why no atoms of atomic number greater 
than 92 are found in nature. The properties of the heaviest 
known nuclei do not suggest that all such atoms would be 
extremely unstable so that some reasonably long-lived 
alpha-emitting transuranic atoms might well be expected. 
Further, the absence of U™* is puzzling, for one would 
expect this nucleus to be as stable as U** and U™*. It is 
striking that the periodic table of the elements ends with 
the only one that gives nuclear fission with thermal 
neutrons. This is possibly more than a coincidence and 
may be the clue to the explanation of the missing atoms. 
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According to the theory of Bohr and Wheeler! the 
probability of fission of any nucleus increases rapidly with 
increase of the energy of the captured neutron above a 
critical lower limit. The nuclei that give fission with slow 
neutrons are thus the ones to give the largest yields for 
bombardment by faster neutrons. Bohr and Wheeler calcu- 
late that of the reasonably stable nuclei (half-lives > 10* yr.) 
only two should give fission with thermal neutrons. These 
are U™ and U**, Pa*! should be so disintegrated by 
neutrons of energy greater than about 0.1 Mev. Atoms of 
lower atomic number have much higher critical energies so 
that the stable nuclei from Bi on down would require 
neutrons of extremely great energy for fission to be possible. 
Proceeding in the opposite direction to higher atomic 
numbers one finds that all reasonably expectable trans- 
uranic atoms should give fission with thermal neutrons. 
Further, U*? should also do this and U™* should require 
neutrons of only 0.1-0.2 Mev of energy. Even if these U 
nuclei and the transuranic nuclei had once been abundant 
they could have been destroyed by neutrons in a later 
phase of the development of the solar system. The presence 
of such neutrons is entirely hypothetical. Von Weiszaecker* 
has concluded that it is unlikely that there is any con- 
siderable concentration of neutrons in stellar atmospheres 
but that the neutron density must have once been high if 
the relative proportions of the different atomic species 
are those of an equilibrium at a high temperature. 

The existence of the long-lived U** (AcU) atoms seems 
at first glance to be in disagreement with this hypothesis. 
The postulated neutron bombardment, however, would 
have produced some U**® by capture. This is known to 
decay with a half-life of 23 min. by the emission of beta- 
rays. The resulting »;EkaRe*® is most probably an ancestor 
of U*5, although the details of the intermediate transitions 
are not yet known. If the half-lives of the intermediate 
nuclei are short enough so that no important number of 
them would have been lost by fission some amount of U™* 
would have been maintained in a balance between this 
production from U** and loss by fission. 

Neutrons would also have been captured by Th®™. The 
Th** so produced has been found to emit beta-rays and 
decay with a half-life of 25 min. The half-life of the re- 
sulting Pa** is not known. The absence in nature of all 
of the atoms of the 4n+1 family to which these belong is 
understandable if the presumably long-lived U**’ was lost 
by fission without the replenishment as with U™*, and if 
neither the Pa** nor any one of its descendants are long- 
lived. If the Pa** decays by beta-emission to U™ this 
would have disappeared either because of being naturally 
short-lived or because of being especially susceptible to 
loss by fission. 

, These considerations will be developed in greater detail 
in a paper to be submitted to The Physical Review for 
publication. 

Louis A. TURNER 

Palmer Physical Laboratory, 

Princeton University, 


Princeton, New Jersey, 
December 31, 1939. 


1N. Bohr and J. A. Wheeler, Phys. Rev. 56, 426 and 1065 (1939). 
?C. F. von Weiszaecker, Physik. Zeits. 39, 633 (1938). 
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A Measurement of Mesotron Lifetime 


To obtain a proof of the mesotron decay hypothesis and 
to determine their lifetime without the necessity of making 
assumptions regarding their energy spectrum and the 
height of the layer where they are produced, the following 
observations have been made. The absorption curves of 
the cosmic radiation of vertical incidence were measured 
in lead from 0 to 15 inches Pb first on top of Mt. Mitchell, 
North Carolina (elevation 2040 meters, mean barometric 
pressure 60 cm Hg) and then at Durham, North Carolina 
(elevation 123.5 meters, mean barometric pressure 75 
cm Hg). In the first series (Mt. Mitchell) the air column of 
15 cm Hg=203 g/cm* was compensated by a layer of 
graphite of 200 g/cm? filling the whole sensitive cone of 
the counter telescope. 

In these measurements we have counted coincidences 
between four G-M counters. The slabs of absorbing lead 
were three times the width of the counters and about 
three times the sensitive length of the counters. 

Some of the observed counting rates per hr. were thus: 


Inches Pb 0 3 8 15 
Mt. Mitchell 46.8 38.2 30.1 25.2 
Durham 32.7 28.0 24.3 22.2 


We estimate the statistical probable errors of the Mt. 
Mitchell data to be of the order of 1.2 percent and that of 
the Durham data to be of the order of 1.5 percent. 

A similar series of observations at Mt. Mitchell with 
10 cm of lead at the sides of one counter shows that the 
background showers were quite negligible in their effect 
on the course of the absorption curve. Furthermore, the 
number of such sidewise showers was small compared to 
the difference in counting rates in the two series of data. 
This large difference between the two series of measure- 
ments at Mt. Mitchell and Durham is therefore real and 
removes any doubt as to the reality of the difference in 
the absorption of the hard component in air as compared 
to dense matter. These measurements give, then, strong 
support to the mesotron decay hypothesis. 

For the evaluation of the lifetime we use the following 
formula for the attenuation factor (reduction of intensity 
at the lower level): ; 


A=(E/(E+Bx)) exp [uc?/crB], 


where E is the energy at the lower level, B the specific 
energy loss in air (2X 10° ev/(g/cm?) = 2.4 108 ev/km air), 
x the distance traveled, u the rest mass of the mesotron, 
and r the lifetime at rest. This formula takes into account 
the change of energy, i.e., the change in the Lorentz time 
factor along the path. 

From the data of the table the lifetime comes out to be 
(uc?/10* ev) X10-* sec., where y=1.2, when A is calcu- 
lated from the values at 8 and 15 inches Pb (corresponding 
to an average mesotron energy E~4.4X 108 ev) and y=2.4 
for the values between 3 and 8 inches Pb (average energy 
~2.8X 108 ev). The latter value might already be falsified 
by a contribution from the soft component. The accuracy 
of the measurements is not yet sufficient to make certain 
the low value of 1.2 for y, but it seems unlikely that it 
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could be larger than 2. There would be no objection against 
a lifetime as small as 1.2 X 10~ sec. from other observations 
(inclination measurements, barometer effect), if the pro- 
duction of the mesotrons takes place at a somewhat lower 
altitude than that of the maximum of the soft component 
or if the mean energy of the mesotrons is somewhat larger 
than heretofore assumed. 

The value of the mesotron lifetime is then of the same 
order of magnitude as found in previous determinations, 
The actual measurements themselves are in reasonable 
agreement with the recent measurements of Rossi and his 
collaborators' insofar as a comparison can be made. From 
the fraction 26/34 of the mesotrons capable of traversing 
12.5 cm or more of lead and which survive after a distance 
of 1900 meters one calculates a mean path of 7 X 10° cm as 
compared to Rossi’s value of 8.5 10° cm. 

The observations are being continued, and a full report 


will be published later. 
W. M. NIELSEN 


C. M. RYERSON 
L. W. NorpHEI™ 


Duke University, 
Durham, North Carolina, 
K. Z. MorGAN 


December 19, 1939. 


1B. Rossi, H. van Norman Hilberry and J. Barton Hoag, Phys. Rev. 
56, 837 (1939). 


Alpha-Particle Groups from the Disintegration of 
Beryllium by Deuterons 


We have been investigating the energies of the alpha- 
particles from the reaction Be*(d,a)Li’, with the variable 
pressure absorption cell described by N. M. Smith, Jr.! 
Since the mean range of these alpha-particles is about 3 cm 
in air at the deuteron energies available to us, and the 
effective length of the absorption cell is 6.699 cm, the 
operating pressures ranged from 15 to 35 cm. 

Thin targets containing beryllium were prepared by 
volatilizing beryllium fluoride from a platinum wire held 
in the cooler parts of an oxygen-gas flame, and then 
inserting polished nickel target disks in the flame for 2 
to 3 seconds. The cell was calibrated with a thin polonium 
source, which could be interchanged with the bombarded 
targets. 

Figure 1 shows results with polonium, beryllium, and 
with a thin carbon target made by smoking one of the 
nickel disks in a gas flame. The scaling circuit was heavily 
biased so that counts were not recorded unless the alpha- 
particles passed through the 3-mm chamber at a velocity 
corresponding to their highest ionizing efficiency. The 
abscissae are pressures in the cell in centimeters of mercury, 
The ordinates for the carbon and beryllium curves are 
counts per 10-* coulomb. The counting was automatically 
stopped by an integrator which tripped at 796 micro- 
coulombs. The peak of the polonium curve corresponds to 
the counting of about 150 particles. The cell temperature 
for the polonium was 31.6°C, for the beryllium, 29°C. 
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Pressure in cell in cm of Hg 
Fic. 1. Counts from Po, Be and thin carbon targets. 


Table I shows the numerical results. Trial 2 of the 
table is that illustrated in Fig. 1. The energy release in this 
reaction may be calculated from the set of masses pre- 
viously established in this laboratory,’ and this calculated 
value is found to be 7.07 Mev, in excellent agreement with 
that obtained from the more energetic alpha-particle 
group. It is seen that there are two distinct energy groups, 
0.31 Mev apart, and of about equal intensity at this 
deuteron energy. The points on the line marked C+D 
were taken with a carbon target, because the alpha- 
particles from the reaction C%(d,@)B" might appear in 
this region if our targets, in spite of our precautions, 


were contaminated with carbon.? This reaction, although 


TABLE I. Results. Energies in Mev. 


ENERGY OF ENERGY OF 


DEUTERON SLOWER FASTER 


TRIAL ENERGY ALPHAS ALPHAS Q SE 
1 _ 0.351 4.34 4.64 7.05 0.30 
2 0.326 4.35 4.67 7.10 0.32 


observed by Cockcroft and Lewis* with 0.55-Mev deuterons 
apparently has too small a cross section to appear under 
our conditions. Visual observation of the oscillograph 
during the carbon counts showed that most of them were 
spurious. A similar negative result has been obtained from 
a target of potassium fluoride. 

These results indicate an excited state of Li’ at 0.48 


- Mev, which is probably the same as that appearing in 


the proton groups from disintegration of Li* by deuterons. 
This work was supported in part by a grant from the 
American Philosophical Society. 
SAMUEL K. ALLISON 
EvLizaBeTH R. GRAVES 


Lester S. SKAGGS 
Ryerson Physical Laboratory, 
University of Chicago, 
Chicago, Illinois, 
January 2, 1940. 


iW. M. oe. I. , Phys. Rev. 56, 548 (1939). 
Hatch hys. Rev. 54, 165 (1938). 


Cockcroft and W. B. Lewis, cz R 
aikeviewed by W. Maurer and J. B. Fisk, 


. Soc. 154, 261 (1936). 
its. f. Physik 112, 436 
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Internal Conversion in 6.4-Hour Cd 


Valley and McCreary! have recently measured the 
internal conversion electron spectrum of 6.4-hour Cd. 
They point out that the ratio of the number of L; to the 


‘number of K conversions is approximately 1 for the 92-kev 


gamma-ray, while the conversion of a gamma-ray of the 
same energy in Ga*’ gives the usual ratio of one to eight. 
The purpose here is to show that this result does not 
appear to be inconsistent with the theory. 

The ratio of one to eight is predicted theoretically for 
gamma-radiation of energy much greater than the ioniza- 
tion potential of the electron.? In the case of Cd, the 
ionization potential of the K shell is about thirty percent - 
of the gamma-ray energy (92 kev), therefore the above 
approximation is no longer valid. The effect is to reduce 
the probability of converting a K electron, as will be 
shown immediately. On the other hand, the L ionization 
potential is only three percent of the gamma-ray energy, 
so the L conversion is not reduced. The ratio of L to K 
conversions is therefore increased. Since the K-ionization 
potential of Ga is three times smaller than that for Cd, 
the reduction in the K-conversion probability for Ga 
should be much smaller. 

The effect of a large K-ionization potential (Jx) is 
twofold: 

(1) It reduces the kinetic energy (ex=hvy—Ix) of the 
ejected K electron and therefore reduces the density of 
energy states into which it can go, since the density is 
proportional to the square root of the kinetic energy. 
This reduction in the statistical weight of the transition 
decreases the transition probability. 

(2) It has been shown by Taylor and Mott? that for 
long wave-length gamma-radiation, internal conversion is 
almost entirely due to direct interaction of the K (or L) 
electrons with the nucleus. If the radiation is assumed to 
be dipole, the potential acting on the electron is of the form 
P/r, where P is the effective dipole moment of the nucleus. 
Therefore the matrix element for the electronic transition 
is Pf =Pf For dipole radiation, the 
final electronic wave function, y¥, is a P function, which 
has a zero of order r/A, A=h/(2me)!, at the nucleus. The 
matrix element is therefore proportional to (1/A) f° Yoerdr, 
and, since Yo and ¢(=Ay/r) are slowly varying functions 
of r that simply determine the effective limits of integra- 
tion, this is proportional to (1/A)So7rdr = x*/2d. The upper 
limit, x, is equal to the smaller of ao and 4/2, where ap is 
the radius of the K or L shell as the case may be. The 
transition probability is proportional to the square of this 
matrix element, i.e., proportional to x*/)*. 

The ratio of L to K conversion is given by 


wr/wr = = (4) 
where the factor } is the value of the probability for Az =A« 
(i.e., for hy>>Ix). For Cd) ~4.7, 
~4.4, \,~4 in units of cm. Since ax <Ax/2 
and az>dAz/2, the ratio becomes 
(wr /wr)ca = (1/128) (An*AL/ax*) = 2. 


On the other hand, for Ga, ax~1.7, a, ~6.8, 
A, in the same units, and 
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Since the estimates are extremely rough, it is not sur- 
prising that these values differ from the experimentally 
observed values! of 1 and } by a factor of two. However, 
the ratio of the ratios, (wz/wx)Ga/(wz/wx)ca, might be 


expected to be somewhat independent of the systematic. 


errors in the rough estimate, and it turns out to be approxi- 
mately }, in good agreement with the experimental 


value. 
R. G. Sacus 
The George Washington University, 


Washington, D. C., 
December 21, 1939. 
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On the Neutron’s Magnetic Moment 


As announced previously,' we have calculated in detail 
the influence of magnetization on the transmission of 
neutrons through polycrystalline iron. Since the coherent 
scattering is concentrated in Debye-Scherrer rings and 
since the transmission effect depends in a sensitive manner 
on the scattering angle, the final result differs radically 
from that previously calculated? for randomly oriented 
atoms, 

The only essential uncertainties in the theoretical 
evaluation occur from minute deviations from magnetic 
saturation and from our incomplete knowledge of the 
distribution of the magnetically active shell in the iron 
atom. In these calculations the neutron’s magnetic moment 
was assumed to be two nuclear Bohr magnetons. The 
depolarization effects due to incomplete saturation can 
only diminish the transmission effect. 

To determine the extension of the magnetically active 
shell we have used various density functions presented in 
the literature. Slater’s* hydrogenic functions give so small 
a form factor for the magnetic scattering that no effect 
could be expected. Adopting the analytic representation‘ 
for the Hartree distribution of the 3d electrons in Cu, to 
Fe, we obtained values for the transmission effects which 
are very much smaller than any observed previously and 
differ from those given by Bloch and Alvarez’ by a factor 
of approximately }. The discrepancy is considerably 
enlarged due to the presence of depolarization effects 
which manifested themselves experimentally in a deviation 
from the law: added transmission proportional to the 
square of the thickness. 

Attempts to explain the discrepancy can in our opinion 
only be based upon the two following possibilities: (1) The 
form factor used in these calculations is too small because 
either the charge distribution is too extended or due to 
some kind of spin orbit coupling the distribution of 
magnetic moment is contracted with respect to the dis- 
tribution of charge. (2) Since only the gyromagnetic ratio 
for the neutron has been determined experimentally,° it is 
still possible to assume that the neutron’s spin as well as 
its moment may be three times® as large. In this case the 
transmission effect observed would agree with the calcu- 
lated value if a reasonable correction for the depolarization 
effect is introduced. 
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We are conscious that a value for the neutron’s spin 
larger than } would seriously change all theoretical con- 
siderations concerning neutron-proton-scattering.’? We 
shall enter into this question more fully in a subsequent 
paper in which the theoretical aspects of the new hypothesis 
will be discussed. Here it may suffice to say that present 
experimental evidence® is far from being completely in 
favor of the commonly presented theoretical view. 

Observations on the paramagnetic scattering’ of slow 
neutrons allow the determination of the neutron’s magnetic 
moment directly with comparative case. Assuming a spin 
of 3 and using the gyromagnetic ratio as given by Bloch 
and Alvarez,* the scatte1ing cross section of free Mn** for 
a very long neutron wave-length should be approximately 
100 x 10-** cm?. ObServations on differential cross sections 
under small angles would clearly decide whether the 
discrepancy here discussed is due to cause (1) or (2) men- 
tioned above. 

O. HALPERN 
M. H. JouNson 


New York University, 
University Heights, New York, 
November 14, 1939. 
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On the Magnetic Moments of Light Nuclei 


That the discrepancies! between theoretical and experi- 
mental values of nuclear magnetic moments may be due 
to oversimplifying assumptions in the theory concerning 
the ground states of the various nuclei has often been 
pointed out. A detailed investigation of an added inter- 
action of the form 


r2 


in the case of the deuteron? led to consideration of this 
coupling as an, explanation. Its possible effect on the 


magnetic moments of several light nuclei has been sug- . 


gested® and the cases of Li and N"™ have recently been 
discussed.‘ 

V has nonvanishing matrix elements between states 
having the same total angular momentum and the same 
spatial symmetry (i.e., characterized by the same “‘parti- 
tion” quantum number). The sequence of observed mo- 
ments for the stable nuclei having an odd number of 
protons and the same number of neutrons (H?, Li*, BY, N") 
suggests an increasing proprtion of *D, in the ground state, 
as provided by this coupling term, and calculations made 
on the basis of the Hartree model qualitatively confirm 
this suggestion. Furthermore, the g=1.788 and probable 
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J=} values for B" can be understood in the same way. 
Spherically symmetric forces alone lead to a ground state® 
of J=}. There are, however, two nearby states with 
J=}3, and a rough calculation makes it seem reasonable 
that they interact sufficiently for one of them to become 
the ground state. In the language of L.S coupling the ground 
level is then a mixture of ?P 3/2 and *Ds3;2, and the correct 
magnetic moment is obtained theoretically if there is 
nearly one-half as much 2D as ?P, 

With the modification produced by the interaction V the 
Hartree model is then consistently correct in its predictions 
of spins for the light nuclei, giving the observed* J =} for 
C™® in contrast to the a-particle model.’ Quantitative 
agreement cannot be expected, and there is as a matter of 
fact no justification for detailed calculation with single 
particle wave functions, but the signs and rough magni- 
tudes of the theoretical magnetic moments are in accord 
with observation. One serious discrepancy has not been so 
explained. The experimental value of 4=3.25 for Li’ is 
larger than the theoretical value as given by the Hartree 
or any other simple model. The angle dependent coupling 
has no appreciable effect, since there is only one low 
energy term with J=}. Until this case is resolved one 
cannot be too optimistic about understanding observed 
magnetic moments’ in terms of the interactions hitherto 
investigated in connection with the problem. 

It is unfortunate that the Alvarez-Bloch result for the 
neutron moment is not sufficiently precise to constitute a 
criterion for the validity of these considerations, The 
coupling V, in accounting for the quadripole moment of 
the deuteron, depresses the value of uy, while their result 
gives no evidence of departure from pure additivity of 
proton and neutron moments in the deuteron. 


M. PHILLIPS 
Department of Physics, 
Brooklyn College, 
Brooklyn, New York, 
December 1, 1939. 
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The Neutron-Deuteron Scattering Cross Section* 


We have calculated the elastic-scattering cross section of 
neutrons by deuterons for energies so small (thermal 
neutrons) that only the s part of the incident wave is 
appreciably different from the corresponding part of a 
plane wave. 

In the coordinate system of the center of mass of the 


three particles, defined by the transformation, 


r=r,— }(f2+fs), 
R= 


(1 refers to the incident neutron and 3 to the proton) the 
wave equation is 


[—(h?/2M)(2V,2+ 39,2) + Vist Ves— E]¥s(r,p) =0. 
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The subscript S is the spin quantum number, and 
E=E,+E’, where Ep is the deuteron binding energy and 
E’ is the kinetic energy of the incident neutron in the 
center of mass system. For the interactions V;; we have 
taken the general form! 


Vig = —C1—g—g1—g2) Pig + 
exp [—ri7?/a?], 


in which Qi;=}(1+@;-@;) is the spin exchange operator 
(o;, 0; are the Pauli matrices of amplitude 1) and Pj; is 
the, coordinate exchange operator. In the present calcula- 
tion V» was taken equal to 72 mc* and a, the range equal to 
0.254/(Mm)!C cm. Because of the complexity of the 
problem we have neglected any polarization of the deu- 
teron by the neutron. 
The following expansions were used: 


WV s(r,p) = 
+xs(123) So" F(r,E)y(p,E)dE, 


P12Q12¥s(r,p) =Go(r) ¥o(p)xs(123) 
+xs(123)So°G(r,E) Wp, E)dE. 


¥o(p) is the ground state deuteron wave function, xs(123) 


. is the total spin wave function for a definite value S of the 


spin, and ¥(p,E£) is the excited deuteron wave function. 

On substituting these expansions in the wave equation 
and integrating over the deuteron coordinates, one ob- 
tains two equations which can be combined by subtraction 
to give 


3h 


In this equation the approximation 


Ws(r,p) — (PQ) 12¥s =[Fo(r) —Golr) xs 
— (PQ) Fo(r) —Go(r) 


(PQ) Pole) xs 
was used. 

Since polarization has been neglected, ¢(r) is correctly 
symmetrized to take account of exchange effects. The 
scattering cross section can therefore be found from that 
solution ¢(r) which has the asymptotic form 


The total scattering cross section is then given by 
o(0)dw=2xf | f(0)|* sin 


If we substitute the general expression for V,; in the 
equation for ¢(r), the sum over the spins can be easily 
performed, and the differential equation can be replaced 
by an integral equation with the correct asymptotic form.* 
This integral equation has to be solved numerically. 
The calculation is greatly simplified by the use of Gaussian 
potentials and a Gaussian function for Yo(p). In the present 
calculation we chose 


which gives the binding energy of the deuteron to within 
about 5 percent. 
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The presence of the Gaussian functions in the integrand 
insures a rapid enough convergence so that the integral 
equation can be replaced by a finite set of linear algebraic 
equations. It was found in this calculation that for r=13a 
the contributions to the integrals were practically zero. 

Since both doublet (d) scattering, S=4 and quartet (q), 
S= 4 scattering are present, we had to solve two sets of 
equations for the two values of the total spin. For thermal 
neutrons (k=0) and for the special choice of potential 
given by g=0.2, g:=g2=0,' we obtained 


¢a(13a) =1.0546; (13a) = 0.749 
and 
|fa(@) | = |a(13a)—1 | (13a) = 1.59 cm, 
| = —1| (13a) = 7.302 X cm. 


From this the scattering cross section was found to be 


4 
=4.57 


The recent experiments of Carroll and Dunning (to be 
published) give the value 5.7X10~ for this scattering 
cross section. If one takes molecular binding into account, 
this value has to be reduced by a factor which lies between 
1.5 and 2.3 At best, then, the theoretical value is about 20 
percent too large. This discrepancy may be accounted for 
by polarization effects. 

We should like to express our thanks to Mr. Jerome 
Rothstein of Columbia University for aiding in the 


numerical work. 
Lioyp Motz 
Columbia 
New York, New York, 
JULIAN SCHWINGER** 
The University of California, 


Berkeley, California, 
December 18, 1939. 


* Publication assisted by the Ernest Kempton Adams Fund for 
Physical Research of Columbia University. 

** National Research Fellow. 

1G. Breit and E. Feenbe oH Phys. Rev. 50, 850 (1936). 

2N. E. Mott and H. S. Massey, The Theor: “wy e Atomic Collisions 
(Oxford Univ. Press, New York, 1933), Chap, I 

3L. I. Schiff, Phys. Rev. 52, 149 (1937). 


Absorption Spectrum of Single-Crystal ZnS Phosphors 
_ A knowledge of the absorption spectrum of a phosphor is 
necessary if one is to understand the process by means of 
which radiation is emitted by the phosphor. In view of the 
experimental difficulties encountered in obtaining the 
absorption spectrum of a powder, it seems desirable to 
study the absorption spectrum of a ZnS phosphor in the 
form of a single crystal; especially so, since there seems to 
be some uncertainty as to the nature of the absorption 
process in activated ZnS." 

‘Several ZnS single crystals, approximately one millimeter 
square and a few tenths of a millimeter in thickness, were 
obtained through the kindness of Dr. R. P. Johnson and 
Mr. F. B. Quinlan of the General Electric Company. 
The crystals were made by firing in vacuum a precipitated 
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Fic. 1. Relative intensity of the light transmitted by a ZnS-Cu single- 
crystal phosphor vs. the wave-length, in angstrom units. 


ZnS with about 0.004 percent Cu. Each of three crystals 
was mounted behind a small hole in black paper and then 
placed in front of the slit of a small Littrow quartz spectro- 
graph. Light from a condensed iron spark was focused on 
the crystal by means of a quartz lens. The spectra were 
photographed on Eastman Spectroscopic Plates, Type 
IV-O. Three identical two-minute exposures were taken 
with the ZnS crystal in position. The crystal was then 


removed and six comparison spectra of 1, 2, 4, 8, 16, and 32 © 


seconds exposure were placed on the plate so as to bracket 
each of the original spectra. The intensities of the lines in 
the original spectra were determined visually by matching 
densities with one of the comparison spectra, whose 
exposure time was then taken as a measure of the intensity 
in question, 

The absorption curves obtained in this way for three of 
the ZnS single-crystal phosphors were quite similar; a 
typical curve, shown in Fig. 1, indicates that for this 
phosphor, absorption is partial for wave-lengths less than 
4250A and is complete for wave-lengths less than 3650A. 
This absorption curve is in good agreement with the curve 
obtained by Kitchener! from a ZnS-Cu powder, using the 
diffuse reflection method. Gisolf has reported? a secondary 
absorption peak on the long wave-length side of the 
absorption edge of ZnS phosphors; nothing of this nature 
is clearly shown by the accompanying curve, although 
there seems to be a slight asymmetry in the absorption 
edge which may be due to an unresolved secondary 
absorption peak. Nor was any evidence found of absorption 
bands in the region 4200—5000A, such as have been reported 
by Kroger? for ZnS-MnS mixed phosphors. Any such 
partial absorption which would reduce the intensity of the 
transmitted light by as much as one-fourth should have 
been detected. 

CHARLES J. BRASEFIELD 


U. S. Coast Guard Academy, 
New London, Connecticut, 
December 22, 1939. 
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